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Automated Longwall Guidance and Control System, DRL Line Item MA03, Phase I 
Report. 


CONTENTS 


Page 

Foreword i ■* 

Contents 1 1 1 

1. INTRODUCTION 1-1 

2. SUfflARY AND CONCLUSIONS 2-1 

2.1 Vertical Control System (VCS) ........... 2-1 

2.1.1 Baseline System Parameter Determination 2-6 

2.1.2 Baseline VCS Performance Characteristics 2-8 

2.2 VCS Performance with Optimum Filtering/Predictive 

Algorithms 2-11 

2.3 Use of Impact Penetrometei/Reflectometer for VCS 

Control ..... 2-14 

2.4 Use of Natural Radiation DID for VCS Control. . . . 2-15 

2.5 Slaving of Bottom Drum. 2-16 

2.6 Recommended VCS System 2-18 

2.7 Performance of Recommended System 2-19 

2.8 Roll Control System . 2-21 

2.8.1 Roll Control Loop Bandwdth Determination .... 2-22 

2.8.2 Roll Control Loop Perfomiance Characteristics . . 2-23 

2.9 Yaw Alignment System 2-24 

2.9.1 Angle Cart Performance Characteristics 2-26 

2.9.2 Directional Gyro Performance Characteristics. . . 2-28 

3. VERTICAL CONTROL SYSTEM 3-1 

3.1 Functional Description 3-1 

3.2 Detailed Mathematical Modiil of the VCS 3-3 

3.2.1 Overall System Description. 3-3 

3.2.2 Chassis Location Model 3-5 

3.2.3 Ranging Arm Geometry 3-11 

3.2.4 CID Sensors 3-19 

3.2.5 Actuator System 3-32 

3.2.6 VCS Control Lawi 3-42 

3.2.7 Mine Model 3-53 

3.3 Simulation Details. 3-70 

3.3.1 Initialization 3-70 

3.3.2 Cut Surface Determination 3-71 

3.3.3 Performance Eva'uation. ............. 3-72 

3.4 Baseline System Performance 3-74 

3.4.1 Philosophy and :>pecification of Baseline System . 3-74 

3.4.2 Performance of Baseline System. 3-74 

3.5 Filtering and Predictive Algorithms for VCS .... 3-119 

3.5.1 Philosophy Behind the Developed Algorithms. . . . 3-119 

3.5.2 Schemes Using P resent Cut Data. . - 3-119 


CONTENTS (Continued) 


Paae 

3.5,3 Schemes Using Present and Past Cut Data 3-138 

3.6 VCS Performance with Predictive Algorithms 3-151 

3.6.1 Predictive Algorithm Simulation 3-151 

3.6.2 Predictive Algorithm Performance with Simplified 

Simulation 3-152 

3.6.3 Predictive Algorithm Performance with VCS 

Simulation 3-154 

3.7 Control of Bottom Drum for Baseline System 3-175 

3.7.1 Drum Control Procedures .... 3-175 

3.7.2 Performance with Bottom Drum Locked . 3-177 

3.7.3 Performance with Bottom Drum Slaved 3-179 

3.7.4 Performance with Bottom Drum Delayed Slaving. . . 3-184 

3.7.5 Performance with Rear Drum Top Follower 3-184 

3.7.6 Performance with Slaving and Sensitized Pick. . . 3-189 

3.8 Conclusions and Recommendations 3-195 

3.8.1 Recommended Sensor Configurations 3-195 

3.8.2 Recommended VCS Systems for Various Mine Appli- 

cations .... 3-199 

4. YAW ALIGNMENT SYSTEM 4-1 

4.1 Functional Description. . «... 4-1 

4.2 Description of the Angle Cart Measurement System. . 4-6 

4.2.1 Basic Angle Cart Measurement System 4-8 

4.2.2 Derivation of Weighting and Bias Estimation 

System 4-11 

4.2.3 Derivation of Angle Cart Measurement Algorithm 

Assuming that a Number of End Conveyer Sections 
Cannot be Measured 4-16 

4.3 Description of the Monte Carlo Technique to Deter- 

mine the Effect of Angle Cart Measurement Errors. 4-16 

4.3.1 Basic Measurement System 4-19 

4.3.2 Performance of Weighting and Bias Estimation 

System 4-21 

4.3.3 Performance of Basic Measurement System When a 

Number of Conveyer Sections Cannot be Measured. 4-23 

4.4 Description of the Yaw Advancement System Simula- 

tion Using the Angle Cart . 4-26 

4.4.1 Overall System Performance Using the Angle Cart . 4-33 

4.4.2 Recommended Measurement Algorithm for the Angle 

Cart 4-40 

4.5 Description of the Directional Gyro Measurement 

System 4-47 

4.5.1 Mathematical Model of Directional Gyro Error. . . 4-47 

4.5.2 Derivation of the Directional Gyro Measurement 

Algorithm Without Averaging .......... 4-50 


CO NTENTS (Continued) 


Paai 

4.5.3 Derivation of the Dire<;tional Gyro Measurement 

Algorithm with Averaging . 4-52 

4.5.4 Derivation of the Directional Gyro Algorithm with 

Drifting Bias Compensation 4-54 

4.5.5 Derivation of the Directional Gyro Algorithm with 

Constant Bias Compensation 4-55 

4.6 Description of the Monte Carlo Technique to Deter- 

mine the Measurement Accuracy of the Directional 

Gyro 4-56 

4.6.1 Directional Gyro Without Averaging 4-57 

4.6.2 Directional Gyro with Averaging . 4-65 

4.6.3 Directional Gyro with Drifting Bias Compensation. 4-68 

4.6.4 Directional Gyro with Constant Bias Compensation. 4-72 

4.7 Description of the Yaw Advancement System Simula- 

tion Using the Directional Gyro 4-76 

4.7.1 Overall System Performance Using the Directional 

Gyro. 4-78 

4.7.2 Recommended Measurement Algorithm for Directional 

Gyro. 4-81 

4.8 Conclusions and Recommendations 4-84 

5. ROLL CONTROL SYSTEM 5-1 

5.1 Functional Description 5-1 

5.2 Active Roll Control System Implementation 5-1 

5.3 Detailed Mathematical Model for Roll Control 

System 5-3 

5.3.1 Inclinometer Sensor Model ..... 5-3 

5.3.2 Vibration Environment 5-6 

5.3.3 Actuator > Shearer Model 5-7 

5.4 Roll Loop Design 5-9 

5.4.1 Open Actuator Loop Design 5-9 

5.4.2 Closed Actuator Loop Design 5-13 

5.5 Roll System Performance with Cross-Axis Accelera- 

tion 5-15 

5.5.1 Performance with Open Actuator Loop 5-15 

5.5.2 Performance with Closed Actuator Loop . 5-28 

5.5.3 Effects of Reduced Flow Rate 5-28 

5.6 Conclusions and Recommendations . . . ....... 5-35 

6. PRELIMINARY DEFINITION OF CONTROLS AND DISPLAYS ... 6-1 

6.1 Control & Display Philosophy for Longwall ..... 6-1 

6.1.1 Control & Display Location 6-3 

6.1.2 Communications & Interlocks 6-3 

6.2 Control and Display Definition. .......... 6-4 


CONTENTS (Continued) 


Page 


6.2.1 Processor Accessed Displays 6-12 

6.2.2. Processor Accessed Cornnands . 6-12 

6.3 Preliminary Control and Display Concepts. ..... 6-14 

6.3.1 Headgate Control Station Panel Layout . 6-14 

6.3.2 Headgate Control Station Console Configuration. . 6-18 

6.3.3 Shearer Control Panel 6-22 

6.3.4 Tailgate Control Station 6-22 

6.3.5 Chock Control Panel ..... 6-22 

6.3.6 Control and Display Component Selection ..... 6-24 

6.4 Proposed Phase II Trade Studies 6-24 

7. OVERALL L0N6WALL SYSTEM REQUIREMENTS 7-1 

7.1 Vertical Control System 7-1 

7.2 Roll Control Loop 7-4 

7.3 Yaw Alignment System 7-4 

7.4 Controls and Displays 7-5 


Figure 


3-1 VCS Simulation Model Block Diagram 

3-2 Mine and Shearer Coordinate System and Parameter Def 

inition 

3-3 Shearer Velocity Profile 

3-4 Shearer Chassis Geometry 

3-5 Ranging Arm Diagram 

3-6 Ranging Arm Elevation Geometry 

3-7 CID Sensor Geometry for Top Drum . 

3-8 Last Cut Follower Geometry for Top Drum 

3-9 CID Sensor Geometry for Bottom Drum 

3-10 Last Cut Follower Geometry for Bottom Drum. .... 

3-11 Nucleonic CID Calibration Curve 

3-12 Natural Radiaticn CID Calibration Curve 

3-13 Sensitized Pick Geometry 

3-14 Penetrometer - Ref 1 ectometer Sensor 

3-15 VCS Hydraulic System 

3-16 Actuator and Ranging Arm Simulation Model ..... 

3-17 Flow Split Gains 

3-18 Open Loop Actuator Displacement Response. ..... 

3-19 Open Loop Ranging Arm Rate Response 

3-20 Open Loop Drum Height Response 

3-21 Closed Loop Actuator System ...... 

3-22 Closed Loop Actuator Displacement Response 

3-23 Closed Loop Ranging Arm Rate Response 

3-24 Closed Loop Drum Height Response. 

3-25 CID Control Law 


3-4 


3-6 

3-8 

3-9 

3-13 

3-15 

3-18 

3-20 

3-21 

3-22 

3-23 

3-27 

3-29 

3-31 

3-33 

3-34 

3-36 

3-39 

3-40 

3-41 

3-43 

3-44 

3-45 

3-46 

3-47 


reproducibility 
original page is 


OF THE 
POOR 


3-26 

3-27 

3-28 

3-29 

3-30 

3-31 

3-32 

3-33 

3-34 

3-35 

3-36 

3-37 

3-38 

3-39 

3-40 

3-41 

3-42 

3-43 

3-44 

3-45 

3-46 

3-47 

3-48 

• 3-49 

3-60 

3-51 

3-52 

3-53 

3-54 

3-55 

3-56 

3-57 

3-58 


CONTENTS (Continued) 


Paai 


CID and LCF Control Law 3-48 

Pick Control Law 3-49 

CID and Pick Control Law 3-50 

CID, Pick and LCF Control Law 3-51 

Autocorrelation of York Canyon Mine 3-54 

Autocorrelation of Robinson Run Mine, 3-56 

Power Spectral Density of York Canyon Mine 3-58 

Power Spectral Density of Robinson Run Mine . . . . . 3-60 

Cross Correlation Between Top and Bottom of York 

Canyon Mine .... 3-63 

Cross Correlation Between Top and Bottom of Robinson 

Run Mine. 3-64 

Generated Profiles for the Robinson Run Top 3-67 

Generated Profiles for the Robinson Run Bottom. . . . 3-68 

3D View of the Generated Bottom for the Robinson Run 

Mine. 3-69 

Drum Cut Geometry 3-73 

VCS Performance with 0.75 Feet CID Displacement ... 3-76 

VCS Performance with 1.5 Feet CID Displacement. . . . 3-77 

VCS Performance with 2.5 Feet CID Displacement. . . . 3-78 

VCS Performance with 1 Sec. CID Output Interval ... 3-81 

VCS Performance with 5 Sec. CID Output Interval . . . 3-82 

VCS Performance with 10 Sec. CID Output Interval. . . 3-33 

Performance with Nucleonic Sensor at 0.25 Sec. Output 

Interval with Noise . 3-86 

Performance with Nucleonic Sensor at 1.0 Sec. Output 

Interval with Noise 3-87 

Performance with Nucleonic Sensor at 2.0 Sec. Output 

Interval with Noise 3-88 

Performance with 3 In. Natural Radiation Sensor at 

a 2 Sec. Output Interval with Noise 3-89 

Performance with 5 In. Natural Radiation Sensor at 

a 1 Sec. Output Interval with Noise 3-91 

VCS Performance in the York Canyon Mine 3-94 

VCS Performance in the Robinson Run Mine 3-95 

VCS Performance with a 6 In. Bias . . . 3-98 

VCS Performance with CID Control Law Leaving a 2 In. 

Bias. 3-100 

VCS Performance with CID & LCF Control Law Leaving a 

2 In. Bias 3-101 

VCS Performance with CID & Pick Control Law Leaving 

a 2 In. Bias. 3-102 

VCS Performance with CD, Pick & LCF Control Law 

Leaving a 2 In. Bias 3-103 

VCS Performance with CID, Penetrometer & LCF Control 
Law Leaving a 2 In. Bias. . 3-104 


CONTENTS (Continued) 


3-59 

3-60 

3-61 

3-62 

3-63 

3-64 

3-65 

3-66 

3-67 

3-68 

3-69 

3-70 

3-71 

3-72 

3-73 

3-74 

3-75 

3-76 

3-77 

3-78 

3-79 

3-80 

3-81 

3-82 

3-83 


VCS Performance for CID, Pick Control Law with 90% 

Pick Errors for a 2 In. Bias 

VCS Performance for CID, Pick Control Law with 80% 

Pick Errors for a 2 In, Bias 

VCS Performance for CID, Pick and LCF Control Law 

with 90% Pick Error for a 2 In. Bias 

VCS Performance for CID Pick and LCF Control Law with 

80% Pick Errors for a 2 In. Bias 

VCS Performance with 100% Correct Pick Control. . . . 
VCS Performance with 90% Correct Pick Control .... 
VCS Performance with 100% Correct Penetrometer-Re- 

flectometer, LCF Control 

VCS Performance with 90% Correct Penetrometer-Reflec- 

tometer, LCF Control 

Block Diagram of Scheme 1 . . 

Scheme III 

Scheme IV 

VCS Performance with Scheme 1 Estimation Using Total 

Measurements and Flat Bottom 

VCS Performance with Scheme 1 Prediction Using Total 

Measurements and Flat Bottom 

VCS Performance with Scheme 2 Estimation Using Total 

Measurements and Flat Bottom. .... 

VCS Performance with Scheme 2 Prediction Using Total 

Measurements and Flat Bottom. . 

VCS Performance with Scheme 1 Estimation, with CID 

Measurements and Bottom Undulations ... 

VCS Performance with Scheme 1 Prediction, with CID 
Measurements and Bottom Undulations ........ 

VCS Performance with Scheme 2 Estimation, with CID 
Measurements and Bottom Undulations ........ 

VCS Performance with Scheme 2 Prediction, with CID 

Measurements and Bottom Undulations . 

VCS Performance with Scheme 1 Estimation with Total 

Measurements and Bottom Undulations 

VCS Performance with Scheme 1 Prediction with Total 

Measurements and Bottom Undulations 

VCS Performance with Scheme 2 Estimation with Total 

Measurements and Bottom Undulations 

VCS Performance with Scheme 2 Prediction with Total 

Measurements and Bottom Undulations 

VCS Performance with Scheme 1 Estimation Under Actual 
Conditions. .................... 

VCS Performance with Scheme 1 Prediction Under Actual 
Conditions 


Pafli 

3-107 

3-108 

3-109 

3-no 

3-112 

3-113 

3-114 

3-115 

3-122 

3-140 

3-144 

3-157 

3-158 

3-159 

3-160 

3-162 

3-163 

3-164 

3-165 

3-166 

3-167 

3-168 

3-169 

3-171 

3-172 


CONTENTS (Continued) 


Page_ 

3“84 VCS Performance with Scheme 2 Estimation Under Actual 

Conditions 3-173 

3-85 VCS Performance with Scheme 2 Prediction Under Actual 

Conditions. ...... 3-174 

3-86 Slaving Sensor Geometry 3-176 

3-87 VCS Performance in Robinson Run Mine with Flat Bottom 

and Rear Drum Locked. ..... 3-180 

3-88 VCS Performance in York Canyon Mine with Flat Bottom 

and Rear Drum Locked. 3-181 

3-89 VCS Performance in Robinson Run Mine with Undulating 

Bottom and Drum Locked 3-182 

3-90 VCS Performance ln York Canyon Mine with Undulating 

Bottom and Drum Locked 3-183 

3-91 VCS Performance in Robinson Run Mine with Drum Slav- 
ing , 3-185 

3-92 VCS Performance in York Canyon Mine with Drum Slav- 
ing 3-186 

3-93 VCS Performance in Robinson Run Mine with Delayed 

Slaving and Floor Undulations 3-187 

3-94 VCS Performance in York Canyon Mine with Delayed 

Slaving and Floor Undulations 3-188 

3-95 VCS Performance with Top Follower in Robinson Run 

Mine 3-190 

3-96 VCS Performance with Top Follower in York Canyon 

Mine 3-192 

3-97 VCS Drum Slaving Performance with Front LCF 3-196 

3-98 VCS Drum Slaving Performance with Middle LCF. .... 3-197 

3-99 VCS Performance in Robinson Run Mine with Recommended 
CID, Pick and LCF Configuration Leaving a 4 In. 

Bias. 3-201 

3- 100 VCS Performance in York Canyon Mine with Recommended 

CIDj Pick and LCF Configuration Leaving a 4 In. 

Bias 3-206 

4- 1 Conveyor Advance-Sequence of Operations 4 4-2 

4-2 Conveyor Advance-Sequence of Operations 3 ..... . 4-3 

4-3 Conveyor Advance-Sequence of Operations 2 . ... . . 4-4 

4-4 Conveyor Advance-Sequence of Operations 1 4-5 

4-5 Angle Cart Measurement Device 4-7 

4-6 Conveyor Geometry 4-10 

4-7 Monte Carlo Simulation. ............... 4-17 

4-8 Basic Angle Cart System - Monte Carlo Results .... 4-20 

4-9 Angle Cart Systems - Monte Carl Results 4-22 

4-10 Basic Angle Cart System - Monte Carlo Results .... 4-24 

4-11 Basic Angle Cart System - Monte Carlo Results .... 4-25 


CONTENTS (Continued) 


4-12 

4-13 

4-14 

4-15 

4-16 

4-17 

4-18 

4-19 

4-20 

4-21 

4-22 

4-23 

4-24 

4-25 

4-26 

4-27 

4-28 

4-29 

4-30 

4-31 

4-32 

4-33 

4-34 

4-35 

4 - 36 

5 - 1 
5-2 


Yaw Advancement Simulation. ... 

Geometry for Defining Placement Commands. ..... 
Basic Angle Cart System - Yaw Advancement Simulation 

Results 

Basic Angle Cart System - Yaw Advancement Simulation 

Results .......... 

Basic Angle Cart System - Ya^ -Advancement Simulation 

Results . 

Basic Angle Cart System - Yaw Advancement Simulation 

Results . 

Weighting and Bias Estimation System - Yaw Advance- 
ment Simulation Results 

Weighting and Bias Estimation System - Yaw Advance- 
ment Simulation Results ..... 

Yaw Profiles - Basic Angle Cart System 

Efficiency of Basic Angle Cart System 

Efficiency of Basic Angle Cart System 

Directional Gyro NDlse Model 

Directional Gyro - Single Measurement - Monte Carlo 

Results 

Directional Gyro - Single Measurement - Monte Carlo 

Results ....... 

Directional Gyro - Single Measurement - Monte Carlo 

Results . ....... 

Directional Gyro - Single Measurement - Monte Carlo 


Results 

Directional Gyro - Single Measurement - Monte Carlo 

Results . . 

Directional Gyro - Averaging - Monte Carlo Results. . 
Directional Gyro - Drifting Bias Compensation - Monte 

Carlo Results . , 

Directional Gyro - Constant Bias Compensation - Monte 

Carlo Results 

Directional Gyro - Constant Bias Compensation - Monte 


Carlo Results , 

Directional Gyro - Single Measurement - Yaw Advance- 
ment Simulation Results 

Directional Gyro - Constant Bias Compensation - Yaw 


Advancement Simulation Results 

Efficiency of the Directional Gyro - Single Measure- 
ment System 

Efficiency of the Directional Gyro - Single Measure- 
ment System ........ . . 

Active Roll System Control Loop Implementations . . . 
Inclinometer Model Block Diagram. .... 


Page 

4-27 

4-34 

4-36 

4-37 

4-38 

4-39 

4-41 

4-42 

4-44 

4-45 

4-46 

4-48 

4-60 

4-61 

4-62 

4-63 

4-64 

4-69 

4-71 

4-73 

4-74 

4-79 

4-80 

4-82 

4 - 83 

5 - 2 
5-5 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


CONTENTS (Continued) 


Paae 

5-3 Shearer Geometry for Roll Control ..... 5-8 

5-4 Open Actuator Loop Performance 5-10 

5-5 Open Actuator Loop Performance 5-11 

5-6 Open Actuator Loop Performance 5-12 

5-7 Open Actuator Loop Time History with No Noise .... 5-14 

5-8 Closed Actuator Loop Performance 5-16 

5-9 Closed Actuator Loop Perfonnance 5-17 

5-10 Closed Actuator Loop Performance. . 5-18 

5-11 Closed Actuator Loop Time History with No Noise . . . 5-19 

5-12 Open Actuator Loop Response with 0,24 Hz Bandwidth. . 5-22 

5-13 Open Actuator Loop Time History with 2.4 Hz Inclinom- 
eter Bandwidth, 0.05 g Cross-Axis Acceleration and 

0.05 Hz Input 5-23 

5-14 Open Actuator Loop Time History with 0.24 Hz Inclinom- 
eter Bandwidth, 0.05 g Cross-Axis Acceleration and 

0.05 Hz Input 5-24 

5-15 Open Actuator Loop Time History with 2.4 Hz Inclinom- 
eter Bandwidth, 0.1 g Cross-Axis Acceleration and 

0.025 Hz Input 5-26 

5-16 Open Actuator Loop Time History with 0.24 Hz Inclinom- 
eter Bandwidth, 0.1 g Cross-Axis Acceleration and 

0,025 Hz Input 5-27 

5-17 Closed Actuator Time History with 2.4 Hz Inclinometer 
Bandwidth, 0.05 g Cross-Axis Acceleration and 0.05 

Hz Input 5-29 

5-18 Closed Actuator Time History with 0.24 Hz Inclinometer 
Bandwidth, 0.05 g Cross-Axis Acceleration and 0.05 

Hz Input. .. 5-30 

5-19 Closed Actuator Loop Time History with 2.4 Hz Incli- 
nometer Bandwidth, 0.1 g Cross-Axis Acceleration 

and 0.025 Hz Input. 5-31 

5-20 Closed Actuator Loop Time History with 0.24 Hz Incli- 
nometer Bandwidth, 0.1 g Cross-Axis Acceleration 

and 0.025 Hz Input. 5-32 

5-21 Open Actuator Loop Performance 5-33 

5-22 Closed Actuator Loop Performance 5-34 

5-23 Open Actuator Loop Time History with 2.4 Hz Inclinom- 
eter Bandwidth, 0.1 g Cross-Axis Acceleration, 4 

GPM Flow Rate 5-37 

5-24 Open Actuator Loop Time History with 0.24 Hz Incli- 
nometer Bandwidth, 0.1 g Cross-Axis Acceleration, 

4 GPM Flow Rate 5-38 

5-25 Closed Actuator Loop Time History with 2.4 Hz Incli- 
nometer Bandwidth, 0.1 g Cross-Axis Acceleration, 

4 GPM Flow Rate 5-39 


CON TENTS (Concluded) 

Pafli 

5- 26 Closed Actuator Loop Time History with 0.24 Hz Incli- 

nometer Bandwidth, 0.1 g Cross-Axis Acceleration, 

4 GPM Flow Rate . 5-40 

6- 1 Joy LW300 Shearer Control Panel ... 6-2 

6-2 Preliminary Headgate Control Panel 6-15 

6-3 Headgate Control Panel with Alpha-Numeric Keyboard. . 6-19 

6-4 Portable Panel Concept. 6-20 

6-5 Mobile Panel Concept. . 6-21 

6-6 Shearer Control Panel - Preliminary Concept ..... 6-23 

Table 


2- 1 VCS Performance SunKtiary 2-20 

3- 1 VCS Parameter Values. 3-12 

3-2 Mine Statistics 3-62 

3-3 Effects of CIU Location 3-75 

3-4 Effects of CID Output Interval ............ 3-80 

3-5 Effects of CID Noise . 3-84 

3-6 Comparison of Performance on Different Mines 3-93 

3-7 Effects of Nonlinear CID Calibration Curve 3-96 

3-8 VCS Control Law Comparison 3-99 

3-9 Effects of Sensitized Pick Errors ..... 3-105 

3-10 VCS Performance of Rock Presence Sensors Only .... 3-111 

3-11 Effect of Actuator Dynamics on VCS Performance. ... 3-117 

3-12 VCS Performance When Varying Cutting Speeds 3-118 

3-13 Comparison of Prediction Error Variances 3-137 

3-14 Predictive Algorithms - Effects of Parameter Varia- 
tions Simplified Simulation 3-153 

3-15 Prediction Algorithms - Effects of Modeling Errors 

VCS Simulation 3-156 

3-16 VCS Performance with Rear Drum Slaving 3-178 

3-17 VCS Performance with Drum Slaving and Sensitized 

Pick . 3-194 

3-18 Recommended VCS Sensor Configuration 3-198 

3-19 Results with Recommended VCS Sensor Configuration . . 3-200 

3-20 Recommended Sensor Configuration for Various Applica- 
tions 3-211 

5-1 Perfomance with Cross-Axis Disturbances 0.05 Hz 

Input and 8 GPM Flow Rate 5-20 

5-2 Performance with Cross-Axis Disturbances, 0.025 Hz 

Input and 8 GPM Flow Rate . 5-25 

5- 3 Performance with Cross-Axis Disturbances, 0.025 Hz 

Input and 4 GPM Flow Rate . 5-36 

6- 1 C & D Functional Requirements 6-5 

6-2 Preliminary List of Processor Accessed Displays . . . 6-13 


REPRODlTCn^ILITY OF THE 
ORIGIiYAL PAGE IS POOR 


INTRODUCTION 


Dwindling domestic and world oil and gas supplies with the inevitable 
rise in the price of these fuels have generated deep concern within govern- 
ment, Industry, and the general public as to how the energy needs of 
the country can reasonably and economically be met in the future. The 
increased use of coal, our most abundant fossil fuel reserve, during 
the next quarter century will necessitate increasing the efficiency of 
underground coal extraction. Longwall mining techniques have tne potential 
of greatly increasing the coal yield per acre and coal production per 
man pe*' shift since it is essentially a continuous mining process. In 
addition, since longwall is a continuous mining process employing continuous 
haulage, It is extremely well suited for automation which is the subject 
of the present study. Automating longwall coal extraction will not only 
Increase production but also minimize the artiount of foreign material 
taken along with the coal thus reducing sorting time and cutter bit wear. 

In additions automating or remoting the longwall mining process will 
increase operator health and safety by removing the miner from the shearer 
and thus the hazards encountered in the immediate cutting area. 

The present study has been divided into two phases. Phase I 
was primarily concerned with the analyses and simulation of candidate 
Vertical Control Systems (VCS) and Face Advancement Systems (FAS) re- 
quired to satisf'sHorily automate the longwall system. The purpose of 
these studies were to specify the desired overall longwall system con- 
figuration for preliminary design which will be performed during Phase 
II of the study. This report outlines the analyses, and simulations 
that were performed during Phase I of the subject study which has led 
to the satisfactory specification of the overall longwall system on which 
a preliminary design will be performed during Phase II. 
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2 . 


SUMMARY AND CONCLUSIONS 


In the paragraphs that follow, a summary of the results obtained 
and the conclusions drawn during the Phase I portion of the study will 
be outlined by major longwall subsystem. These major subsystems are 
the Vertical Control System (VCS), and the Face Advance System (FAS) 
which consists of the yaw and roll control subsystems. 


In a study of this nature, a procedure 1s needed by which the perform' 
ance of various control system configurations could be compared in order 
to determine which one is "best" for a particular application. In order 
to accomplish this end, "baseline" system configurations for all of the 
control subsystems required for longwall automation were defined. The 
baseline definitions used the work performed by MSFC as a starting point 
and were the "simplest" (and hence least expensive and probably most 
reliable) subsystem configuration that would satisfactorily perform a 
particular control function. It was against the performance achieved 
by the "baseline" systems that various other control system configura- 
tions/additions/modifications were compared. 

. Proposed modifications/additions to the baseline system definition 
were evaluated by considering the improved performance realized vs. the 
additional complexity and relative cost introduced by incorporation of 
the modification, and only those modifications/additions to the baseline 
systems that were cost effective were incorporated. Following this pro- 
cedure, a highly reliable and economically advantageous longwall system 
design emerged, 

2.1 Vertical Control System (VCS) 


The function of the VCS is to control the two cutting drums 
of the double ended ranging arm shearer so that they remain within the 
coal seam. The VCS should be capable of taking all the coal on the roof 
and floor, or if desired, leave a given amount (i.e. , bias) of head and/or 
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floor coal while taking a minimum amount of rock. In addition, the dis- 
tance between two successive cuts would have to remain within given limits 
in order that the roof supports advance satisfactorily. 

In order to analyze and design the VCS, a typical longwall 
shearer had to be chosen for computer simulation, The longwall shearer 
chosen was the Joy LW 300 and all control system studies Involving the 
shearer were performed with this model. In addition, a reasonable repre- 
sentation of the coal/shale interface that would be encountered during 
1n-m1ne operations had to be simulated. It became apparent early In 
the study that using artificial “worst case" profiles would not be very 
satisfactory from two viewpoints: 

a. It Is very difficult if not impossible to define a 
"worst case" coal shale interface without having access 
to the actual interface profiles in most or possibly all 

of the operational longwall faces and those being projected 
for longwall mining. 

b. Even if a "worst case" profile could be defined, it 
would place extreme requirements on the required longwall 
system resulting in an unrealistic and costly system design. 

It was, therefore, decided to model typical coal/shale interface 
profiles that one would commonly encounter during longwall mining operations. 
To accomplish this end, use was made of a survey performed by A. D. Little, 
Inc., the results of which were given in a report titled, "Last Cut Survey," 
published in June 1977, in which the coal/shale interface profiles for 
six mines representing a variety of geologic conditions and seam heights 
were given. Statistical analyses were performed on the twelve profiles 
(i.e., top and bottom) given and it was determined that all of the mine 
profiles were statistically similar. This result certainly gives credence 
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to the hypothesis that coal/shale interface profiles are more similar 
than different when going from mine to mine. Therefore, the York Canyon 
profiles were chosen from the given data set as representing typical 
(I.e., mean) coal/shale interface profiles that would be encountered 
during in-mine operations and most performance runs were made with these 
profiles. Performance evaluations made with other mine profiles had 
little effect on overall performance as will be discussed In Section 3. 

In order to evaluate VCS performance and compare various con- 
figurations with respect to each other, a measure of performance had 
to be defined. There are numerous performance measures that one can 
define for the VCS, no single one of which will be an adequate measure 
of performance for. all operational conditions. The most commonly used 
to characterize the performance of a control system is the RMS of the 
error between the desired and actual control achieved which, in our case, 
is the RMS of the difference between the desired and actual cuts achieved, 
or the cut error. Although this is generally an excellent measure of 
VCS performance, it does not totally specify all aspects of VCS per- 
formance. Other performance characteristics that are of interest are 
the amount of time spent cutting rock, the volume of rock taken, the 
volume of coal erroneously left, and the volume of coal erroneously taken. 
Therefore, in order to augment the RMS cut error, a number of other per- 
formance parameters were also determined which are: a) percent of time 

spent in rock; b) volume of rock taken; c) volume of excess coal taken; 
d) volume of coal erroneously left; and e) the average of the absolute 
value of the cut error. These performance parameters, in conjunction 
with the RMS cut error completely specify VCS performance over its desig- 
nated range of operation. 

A block diagram of the VCS computer simulation is shown in 
Figure 3-1. The major features of the simulation are listed below: 
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a. Detailed modeling of the mine geometry including actual 
Interface profiles and techniques for generating statisti- 
cally similar coal/shale interface profiles from a given 
initial profile which is used to give realistic results 
when multiple passes are made. 

b. Detailed geometrical model of the shearer including 
its location along, and orientation with respect to, the 
face. 

c. Technique for determining the actual track that would 
result when laying five foot conveyor sections along a 
roughly cut bottom. 

d. Detailed non-linear representation of shearer actuator 
dynamics including the hydraulic coupling between the two 
ranging arms. 

e. Detailed representation of the nucleonic and natural 
radiation coal interfat* detectors including the nonlinear 
calibration curves, Poisson distributed sensor noise, and 
air gaps. 

f. Complete representation of coal/shale presence sensors 
including a probabilistic representation of their correct- 
ness. 

g. Complete representation of last and present cut followers. 

h. Modular hierarchical control law algorithm that is able 
to accommodate a variety of sensor complements and control 
law implementations. 
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1. The effect of drum filtering on the cut obtained. 


For a more detailed discussion of the VCS simulation see sections 
3.2 and 3.3 

The "baseline” system configuration for the VCS consisted 
of the following: 


a. An active nucleonic coal interface detector (CID) for 
measuring coal thickness on the roof. 

b. Sensitized picks incorporated into each of the two 
cutting drums in order to detect the presence of coal or 
rock . 

c. A last cut follower mechanism in order to measure the 
present cut relative to the last one made. 

d. The bottom cutting drum is slaved to the top one 

if a given amount of floor coal is to be left. The reason 
for not using a CID for proportional floor cutting drum 
control is that it is not feasible to mount such a detector 
on the floor given all the debris and water found there. 

e. Use presently obtained CID data in a measure and cut 
algorithm as opposed to optimal filtering/prediction algo- 
rithms which can be used to compensate for CID noise and 
that the measurement point is some 2.5 feet behind the 
actual drum cutting point. 


f. Use of the nonlinear hydraulic actuation system commonly 
found on longwall shearers. 


EEPRODUCIBITJTY of the 
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Modifications/additloriS to the defined VCS "baseline" con- 
sidered are listed below. 

1. Use of optimal filtering/predictive control algorithms 
for cutting drum control. 

2. Use of a natural radiation CID. 

3. Use of an impact penetrometer/ref lactometer for a coal/ 
shale presence sensor. 

A number of techniques for bottom drum slaving were considered. 

They are: 


1. Maintaining constant height between both cutting drums. 

2. Using CID measurements instantaneously and in a delayed 
mode. 

3. Use of the last cut follower. 

4. Use of a last cut follower located in the middle of 
the shearer. 

5. Use of a separate present cut follower located above 
the rear drum. 

Figure 3-86 shows the shearer slaving geometry and the various 
slaving techniques listed above. 

Baseline System Parameter Determination - For the 


baseline VCS defined, system parameters (i.e., gains, sensor locations, 
sampling and averaging times, etc.) had to be determined in order to 


obtain optimum performance for this configuration. One of the first 
parameters to be determined was the loop gain/compensation to be used 
for the position loop closed around the actuators which control cutting 
drum height. Various loop gains were investigated and it was determined 
that a gain of 170 (see Figure 3-21) yielded best overall performance 
even though the actuator did have a limit cycle at this gain value. 

The only problem that can occur if the actuator is operated in the limit 
cycle region is that it may heat excessively. Should that be the case, 
a heat exchanger could easily be designed which will maintain tempera- 
tures within acceptable limits. Therefore, a loop gain of 170 was speci- 
fied and used for the gain of the actuator loop in all of the control 
studies subsequently performed. 

Another parameter to be determined was the mounting location 
of CIO sensors. It was clearly shown (see Table 3-3) that system performance 
improves aS the CIO is moved closer to the cutting drum (and hence the 
cutting point) since this minimizes the distance between the point measured 
and the point being controlled. The closer the CID can be moved to the 
cutting point the better is the probability that the CID measurement 
reflects the coal/shale interface at the cutting locatiooj hence improv- 
ing the measurement. Therefore, the CIO's were located 2.5 feet behind 
the cutting drum which was as close to the cutting drum as the CID 
could be mounted without interference. 

Another critical parameter to be determined was the CID 
output interval or averaging time. Again, it was clearly shown (see 
Table 3-4) that in the absence of CID measurement noise, increasing the 
CID output interval degraded system performance. The reason for this 
is that as the output interval is increased, greater measurement lags 
are introduced into the system degrading its performance. However, the 
longer the output interval or averaging time, the smaller the effect 
of CID measurement noise and hence the more accurate the measurements 
obtained. Therefore, there should be an output interval which will 


balance the detrimental effect of increased output interval due to added 
lag against the benefit gained by reducing the effect on CIO measurement 
noise, thus yielding optimum system performance. Such an output interval 
was indeed found (see Table 3-5), and was 0.25 second for the active nucleonic 
CIO. When a natural radiation CIO was employed the optimum output inter- 
vals were 1 and 2 seconds for 5 and 3 inch detection crystals respectively. 

2.1.2 Baseline VCS Performance Characteristics 


2. 1,2.1 Effect of Nonlinear CID Calibration Curve 


The active nucleonic or passive natural radiation 
CID‘s have calibration curves that are nonlinear (see Figures 3-11 and 
3-12) and saturate as the thickness of coal increases. Therefore, for 
the same CID measurement noise, the error in the sensor measurement 
would increase as the desired thickness of coal to be left approaches 
the saturation characteristics of the sensor thus giving degraded system 
performance. Such an effect has been found and is outlined in Table 
3-7 for the active nucleonic CID. Examination of this table indicates 
that for a one second output interval the RMS cutting error went from 
1.61 inch to 2.39 inch, an increases of 48.45 percent, when the coal 
thickness to be left increased from 4 to 6 inches. However, it should 
be noted that cut error will remain invariant as the desired coal to 
be left is 5 inches or less since in this range the CID calibration curve 
is essentially linear. 

2. 1.2.2 Effect of Adding Sensitized Pick and Last Cut 
Follower (LCF) to CID 

The effect of adding a sensitized pick and using 
it in conjunction with the CID in a hierarchical manner is to improve 
cutting performance. The technique that is used is to command the cutting 
drums on the basis of CIO information and when the pick indicates that 
rock is being mined to override the CID coimiand and depress the cutting 
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head a given amount. The improvement realized by the addition of the 
sensitized pick is most pronounced when attempting to leave a small coal 
thickness (i.e,, between zero and 2 inches), since it is during this 
mode of operation that using CID information only has a high probability 
of coninanding the drum to cut the roof due to CID measurement errors. 

Once the drum has cut into the roof the CID cannot say how far into the 
roof the drum has cut and gives an indication of zero coal thickness 
regardless of the depth the cutting drum has penetrated the rock. 
Therefore, once the roof is penetrated the cutting drum would only be 
comnanded to depress the coal bias amount even if the roof were pene- 
trated considerably more than this value, hence degrading cutting per- 
formance. The addition of a sensitized pick does not allow the cutting 
drum to appreciably penetrate the roof, thereby keeping the CID in its 
linear range, thus improving cutting performance. The addition of a 
last cut follower to the CIO and sensitized pick and employing this 
sensor configuration in a hierarchical manner where the last cut follower 
is given highest priority, the sensitized pick next, with the CID having 
lowest priority, further improves cutting performance. The reason for 
the improved performance when introducing a last cut follower is that 
it filters even further errors in the CIO measurements since the present 
cut cannot deviate more than a prescribed amount (i.e., + 2 inches) from 
the previous cut. 


It should also be noted that the cutting accuracy achieved when 
using a sensitized pick in conjunction with the CID or CID and last cut 
follower, is not appreciably affected as the accuracy of the pick in- 
formation degrades to 80 percent. The reason for this is that the pick 
is one of two or three sensors controlling the VCS and hence a wrong 
piece of information from it does not have a one-to-one effect on system 
performance. The results discussed above are outlined in Tables 3-8 
and 3-9. 




2. 1.2.3 VCS Performance with Sensitized Pick Only 


When attempting to take all of the coal the VCS 
will be controlled by the sensitized pick only, or possibly acting in 
conjunction with a last cut follower. In this mode of operation error 
In the sensitized pick information does have a marked effect on system 
cutting performance. When the pick information is 100 percent correct 
an RMS cut error of 0.68 in, is achieved. However, if the pick information 
is 90 percent correct the RMS cut e'"ror is 1.24 in., when it is only 
80 percent correct the RMS cut error grows to 2.33 in. The reason for 
this sensitivity is that when the VCS is under pick control, errors in 
the pick information cause an immediate error in the cutting drum command, 
hence directly affecting system performance. It is, therefore, concluded 
that the sensitized pick should be between 80 to 90 percent accurate 
to enable taking all of the coal, while maintaining the rock taken to 
a minimum. These results are given in Table 3-10. 

2 • 1 • 2 . 4 Effects of Actuator P.ynamics of VCS Performance 

Inserting the nonlinear actuator dynamics into 
the VCS control system had little effect on system performance when a 
given amount of coal is to be left, and hence a CID is employed. However, 
if all of the coal is to be taken, and the VCS is under pick control, 
then the insertion of actuator dynamics has a more marked effect on the 
performance achieved with the RMS cut error going from 0.68 inches without 
actuators to 2.06 inches when the effect of actuator dynamics are con- 
sidered. The reason for this sensitivity is that when the sensitized 
pick is employed it actually measures the presence of coal or shale at 
the cutting point. Therefore, if the actuators respond quickly, a very 
small RMS cut error is realized (i.e., 0.68 inch RMS). Adding the ac- 
tuator dynamics has the effect of introducing lag into the system and 
hence performance degrades. When the VCS is under CID control the errors 
introduced by CID noise and its being mounted behind the cutting drum 
overshadow the effect of the actuator dynamics thus minimizing its effect 
on overall system performance. These results are outlined in Table 3-11. 
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2. 1.2.5 Effects of Varying VCS Traverse Speed Across 
the Face 


The nominal traverse speed at which the baseline 
system was designed and its performance evaluated was 30 ft/min, In- 
creasing this speed to 50 ft/min or decreasing it to 20 ft/min had no 
effect on system performance. The reason for this result is that a 0.25 
second output inter»^al is being used for the CIO which is quite fast 
and within the speed range of 20-50 ft/min the shearer does not travel 
an appreciable distance between CIO information updates. This would 
imply that if the CIO output interval were increased one should see an 
effect of varying shearer traverse speed on overall system performance 
with performance degrading the faster the shearer traverses the face. 
This effect was observed when '^he CIO output interval was increased to 
one second and the results are tabulated in Table 3-12. 


2.2 VCS Performance with Optimum Filtering/Predictive Algorithms 


Since the CIO is mounted behind the cutting drum it is not 
measuring the coal/stale interface at the cutting point. Using the 
measured depth of the coal/shale interface at the CID location and saying 
that this is the depth at the cutting drum location obviously is in error 
since the coal/shale interface changes as one proceeds down the face. 
Therefore, one would suspect that using some form of filtering/predictive 
algorithm syste^'n performance could be improved, since a more accurate 
estimate of the coal/shale interface both at the point of measurement 
and at the cutting drum could be determined. There are many types of 
f ntering/predictive algorithms one can formulate, all of which to 
varying degrees make use of past as well as the present CID measurements 
in the determination of the coal/shale interface. At first glance one 
would be tempted to make use of all or many past CIO measurements in 
order to obtain a “best estimate" of the coal/shale interface at the 
present measurement location and cutting point. Implied in such a scheme 
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would be a method by which all of the CIO measurements could be referenced 
with respect to each other through some common reference frame. However, 
no such frame exists in the longwaVI system and there doesn't appear 
to be any reasonable way by which such a frame could be obtained. There- 
fore, if an algorithm were formulated that used many past CID measurements 
to aid in the determination of the coal/shale interface at the present 
measurement point and the cutting drum, any contribution in these measure- 
ments due to the floating reference frame in which the measurements 
are taken would go undetected. The frame, with re‘'^pect to which the 
height of the coal/shale interface is measured could be considered the 
shearer skid plane. Therefore, measurement changes due to floor undula- 
tions could not be differentiated from changes due to actual variations 
in the coal /shale interface height above the skid plane. Therefore, the 
contribution due to floor undulations in the measurement obtained would 
be attributed to coal/shale interface variations and the filtering and 
predicive alogorithms would yield answers that would be in error by ap- 
preciable amounts. Optimum filtering and predictive algorithms that 
have extremely good potential of determining and predicting the coal/ 
shale interface by using many past CID measurements are generally quite 
sensitive to unmodeled processes that were not considered in their formu- 
lation which in this case would be a floor undulation. In such a case, 
performance degrades rapidly to where it is actually worse than that 
which would be obtained if these algorithms were not employed. Given 
this situation the filtering/predictive algorithms that will be considered 
for the longwall are those which employ CID measuremer^ts in the immediate 
vicinity of the present measurement and cutting points. These types 
of algorithms would require a priori knowledge of the coal/shale inter- 
face statistics in terms of correlation function, mean and variance as 
well as a description of the sensor noise in terms u* its mean and vari- 
ance. The reason for this a priori knowledge requireme. t is that the 
statistics of the coal/shale interface could not be derived from past 
CID measurements due to the lack of a common reference frame. 






?-12 


Two types of optimal filtering/predictive algorithms employing 
CIO measurements In the immediate vicinity of the present measurement 
and cutting points were formulated. One algorithm (Scheme 2) uses the 
present CIO measurements in conjunction with CIO data In the immediate 
vicinity of the present measurement to estimate and predict the coal/shale 
interface at the measurement point and cutting drum respectively. The 
other technique (Scheme 1) uses the present CIO measurement in conjunc- 
tion with estimates of the coal/shale interface in the immediate vicinity 
of the present measurement point for estimation and prediction. These 
two types of algorithms were formulated in a one-dimensional sense using 
present cut CIO data/estimates, and in a two-dimensional sense using 
present and past cut CID data/estimates yielding four different optimal 
filtering algorithms for possible use in the VCS control system. 

Analyses performed on these algorithms indicated that both 
should give appreciable performance improvement if the floor didn't un- 
dulate, with Scheme 1 being better than Scheme 2. In addition, it was 
shown that the algorithm employing CID measurements only would be less 
sensitive to the unmodeled floor undulations than the scheme employing 
CID data in conjunction with coal/shale interface estimates. These re- 
sults were verified by simulation for the two one-dimensional algorithms, 
and the data indicated that cutting performance degraded due to floor 
undulations to the point where it wasn't any better, and in most cases, 
worse than that obtained with the baseline measure and cut algorithm 
(see Table 3-15). The degradation in performance due to floor undula- 
tions would be even worse for the two-dimensional version of these al- 
gorithms and hence were not simulated. Therefore, the conclusion of 
this investigation was that there is no advantage to be gained by using 
optimal filtering and predictive algorithms for VCS control. Thus, the 
simple measure and cut control algorithm defined in the baseline system 
should be used for cutting drum control. 


reproducibility op the 
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2.3 Use of Impact Penetrometer/Reflectometer for VCS Control 

Investigations were conduct ,d into the use of an impact pene*- 
trometer/ref lactometer as the coal/shale presence sensor, i.e., to replace 
the sensitized pick, The prime difference between the two types of coal 
presence sensors from an overall systems viewpoint is the manner in which 
they are integrated into the VCS. The sensitized pick, which is mounted 
on the cutting drum measures whether the drum is presen''iy cutting coal 
or rock, while ti,^ impact penetrometer/ref lectonieter is mounted at the 
approximate CID location, and hence can only measure whether the drum 
has cut coal or rock. This difference is basic in determining the VCS 
performance obtained when using these two coal presence sensors. The 
sensitized pick yielded excellent cutting performance (0.68 inch RMS) 
and YCS seam following capability when used by itself since it does 
measure what the drum is cutting virtually instantaneously. However, 
the impact penetrometer/ref lactometer being mounted some 2.5 feet behind 
the cutting drum cannot measure whether the drum is cutting rock or coal 
until the shearer has traversed 2.5 feet to the present drum location, 
which at 30 ft/min shearer traverse speed takes 5 seconds. During this 
total time, if the impact penetrometer/reflectometer is measuring rock, 
it is continuously depressing the cutting drum or if it is measuring 
coal it is continuously commanding the cutting drum to rise, thus giving 
rise to severe system limit cycles and an RMS cut error of 21.7 inches. 
When the impact penetrometer/reflectometer is augmented by a last cut 
follower, system limit cycles are limited by the deviation allowed between 
the present and previous cuts. If, as in our case, this value is set 
at +2 inches the ensuing RMS cut error is 1,83 inches which is essen- 
tially the RMS of a limit cycle ranging between + 2 inches. These results 
are outlined in Table 3-10 with the actual system cutting performance 
shown in Figures 3-65 and 3-66. Therefore, the conclusion of these 
studies are that the impact penetrometer/reflectometer is not a very 
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satisfactory coal presence sensor for our longwall application, parti- 
cularly when all the coal is to be taken, and hence is not presently 
recomnended, 

2.4 Use of Natural Radiation CIO for VCS Control 


The use of a natural radiation CIO instead of an active 
nucleonic one was investigated. The prime difference between these two 
sensors that affects system performance is the number of counts per unit 
time for a given coal thickness. The active nucleonic CID has many more 
counts per second (approximately a factor of 20) for a given coal thick- 
ness than a reasonably sized (i.e., 3 to 5 inch detection crystal) natural 
radiation ICO. The counts per unit time as a function of coal depth 
for both types of CIDs are shown in Figures 3-11 and 3-12. Therefore, 
the noise and hence the coal depth measurement uncertainty is worse with 
the natural radiation than with the active nucleonic CID and degraded 
system performance when using the natural radiation CID should be ex- 
pected. This indeed has been verified by simulation with RMS cutting 
error increasing from 1.64 inches to 2.25 inches when the active nucleonic 
CID was replaced by the natural radiation sensor. These results are 
documented in Tables 3-5 and 3-11. Therefore, from a performance view- 
point the active nucleonic CID is recommended for system implementation. 
However, there are drawbacks to the active nucleonic sonsor ranging from 
sociological, (i.e., active radiation source) to technical (i.e., sensi- 
tivity to air gaps) to potentially severe integration problems (i.e., 
sophisticated suspension that must keep the sensor against roof with 
a minimum of air gaps, deployment and stowage mechanisms, etc.) that 
the natural radiation CID does not have. During preliminary design (i.e.. 
Phase II) these problems will be addressed in detail. Should these in- 
vestigations indicate that the problems related to the active nucleonic 
CIO precludes its use or detracts significantly from the merits of its 
use, the natural radiation CID would be used instead. Although, as 
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pointed out above, system cutting performance degrades to a degree by 
its use, however, the performance achieved is still quite good when com- 
pared to what is presently being obtained manually. In addition, all 
of the benefits with respect to the health and safety that one gains 
by longwall automation would remain intact when using a natural radiation 
CIO. 


2.5 Slaving of Bottom Drum 

Since it became clear early in the study that it would not 
he feasible to mount a CIO on the floor for controlling the drum cutting 
the bottom due to debris and water that will be there, an alternate tech- 
nique of control for the bottom drum had to be devised for the condition 
where a given amount of coal is to be left on the bottom. (When all 
the coal is to be taken the bottom cutting drum can be controlled by 
a sensitized pick in the same manner that roof cutting drum control is 
achieved.) The technique used was to slave the bottom cutting drum 
to the top one by maintaining the instantaneous height of the cut constant. 
A number of ways of accomplishing drum slaving were investigated. The 
first of these was to maintain the height between the two cutting drums 
constant by using the readouts of the linear actuators controlling both 
drums. This implementation was quite sensitive to floor undulations 
and did not yield satisfactory system performance. Another implementation 
that was investigated was to us«- CID measurements to determine the height 
of the roof and command the bottom drum to maintain the instantaneous 
cut height constant. This implementation, although better than maintain- 
ing the instantaneous height between the cutting drums constant, still 
yielded performance that wasn't any better than locking the drum for 
a flat floor, and degraded when the actual floor undulations were in- 
troduced. Examination of Table 3-2 indicates that four of the six mines 
surveyed exhibited a high degree of correlation between the top and 
bottom coal/shale interface profiles. This phenomenon suggested that 
if the CID measurements were delayed until the bottom cutting drum was 
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directly below the CID measurement point and then used to control the 
bottom drum, better slaving performance should be realized. This was 
indeed the case (for a flat floor) with the improvement being quite ap- 
parent when the top and bottom interfaces were highly correlated such 
as in the Robinson Run Mine. It was also true to a lesser degree in 
a mine that did not have a high degree of correlation betv/een top and 
bottom such as in York Canyon. However, when the floor undulations were 
introduced, significant performance degradation resulted and better system 
performance would be achieved by locking the bottom cutting drum. The 
results discussed above are tabulated in Table 3-16. 

# 

One basic problem (besides performance) that the slaving im- 
plementations discussed above have in common is that they do not tend 
to maintain the difference between two successive bottom cuts within 
a prescribed value. Slaving implementations that would tend to maintain 
two successive bottom cuts within specified limits would be to slave 
the bottom drum to the last cut follower or to a present top cut fol- 
lowing mechanism. Slaving the bottom drum to the last cut follower 
mechanism was tried and it was found that the system performance was 
sensitive to floor undulations. An attempt was made to slave the bottom 
drum to a last cut following mechanism located in the middle of the 
shearer in order to minimize the effect of floor undulations. This 
Implementation did indeed improve slaving performance to acceptable 
levels; however, roof cutting performance degraded appreciably. Slaving 
the bottom drum to a separate present top cut following mechanism was 
tried and the performance obtained was the ''best'' of all of the imple- 
mentations evaluated. Having a separate present top cut following mech- 
anism above the floor cutting drum essentially eliminates the problems 
due to floor undulations and is presently the recommended technique for 
bottom drum slaving. These results are outlined in Table 3-17. 


2.6 Recontnended VCS System 


As a result of the studies performed, the recommended VCS im 
plementatlon is outlined below as a function of the mining capabilities 
desired. 


a. Leaving a prescribed amount of coal on top and bottom. 

1) Top Drum Control 
Active Nucleonic CIO 

Two sensitized picks mounted 180 degrees apart on 
cutting drum 

Last Cut Follower 

2) Bottom Drum Control 
Present Top Cut Follower 

Two Sensitized Picks Mounted in the Cutting Drum 180 . 
Degrees Apart 

b) Leaving a Prescribed Amount of Coal on Top Taking all Coal 
on Bottom. 

1) Top Drum Control 

Active Nucleonic CIO 
Two Sensitized Picks 
Last Cut follower 

2) Bottom Drum Control 
Two Sensitized Picks 








2-18 


c) Taking all of the Coal on Top and Leaving a Prescribed 
Amount of Coal on Bottom 

1) Top Drum Control 

Two Sensitized Picks 
Last Cut Follower 

2) Bottom Drum control 

Two Sensitized Picks 

Present Top Cut Follower Mechanism 

d) Taking all of the Coal on Top and Bottom 

1) Top Drum Control 

Two Sensitized Picks 
Last cut Follower 

2) Bottom Drum Control 
Two Sensitized Picks 

The recomnended systems for VCS control are shown in tabular 
form in Table 3-20. 

2.7 Performance of Recommended System 

The average performance of the recommended VCS configuration 
is shown in Table 2-1. Examination of this table indicates that the 
reconmended VCS configuration will outperform a human operator by at 
least a factor of two at a shearer speed as high as 50 ft/mi n. While 


Table 2-1. VCS Performance Summary 


AVERAGE RMS ROOP AVERAGE R>fS FLOOR 
CUT ERROR (IN) CUT ERROR (IN) 


AUTOMATED VERTICAL CONTROL 

I^CLEOHIC CID 1.36 


2.84 


AUTOMATED VERTICAL CONTROL 

NATURAL RADIATION CID 2.25 4.63 

5 IN. CRYSTiU. 


* MANUALLY OPERATED VERTICAL 

CONTROL SYSTEM N/A 


N/A 


* AVERAGE OP ALL MEASURED PERFORMANCE 
12 PACES IN 10 DIFFERENT MINES 
A.D, LOTLE SURVEY (DEC. 1976-JUNE 1977) 

** 3.1 IN. MEAN ERROR 

NBEARBR VELOCITY "• 30 FT/MIN 

RESULTS ESSENTIALLY UNCHANGED AT 50 FT AlIN 
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AVERAGE RMS 
CUT ERROR (IN) 

2.1 

3.44 

** 5.5 


under manual control the shearer can only move approximately 10 to 12 

ft/mln. These results Indicate the increased productivity and hence 

economic potential of automating the longwall, in addition to the health and safety 

gains. 


2.8 Roll Control System 

The function of the roll control system is to maintain the 
attitude of the longwall s'tearer about its longitudinal (roll) axis at 
a specific value regardless of the twists that may be in the conveyor 
track on which the shearer is riding. This mode of control is necessary 
If the shearer is to be maintained perpendicular to the longwall panel 
being mined and hence within the coal seam as repeated passes are made. 

Two control loop implementations were investigated. One has 
the control loop closed around the inclinometer, i.e., the inclinometer 
is within the closed loop. The advantage of this implementation is that 
there is no need to instrument the hydraulic ram. However, since the 
inclinometer is within the closed loop any filtering that might be added 
to eliminate the effects of cross-axis acceleration will affect overall 
control loop stability and hence is an added constraint on the filter 
design. The other implementation closes the control loop around the 
hydraulic actuator with the inclinometer being external to this loop 
and furnishing what essentially can be looked at as an update signal. 

This Implementation has the advantage of eliminating the dependence of 
control loop stability on the filter characteristics that may be used 
to eliminate the effects of cross-axis accelerations, thus enabling the 
filter design to proceed from these considerations only. However, this 
implementation requires that the linear hydraulic actuator be instrumented. 

[TY of 
1,1 IS 
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A Gomputer simulation for the roll control system was developed 
and consisted of the following: 

a. Full nonlinear model of the linear hydraulic actuation 

system. 

b. Second order Inclinometer dynamics 

c. Cross axis accelerations due to the cutting forces and 

torques as the shearer proceeds along the face. 

This simulation was used to evaluate the performance character- 
istics of the two roll control system Implementations described. 

2.8.1 Roll Control Loop Bandwidth Determination 

In order to choose the parameters for both control 
loop Implementations the desired loop bandwidth had to be determined. 

This was done by estimating the maximum Input frequency that the roll 
control loop would have to respond to. Assuming that two adjacent 5 
foot conveyor sections are twisted In opposite directions, an Input 
frequency of 0.05 Hz results for a shearer traverse speed of 30 ft/mln. 

Hence the desired roll control loop bandwidth was set at 0.05 Hz. 

In order to choose system parameters so that the desired 
0.05 Hz control loop bandwidth would be realized, system frequency responses 
were taken at various amplitude inputs ranging from the maximum commanded 
roll angle of + 5 degree to one quarter this maximum. The reason for 
taking frequency responses at different input amplitudes is due to system 
nonlinearities which make system response characteristics a function 
of Input amplitude. The gains that were finally chosen for both control 
loop Implementations were such that a control loop bandwidth of 0.05 
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Hz would be realized at a 2.5 degree input, or one-half the maximum possi- 
ble command. This was judged to be adequate since it is highly unlikely 
that the roll control loop would be required to command a larger value 
than 2.5 degree. In addition, the loop bandwidth also falls off due 
to flow constraints and attempting to achieve the desired loop bandwidth 
at 5 degree Inputs would probably require a hydraulic modification which 
isn't desirable. 


Frequency responses for the two control loop implemen- 
tations considered are given in Figures 5-4 through 5-6 and 5-8 through 
5-10. 


2.8.2 Roll Control Loop Performance Characteristics 

Both control loop implementations perform well (i.e., 
RMS errors between 0.05 to 0.1 degree) in the absence of cross-axis ac- 
celerations, with the implementation that closes a control loop around 
the actuator being slightly more accurate. However, system performance 
for both implementations degrade significantly even when relatively small 
values of cross-axis accelerations are introduced, as outlined in Table 
5-1. The reason for this sensitivity becomes clear when one examines 
the value of RMS g represented by a 10 percent error in tracking a 1.77 

degree (i.e., + 2.5 degree peak) input signal. This value is approxi- 
-3 “ 

mately 3x10 g RMS which means that if the sensed accelerations are 
in error by approximately 3 milli g's a ten percent tracking error should 
be expected. Or alternately, if the cross-axis acceleration sensed 
by the inclinometer mounted on the shearer has an RMS level ot approxi- 
mately 3 milli g's within the control loop bandwidth at least a ten per- 
cent tracking error will result. In addition, it should be noted that 
if the cross axis RMS g levels within the control loop bandwidth are 
in the order of three milli g's or greater, filtering of the inclinometer 
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output win not yield appreciably better performance if the control loop 
bandwidth Is to be maintained at 0.05 Hz. This result was verified by 
simulation when the Inclinometer bandwidth was reduced from 2.4 to 0.24 
Hz to simulate the effect of filtering. These results are also outlined 
in Table 5-1. 


The conclusion of this portion of the study is that 
since both of the control loop implementations yield almost identical 
performance for all operational conditions, the open-loop actuator Im- 
plementation (i.e., the inclinometer within the feedback loop) is re- 
conmended since it eliminates the requirement for instrumenting the 
hydraulic ram. However, due to the sensitivity of roll control system 
performance to cross-axis accelerations, it is recommended that measure- 
ments be taken of the cross-axis acceleration levels on a longwall shearer 
while mining coal, particularly within and around the expected roll con- 
trol loop bandwidth. It is possible that when these measurements become 
available that it might be more economical to remove the inclinometer 
from the feedback loop and implement a position loop around the actuator 
in order to better filter the inclinometer output. 

2.9 Yaw Alignment Systems 

The function of the yaw alignment system is to automatically 
advance the conveyor and roof supports in such a manner that the face 
"straightness" remains within acceptable limits. The "baseline" system 
defined to accomplish this control function consists of an angle cart, 
diagrammatically shown in Figure 4-5, which is integrated with the shearer 
and measures the angles between respective conveyor sections as the 
shearer proceeds along the face. (The conveyor geometry and the angles 
measured by the angle cart are shown diagrammatically in Figure 4-6.) 

Once a set of measurements are taken the present conveyor profile is 
computed. A set of conveyor advance commands are then computed and sent 


to each roof support, such that if perfectly implemented, would result 
In a totally straight conveyor. This procedure is followed (i.e., angle 
cart measurements taken, conveyor profile computed, and appropriate con- 
veyor advance commands computed) on every advance. 

Various algorithms for computing the conveyor profile from 
the angle cart measurements were investigated. These algorithms were 
a basic angle cart measurement scheme and an optimal weighting scheme. 

In addition the effect of not being able to measure the angle between 
a number of conveyor sections on either end due to the inability of moving 
the shearer far enough into the headgate and tailgate was also investi- 
gated. 


An alternate system implementation to the baseline defined 
was also investigated. This implementation consisted of a directional 
gyro mounted on the shearer which measures the orientation of each con- 
veyor section with respect to a defined reference. These measurements 
are then used to compute the present conveyor profile and a set of con- 
veyor advance commands are subsequently generated that if implemented 
perfectly would straighten the conveyor in much the same manner as in 
the angle cart implementation. 

Two computer simulations were developed in order to evaluate 
the performance of the yaw alignment system. The first was a Monte-Carlo 
simulation which evaluated the accuracy with which the conveyor profile 
could be computed in the presence of measurement error. This was accomp- 
lished by performing 500 measurement runs using particular values of 
measurement errors and computing the conveyor profile for each run. 

The RMS error between the computed and actual conveyor profiles was then 
computed and the average RMS error was determined for the 500 runs made. 
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The average of the RMS error for the BOO runs made, called the "mode of 
RMS", was used as the performance measure of the measurement scheme 
being evaluated In the presence of measurement error. A block diagram 
of this simulation 1s shown In Figure 4-7. 

The other simulation developed was the Yaw Advance Simulation. 
This simulation actually made repeated advances of the conveyor in the 
presence of system measurement, conveyor placement, and roof support 
pull-up errors, and yaw advancement stability evaluated. In order to 
evaluate whether yaw advancement is proceeding stably, twenty advance- 
ments of the conveyor were made in the presence of a particular set of 
system errors. System instability was defined when a particular roof 
support was commanded to move the conveyor backwards, which is impossible, 
within the twenty advances made. If this condition resulted for a par- 
ticular set of system errors, the system error vector was perturbated 
slightly and two more sets of twenty advances made. The average value 
of the error vector for the set of three twenty advances were then used 
to define the stability boundary. The yaw advancement simulation block 
diagram is shown in Figure 4-12. 

2.9.1 Angle Cart Performance Characteristics 

The basic angle cart measurement system uses two sets of measure- 
ments made by the angle cart. One set of measurements are taken when 
the angle cart is totally in a single conveyor section. These measure- 
ments are used to estimate the bias or constant error in the resolvers, 
and are used in correcting the second set of measurements which are the 
angles between contiguous conveyor sections. The angle that the first 
conveyor section makes with the desired reference (i.e., line perpen- 
dicular to headgate and tailgates) is computed via the relative angle 
measurements between the conveyor sections in conjunction with the 
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knowledge where the ends of the conveyor are with respect to each other. 
This is obtained by measuring the distance on both ends of the conveyor 
with respect to surveyed points within the headgate and tailgate. 

The results with the basic measurement scheme indicate that 
the mode of RMS increases linearly with increased measurement error as 
shown in Figure 4-8. From a stability viewpoint the maximum allowable 
one sigma measurement error that can be tolerated in the presence of 
a 0.1 ft (i.e., 1.2 in.) one sigma roof support pullup error and a one 
sigma conveyor placement error of 0.02 ft is 0.03 degree. These results 
remain unchanged even if the angle between one conveyor section on either 
end cannot be measured and hence assumed to be zero. These results are 
outlined in Figures 4-14 through 4-17. It should be noted that the angle 
cart employs eight speed resolvers to make the angular measurements. 

The accuracy of an eight speed resolver is in the order of two arc minutes 
(0.033 degree). Since 200 measurements are taken for each angle between 
the conveyor sections and averaged, and two resolvers are used to deter- 
mine those angles, the basic resolver measurement accuracy needed to 
achieve a 0.03 degree one sigma measurement accuracy on the angle between 
the conveyors is 20 arc minutes (0.33 degree). This is approximately 
an order of magnitude greater than the two arc minute accuracy of the 
eight speed resolver. Hence, a considerable amount of degradation could 
be acconmodated by the angle cart sysem due to vibrational errors as 
the shearer is traversing the face, while still meeting the required 
measurement accuracies for the angles between the conveyor sections. 

Conversations with mining personnel at Old Ben Mining Company 
where they are presently operating longwall faces indicated that the 
roof support pullback error does not exceed one inch. In addition, a 
ram placement error of 0.02 ft (0.24 inch) represents approximately one 
part in 100 for a 30 inch advance which can easily be achieved with state- 
of-the-art instrumentation. Therefore, it is seen that the expected 
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system errors are within those allowable for stable face advancement 
as determined via computer simulation, making the basic angle cart system 
a viable way of implementing yaw face alignment. 

The weighting and bias estimation algorithm for determining 
the conveyor profile using the angle cart measurements requires that 
the angle between the conveyor sections on both ends be known with respect 
to a line perpendicular to the headgate and tailgate. Once these values 
are known the algorithm makes use of the fact that the conveyor profile 
could be computed with the same set of measurements proceeding from either 
end. Therefore, two coordinates could be determined for each joint of 
the conveyor which then could be combined to obtain a minimum variance 
estimate of the conveyor profile. This technique yields better results 
than the basic angle cart measurement scheme by allowing larger measure- 
ment errors, or, for the same measurement errors yielding a more accurate 
determination of the conveyor profile. However, in order to obtain any 
benefit from the bias and estimation algorithm the angle of the end con- 
veyor sections with respect to the desired reference must be known to 
0.2 degrees or better as shown in Figures 4-18 and 4-19. This would 
require that specialized instrumentation to be located in the headgate 
and tailgate to determine these angles to the required accuracy. There- 
fore, the basic angle cart measurement scheme is the technique recommended 
if the angle cart implementation is used for yaw alignment. 

2.9.2 Directional Gyro Performance Characteristics 

In order to evaluate the performance of a directional 
gyro implementation for yaw alignment and to specify the required gyro 
for satisfactory system performance a directional gyro error model had 
to be formulated. This model is shown diagrammatically in Figure 4-23 
where the errors considered are: 
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1) Low frequency random walk error (white noise through 

the integrator) 

2) Constant drift error (ej^) 

3) High frequency error (white noise through first 

order filter) 

4) Gyro reference error (e^) 

Various types of algorithms for processing the direc- 
tional gyro data in order to determine the conveyor profile were investi- 
gated and are discussed in section 4.5. The results of these investigations 
indicated that the best algorithm for processing the directional gyro 
data is one that uses a single measurement as soon as the gyro is aligned 
in a particular conveyor section. The measurements thus obtained are 
then put through an optimum filtering scheme such that a minimum variance 
estimate of the conveyor profile is obtained. 

Investigations into the allowable system errors using 
the above algorithm with a shearer traverse speed of 30 ft/min, indicated 

that a spectral noise intensity of approximately IxlO"^ -could 

be tolerated in the presence of a gyro reference error of 0.2 degrees 
and a roof support pull up error of 0.1 ft (1.2 inch), for stable yaw 
advancement (see Figure 4-33). This value of spectral noise intensity 
corresponds to a random drift error of approximately 0.6 deg/hr which 
is considered to be a good and relatively expensive gyro. In order to 
increase the allowable random drift error that can be tolerated the gyro 
would have to be towed along the conveyor at a greater speed than 30 
ft/min. This implies that the gyro would have to be mounted on a separate 
cart which would either have its own motive power or possibly be towed 
by the conveyor. In either case, the money saved by allowing the use 

repeoducieility op the 
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of a cheaper gyro would have to go into the cost of manufacturing a separate 
gyro cart. In addition, the use of a separate gyro cart would introduce 
operational problems with respect to integrating the measurement with 
standard longwall operations requiring additional time to make the re- 
quired measurements, and would almost certainly preclude the measurement 
of the angles of the end conveyor sections due to interference with the 
shearer. Therefore, the recommended system implementation for yaw face 
alignment is the angle cart system employing the basic measurement scheme 
described above. 
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3. VERTICAL CONTROL SYSTEM 


3.1 FUNCTIONAL DESCRIPTION 

The function of the vertical control system (VCS) is to control 
the position of the shearing drums so that their cut follows the coal/ 
shale interface (CSI) on the top and the bottom of the mine and leaves a 
prescribed amount of top and botto(n coal. The VCS must also maintain 
the present cut within a prescribed amount of the last cut made. The 
shearing drums are positioned by raising and lowering the ranging arms 
on which they are mounted, using hydraulic actuators. The drums are 
controlled independently and are coupled only through the hydraulic 
system, where the pump flow is shared between the two ranging arm 
actuators . 


Commands to the hydraulic actuators are generated from a canple- 
ment of coal interfaces and drum position sensors. These sensors 
include proportional coal interface detectors (CIDs) which measure the 
coal depth remaining, and discrete CIDs which indicate the presence of 
coal or rock. The position sensors include a last cut follower (LCF) 
which indicates the shearing drum's position relative to the last cut 
made, and an actuator position sensor v\^iich allows the position of the 
ranging arms relative to the shearer body to be determined. The propor- 
tional CIDs are radiation detectors, either sensing natural background 
radiation or backscatter from a radiation source. The discrete or coal 
presence sensors use the optical and hardness properties of coal and 
rock to detect the difference between them, hence the presence of 
either. 
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The sensors are combined in a hierarchical manner to develop actu- 
ator commands. The LCF has the highest priority and will never allow 
the present cut to differ by more than a prescribed amount from the last 
cut. The coal presence CIO has the next priority and will move the 
shearing drum away from the interface if rock is detected. The propor- 
tional CID is the lowest priority and will provide the control for the 
cutting drum when it is within the coal beam and the commands are within 
the prescribed limits relative to the last cut. 

Since the shearer operates in either direction, both drums must be 
capable of cutting either the top or the bottom. The lead drum cuts the 
top while the trailing drum cuts the bottom, and the individual drum 
control laws are therefore changed appropriately, depending on whether 
the drum is leading or trailing. The control laws also change, de- 
pending on the desired coal depth to be maintained. If it is desired to 
take all the coal, only the coal presence sensors are used for control. 

It appears impractical to use proportional CID sensors on the 
bottom, due to loose coal, rock, and water which collect there. The 
discrete coal sensors provide the capability of taking all the bottom 
coal, but when it is desired to leave the bottom coal, an alternative 
control scheme is used. The rear drum is slaved to the top drum, thus 
providing the capability of maintaining a constant cut height. This 
leaves a prescribed amount of bottom coal. In addition, slaving of the 
bottom drum will render the system relatively insensitive to undulations 
in the conveyer, thus maintaining the shearer in the coal seam. 
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3.2 DETAILED MATHEMATICAL MODEL OF THf: VCS 


3.2.1 Overall System Description 

The VCS simulation model, shown in Figure 3-1, includes 
the control system, sensors, actuattrs, and shearer dynamics as well as 
the mine geometry. Asterisks indica;e those portions of the figure 
which are part of the simulation, but are not part of the actual shearer 
or control system. Each portion of the diagram will be discussed in 
detail in Sections 3.2.2 through 3.2.7. 

The chassis location module represents those calculations 
necessary to detenaine the position and orientation of the shearer chas- 
sis relative to the mine. This includes keeping track of the chassis 
position along the face and the chassis pitch angle a. The mine model 
is that portion of this simulation which represents the mine geometry. 
Actual coal seam profiles were digitized and used to realistically 
simulate the conditions in which the shearer operates. 

The ranging arm geometry calculations use the pitch angle 
o and the arm angular position relative to the chassis e to calculate 
the shearing drum position in the mine. The quantities computed are the 
coal depth presently being cut, the dr im height relative to the CID, P, 
and the drum height relative to this 1 ist cut, Q. The CID module pro- 
cesses the coal depth measurements, introducing the appropriate errors 
to simulate the various CID sensors. The coal depth output measurement 
is designated c. 

The baseline control compares the measured coal depth with 
the desired coal depth bias Cg. This difference is summed with the drum 
position relative to the CID ser sor, resulting in a command az, the 
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Figure 3-1. VCS Simulation Model Block Diagram 














desired change in height of the cutting drum. This information, along 
with the last cut data and present actuator ram displacement is used 
in the control law computations. The result is X^, the desired dis- 
placement of the hydraulic ram. 

Actuator dynamics and the ranging arm dynamics are modeled 
in some detail. The actuator model includes the control valve with in- 
put thresholds, check valves, and pump flow rates. Ranging arm dynamics 
are coupled with the actuator and develop the ranging arm angle p. The 
actuator is instrumented with a position sensor and used in a closed 
loop so it can be position commanded. 

3.2.2 Chassis Location Model 


Figure 3-2 shows the mine coordinate system and defines 
major components of the chassis and mine simulation. The shearer posi- 
tion along the face is taken to be the +X direction, positive to the 
right as viewed from the gob side. The seam height is the Z direction, 
positive upward, and the cut is in the y direction positive into the 
face (into the plane of the y figure). 

The six profiles shown in Figure 3-2 will be discussed in 
Section 3.2.7, which deals with the mine model. These are the interface 
or CSI, the present cut, and the last cut for both the top and bottom. 
The right drum is designated as drum 1 and the left drum as drum 2. 

Drum hub coordinates are designated as and ^Di for i = 1, 2. The 

CID sensors are located a distance Kj^ behind the drum hub. The sensor 
coordinates are X^^ and . Present cut locations are given by (Xq^, 

Zee) ^^FC' ^FC^’ 



Figure 3-2. Mine and Shearer Coordinate System and Parameter Definition 






The shearer chassis contacts the conveyer track at points 
(Xj» Zj) and (Xg* Z 2 ). The motion of the shearer is computed assuming a 
constant velocity in the X direction. That is 

Xj = VT + Xjj 


where T is the integration interval of 0.05 s, and Xjj is the initial 
value of Xj, and V is the chassis velocity giv6n by the profile in 
Figure 3-3. In the figure, 


Tn = 


Tk ■ Tn * T„, 


Lp = length of face to be cut (ft) 

. maximum velocity PAGE IS POOR 


The quantity T^^ is the time to accelerate to the maximum velocity 
which in most cases was 0.5 ft/s (30 ft/min). The value of T^^^ used was 
6 s, and the deceleration time Tp - Tj^ is equal to The time Tp is 
the total cutting time of a given simulation run. 

The chassis pitch angle a is computed from the difference 
in 2 coordinates at the chassis contact points. The shearer chassis is 
illustrated in more detail in Figure 3-4. Alpha is computed as 

Xg = Xi + C cos a 
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Zj, I 2 • 1 coordinate at Xj and 


•i 

a = sin 



Finally, Z coordinates of the shearer contact points are 
computed from the track. The track simulates the effect of placing the 
5 ft conveyer sections on the bottom just cut. This is necessary 
because the 5 ft conveyer section effectively filters the bottom cut and 
provides a smoother surface than the cut itself. The procedure for 
determining the track is recursive, placing the sections one at a time 
on the bottom cut. The procedure is as follows: 


Let 

Xs> Zg = coordinates of start of a conveyer section 
Xp, Zp = coordinates of the end of the section 


The above points may or may not be on the bottom cut. An initial value 
for Xp and Zp is found by searching the bottom cut profile until 

(Xg - Xp)^ + (Z 5 - Zp)^ > 25 

Once the probable conveyer endpoint is found, it is necessary to deter- 
mine if any point on the cut is higher than the endpoint (Xp, Zp). The 
highest point on the cut between (Xg, Zg) and (Xp, Zp) is found and 
designated (Xj^, Zj^). The following logic is then used: 
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(1) If ^ Zp, then (Xp, Zp) is the actual conveyer endpoint 

(2) If Z^ > Zp 

Y = tan“^ [(Zj^ - Zg)/(X^ - Xg)] 

and Xp = Xg + 5 cos y 
Zp = Z^ + 5 sin Y 


The above procedure computes the conveyer coordinates. If the point 
(X|^, Zj^) is higher than (Xp, Zp) , then it must detennine the slope y and 
new endpoint of the conveyer. The procedure repeats with (Xg “ Xp) and 
(Zj = Zp) until all sections have been placed. 

3.2.3 Ranging Arm Geometry 

Orientation of the ranging arms relative to the shearer 
chassis is shown in Figure 3-4. Pivot points of the anns are located a 
distance 2a (163 in.) apart and a distance b (35.25 in.) above the con- 
veyer pan. The ranging arm angles relative to the shearer are and 
02. Note that 0 ^ is positive in a counterclockwise sense, while 02 
positive clockwise. Chassis and ranging arm parameter values are given 
in Table 3-1. 


Figure 3-5 illustrates details of the right ranging arm. 
The shearing drum hub is located t (74.5 in.) from the arm pivot point, 
at an angle 0 ^ (4.04 deg) from the horizontal. The nominal actuator 
displacement X^ is 23.25 in. The actuator requires an additional 
displacement of 3.255 in. to raise the arm 4.04 deg and set 0 to zero. 
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Table 3-1. VCS Parameter Values 



W *= 28.324 in. 

0 “ 81.0 in. 
c = 35.25 in. 
a " 81.5 in. 
t = 74.5 In. 

Xq » 23.25 in. 
r = 27.0 in. 

^ = 69.3256 deg 

6|^ = 4.041 deg 







Figure 3-5. Ranging Arm Diagram 



Figure 3>6 shows the geometry of arm elevation. The triangle of 
Interest is formed by W, D and Xq From the law of cosines 

(Xo + - 2WD cos (♦ + 0^+0) 

This relationship can be solved for the ranging arm rigid position. 
Thus , 


1 1 - (X^ - X.)‘ 

0 * cos I j - (<ti + 0 ) 

2WD / H 


and 


X^ = y + o'" - 2WD cos ( $ + 0^ + 0) - Xg 


.2 . r .2 


when the arm is horizontal 


9 = -4.041 deg 
X^ = 3.2547 in. 


The range of X^ is between 1.32 and 9.5 in. for 
X^ = 1.32 in. 0 = -14.6442 deg 

X^ = 9.5 in. G = 32.2444 deg 
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Once the angles Pj and 62 known, it is possible to 
compute the coordinates of the hubs of the shearing drums (Xpj. and 
(Xo2> Zq 2)» When the shearer chassis rotates an angle a in pitch, it 
does so by rotating about the track contact which is lowest if about the 
point with the smallest Z value. The coordinate equations are as 
follows. 

Xpi ~ a-bsina+t cos (Bj + Ot) 

Zpj = Zj + Aj sin a + b cos a + *. sin (Bj + o) 

Xp2 ~ “ ^2 0 - b sin tt - 4 cos (o - B2) 

Z02 * Z-j- - Ag sin a + b cos o - 4 sin (o - B2) 

The parameters Xj, Zj, Aj, and vary, depending on which contact point 

is acting as the pivot. This is determined by the algebraic sign of a. 
If 


a > 0 


Aj - a + C /2 
A2 - a - C /2 
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= a - c/2 


and for any a < o, 


Ap = a + c /2 


Xf = X2 


Zt = Z2 


There are several sensors associated with the 
shearing drum and the CID. These sensors indicate the elevation of the 
shearing drum relative to the CID and to the last cut made. The 
geometry of the CID sensor mounted on drum 1 is shown in Figure 3-7. 

The CID IS mounted on an arm extending from (X^j, Z^^) to (Xp., Zp,), a 
distance Kp from the drum hub. This arm is maintained parallel to the 
chassis, so that any pitching is reflected in the orientation of the arm 
relative to the seam. The distance Pj is measured by a displacement 
sensor* Ideally, it is desirable to measure the drum height relative to 
the CID in the vertical direction indicated by (Xp^, Zpj) and (Xgj, 

Zgj). However, when a is nonzero, the distance is measured in the 
direction (Xpj. Zp^) to (X^j, Z^j). p^ is therefore 

I 


The coal depth at any instant is ej which is measured vertically and i- 


=1 = Vo 


^S1 “ ^Cl^ * - Zpi) 


SI "cr 
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j. where * Xgj, Figure 3-8 shows the geometry of the last cut follower 

^ for drum 1. The distance Qj is the height of drum 1 relative to the 

last cut. This measurement is also effected by pitching and is given by 



Figures 3-9 and 3-10 show the sensor geometry for drum 2. 
The equations for the measurements are the same. 

3,2.4 CID Sensors 

The coal interface detectors or CIDs are of two basic 
types: proportional sensors which measure coal depth, and discrete 

sensors which indicate the presence of coal or rock. The proportional 
sensors are the nucleonic sensor which measures depth by counting 
backscattered y radiation from a cesium source, and a natural radiation 
sensor which measures depth by counting natural background radiation 
transmitted through the coal. Both of these sensors are mounted behind 
the shearing drum. The descrete sensors are a sensitized pick and a 
combination impact penetrometer- ref lectometer. The sensitized pick 
measures forces exerted on instrumented cutting picks to determine rock 
presence, while the penetrometer- ref lectometer measures surface 
hardness and reflectance to sense rock. The pick is of course mounted 
on the shearing drum while the penetrometer- ref lectometer is mounted 
behind the drum. 

3. 2. 4.1 Nucleonic CID Model 

The calibration curve for a 20 in. CID sensor is 
shown in Figure 3-11. There is a cesium source emitting ganma radia- 












NUCLEONIC SENSOR; 20 IN. SEPARATION 
DETECTOR: CSICRYSTAL 
SOURCE; 30mCiCS-137 
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Figure 3-11. Nucleonic CID Calibration Curve 



tion, which is backscattered by the coal. Therefore, the thicker the 
coal, the more backscatter. Thus, the number of counts per second 
increase with the coal depth. Two basic sources of error are modeled 
for this sensor: a variation in the number of counts received in a 

given inteval, and the effect of air gap. 

The curve in Figure 3-11 illustrates the rela- 
tion between coal depth and the number of counts acquired in 8 s. The 
counts represent the average number of counts received in this interval. 
The number of events occurring in a given interval is Poisson distri- 
buted, with density function 

f(x) = i = 0, 1. 2, ... 

where x is the parameter of the distribution - in this case, the average 
number of counts per unit time. For a given time interval t, the mean 
and variance of the Poisson distribution are 

ji = Xt 


thus, as the interval t is reduced, ihe average number of counts 
decreases as well as the variance. The measurements become noisier, 
however, because the ratio of standard deviation to mean becomes larger 
as T decreases. 


The air gap error occurs whenever the source or 
detector is not in complete contact with the coal surface. In this 
case, radiation short-circuits the coal and is picked up by the sensor 
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directly, thereby resulting in a larger number than would normally be 
expected for a given depth of coal. When the calibration curve shown in 
Figure 3-11 is used to Interpret this data, the coal depth indicated by 
the sensor is larger than that actually present. 

The simulation of this sensor was as follows: 

(1) Coal depth is computed every 0.05 s and averaged for t 
seconds . 

n = t/0.05 

‘A 'Ti Y. 

i=l 

(2/ Calibration curve (X) is used to determine the 
number of counts/second for inches of coal 


^ " ^CAL 


^AGT * 


(3) Poisson distributed noise is added to C^^j with mean 
and variance of Ct 


= ^ACT 
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( 4 ) Uniformally distributed noise on the interval (0 - 250) 
counts/second is added to simulate an air cjap of 

0 - 0.5 in. 

C,^ = + u(0 - 250) 

(5) The calibration curve is again used to interpret Cj^ as 
inches of coal 

The resultant value of is treated as the CID output, 
available every t second. 

3. 2. 4. 2 Natural Radiation CID Model 

The natural radiation sensor measures background 
radiation transmitted through the coal, from the naturally radioactive 
overburden. Figure 3-12 shows the calibration curve for a 3 in. crystal 
detector without shielding. Since the radiation is transmitted rather 
than backscattered, the number of counts decreases as coal depth in- 
creases. The number of counts is considerably lower than that with the 
nucleonic sensor and, as a result, the natural radiation sensor is in- 
herently noisy. Since there is no radiation source, air gaps have 
little effect on the sensor performance. The simulation is the same as 
with the nucleonic sensor, except that the natural radiation calibration 
curve is used, and air gap noise is omitted. 
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3. 2. 4. 3 Sensitized Pick Sensor 


The sensitized pick sensor is implemented by 
instrumenting picks on the cutting drum. The picks are instrumented 
with force transducers whose outputs are analyzed to detect the signa- 
ture of coal, rock, or air. A minimum of two picks is used, spaced 180 
deg apart. Thus, there is a pick latency of 0.66 s based on a nominal 
drum opening of 45 rpm. That is, new data is available only every 0.66 s. 
Since it is desirable to measure rock presence, it is necessary to 
interrogate the picks only when they are in a region near the interface. 

The simulation model determines the presence of 
rock, coal, or nothing, under the pick in an approximately 15 deg region 
near the top of its rotation. Figure 3-13 shows the point examined on 
the top drum. These points are vertical from the drum hub, 3 and 6 in. 
in front of the hub. That is 


Xpj + 3 in., Zpj + 26.833 in. 

Xpj + 6 in., Zpj + 26.325 in. 

Each of these points is examined to determine if they are in rock or 
coal. It is also possible that this area may have been previously cut 
and the picks are therefore in air. 

The sensitized pick sensor model assumes that a 
fixed force threshold would be established and any output above that 
level would be rock. Also, the lower forces resulting from cutting coal 
or air yield the same output. The simulation model outputs a 0 or 1 
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every 0,66 s, depending on whether rock or coal/air is sensed, respec- 
tively. A 0 is outputted if any of the 3 test points indicate rock. 

The output of the sensor is assumed to be a 2 when it is interrogated 
and the data is not fresh. To simulate errors in detection of coal or 
rock, the 0 or 1 output which is modified to be correct a percent of the 
time is, 90 and 80 percent. 


3. 2.4. 4 Impact Penetrometer-Reflectometer Sensor 

The impact penetrometer- ref lectometer detects 
the presence of coal or rock by sensing the optical and hardness proper- 
ties of the present cut surface. The penetrometer is really not a 
penetration device, but rather taps the cut surface. The deceleration 
of the impact device can be interpreted to indicate the hardness of the 
surface. Since rock is harder than coal , it can be detected by this 
device. The ref 1 ectometer measures light return from the cut surface, 
and since there is a reasonable difference in color and reflectance 
between coal and rock, the two materials can be detected optically. The 
sensor, illustrated in Figure 3-14, contains a penetrometer and two 
ref 1 ectometer which are voted to determine the sensor output. 

This sensor was simulated by determining the 
presence of rock or coal at the CID location, i.e., behind the shearing 
drum. The simulation model would output a 0 for rock and a 1 for coal. 
It was assumed that the sensor could provide data as often as required. 
The percent correctness could also be varied to simulate sensor errors. 



Figure 3-14. Penetrometer - Ref 1 ectometer Sensor 







3.2.5 Actuator System 


A diagram of the VCS hydraulic system is shown in Fig- 
ure 3-15. A single hydraulic pump, regulated to 2200 psi , supplies all 
the hydraulics on the shearer. Each ranging arm is actuated by a 
hydraulic cylinder, which is controlled by a control valve through a 
dual pilot check valve and a needle and check valve. The basic system 
is a rate conmand system where a control valve displacement causes the 
cylinder to move at a constant rate, determined by the amount of valve 
displacement. In manual operation, the ranging arm height is adjusted 
via the control valve, and held in place by the dual check valves, when 
control is stopped. The needle and check valve are used to compensate 
the system when lowering the arm. The cylinder piston has different 
areas on each side because the actuator ram is attached to one side. 

This causes the rate at which the cylinder retracts to be greater than 
the rate at which it extends. The needle valve is adjusted to restrict 
the flow on retraction and equalize the two rates. The needle valve 
also prevents chattering when the arm is lowered. 

Figure 3-16 illustrates the simulation model of the 
hydraulic actuator and ranging arm dynamics. The input deadband repre- 
sents an electrical threshold of the control valve, whose input can 
range from ±8 V, with 2.4 V required to produce any valve motion. The 
control valve is modeled by a 6 Hz second order loop with a 0.9 damping 
ratio. This second order loop is also rate and position limited. The 
rate limit is adjusted so that the valve reaches its maximum travel in 
0.25 s. The maximum valve travel is limited to ±1, with the valve being 
fully open at ±1 and closed at zero. The maximum flow rate available is 
8 gallons per minute (GPM), while the control valve capacity is 10 GPM. 
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Figure 3-15. VCS Hydraulic System 


















Thus, the 1.25 gain, limit, and 8 gain simulate this effect. When the 
valve displacement is 0.8, the flow is 8 GPM and limits at that 
valve. 


The gain consists of three parts, which simulate the 
dual pilot check valve, a flow split between ranging arms, and a supply 
regulator delay. The dual pilot check valve is simulated by a deadband, 
where no flow is experienced until the flow rate reaches 5.6 GPM. At 
this time, the valve opens and flow jumps to 5.6 GPM and is proportional 
from there to 8 GPM. This check valve represents the principal non- 
linearity in the system, and in conjunction with the control valve rate 
limit, creates an effective delay. It takes 0.8 s for the control valve 
to open in response to an 8 V command, such that 5.6 GPM flow is achieved, 
resulting in a delayed response to an input. There is an additional 
delay due to the supply regulator. If no flow is required, the pump out- 
put is diverted by the regulator, and when a flow demand occurs, the 
regulator must respond, and restore flow to the system. This procedure 
takes 0.33 s. If one cylinder is in use or has been in use in the last 
0.5 s, this delay does not occur because the pump is already on line. 

The flow split logic simulates the effect caused by simultaneous opera- 
tion of both ranging arms, when the pump flow must be shared between 
the two cylinders. This sharing of flow reduces the ram extension rate 
by a factor of two. However, the retraction rate remains constant since 
the arm weight aids in this direction. Figure 3-17 shows a diagram of 
the gain Ap which varies as a function of load torque to simulate flow 
split. 


The remaining portion of Figure 3-16 comprises the actuator 
cylinder and ranging arm dynamics. The actuator model is a simplified 
linearized hydraulic model assuming constant supply pressure, and a rigid 
coupling between the actuator and ranging arm. Defining the following: 




Figure 3-17. Flow Split Gains 
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C = |a 

X 3x 


C = ^ 
\ 3Pi 


Kq = fluid bulk modulus 


X = valve position 


X^ = ram position 


q = flow rate 


P = load pressure 


A = piston area 


V = volume of fluid 


Dp = ranging arm lever arm 
= leakage coefficient 
T|^ = load torque 


The actuator equation of motion is given as follows 


'' 'a = - ifc ^ - 'S * *1> ^ 


'B P 

The load equation (ranging arm) is: 


T = ADpP = J0 + BO + Tl 


reproducibility of Tm? 
original page is POOR^ 
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Rearranging: 


^ P = CxX - (Cp H. P - A 


J0 = -B6 + ADpP - Tl 


To find the relation between and 0, the expression from Section 3.2.3 
is used here: 



+ 


- ?WDcos( + 0(^ 



X 


0 


or 


(Xa + X^)^ = - 2WDcos(<t) i- 6,^ + 6) 

Differentiating with respect to time yields 

WDsin{(|) + 0j^ + 6) . 

’'a” — 

If the leakage terms are neglected, the resulting equa- 
tions can be diagrammed as shown in Figure 3-16. There is an additional 
second order loop shown in this figure, whose output is subtracted from 
0 to yield the ranging arm elevation angle 6. This loop is a 4.7 Hz 
0.05 damped second order loop used to simulate arm flexibility. The 
amplitude of this added term was chosen so that the simulation output 
agreed with actual test data. 

Figures 3-18 through 3-20 show the open loop response of 
the actuator to step inputs. The input was an 8 V positive step for 
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3 s, zero V for 3 s, a -8 V step for 3 s, then zero again. This is the 
response of a single actuator so the supply regular delay is very evident. 
The total of this and the effective check valve delay is approximately 
0.5 s as shown in the figures. The rate oscillations shown in Figure 
3>19 are due to the flexibility model described above. 

The actuator is used in a closed loop configuration in 
the VCS, such that its displacement can be commanded. The ram is instru- 
mented with a position sensor and a loop closed around it with a gain 
of 170 as shown in Figure 3-21. This gain is such that a 0.15 in ram 
error will generate 2.4 V and thus exceed the control valve threshold. 
Figures 3-22 through 3-24 show the actuator response to a step command 
to extend to 7.5 in. It is important to note the pronounced limit cycle 
behavior present in this response, which is clearly seen in the rate 
time history. Lowering the loop gain will eliminate this behavior at 
the expense of positioning accuracy. Despite this limit cycle behavior, 
the loop was operated as shown in the figures because studies showed 
better mining performance than with Tower gains. 

3.2.6 VCS Control Laws 

The basic VCS control laws, shown in Figures 3-25 through 
3-29, include various combinations of CIO and position sensors. Through- 
out these control laws, a sensor hierarchy was employed. The discrete 
or rock detecting CIDs were given the highest priority of the CID sensors, 
and the shearing drum would be moved away from the CSI whenever rock 
was detected. The proportional CIDs had the lowest priority and provided 
the control for the drum when it is in the seam and within prescribed 
limits of the last cut. The highest priority sensor of all was the last 
cut follower, which maintained the present cut to within +2 in. of the 
last cut. If either the discrete or proportional CID data generated 
drum commands which would cause the cut to deviate more than +^2 in., 
the drum displacement was limited to +2 in. until CID data again com- 
manded the cut within these limits. 
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The control law for the proportional CID alone, shown in 
Figure 3-25, is used for both the nucleonic and natural radiation CIDs. 

The desired drum height change command is simply the sum of the 
measured coal depth e and the drum position relative to the CID P minus 
the desired bias Cg. The P sensor maintains the drum at the same dis- 
tance from the skid plane as the CID, and the coal depth minus bias 
provides the cut error, A^ is computed only when there is current CID 
data available. 

Figure 3-26 is a flow chart of the proportion CID and LCF 
control law. The LCF measurement Q is assumed to be available as often 
as necessary; in this case, every 0.05 s, the iteration rate of the simu- 
lation. If CID data were not current, the LCF was checked for +2 in. 
and the drum commanded, if necessary, to maintain this limit. If CID 
data were available, a A^ was computed and used if lQ| was not greater 
than 2 in., and would not be greater than 2 in. it A^ was used. If Q 
exceeded limits, the A^, was modified to maintain the 2 in. limit in the 
direction determined by CID. 

Figure 3-27 shows the control law for the pick sensor only. 

To use this sensor, the drum is raised until rock is detected, withdrawn, 
and alternately raised until rock is encountered again and the cycle 
is repeated. Controlling in this limit cycle manner allows the shearing 
drum to follow the CSI. The rate of drum rise and withdrawal were adjusted 
to yield a good cutting performance. The control scheme used a nominal 
depression of 0.25 in. and a nominal withdrawal of 0.5 in. If current 
pick data is available, a A^ is computed, depending on whether coal or 
rock is sensed. If the same sensor indication occurs several times in 
a row, the depression command is integrated until the sensor indication 
changes. That is, if f'e pick indicates rock the drum is withdrawn 0.5 
in. If the next reac'^-g is a rock indication, A^ is then set to 1 in., 
then 1.5 in., etc,, until the sensor indication changes state. 
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The control laws 'or proportional CID and pick, and CID, 
pick, and LCF are shown in Figures 3-28 and 3"-29 respectively. These 
are essentially combinations of the previous control laws. When the 
pick is used with other sensors, it determines only the withdrawal of 
the drum; the depression is controlled oy the proportional CID. The 
last cut follower, when used with other sensors, has the highest priority 
and overrides the other commands. 

These basic control laws were used with other sensors as 
well. The same CID schemes were used with the natural radiation CID 
providing proportional information. The pick control law was also used 
with a penetrometer, and the penetrometer with an LCF. Drum slaving 
control, which will be discussed later, uses a control law essentially 
like the CID and pick law with the slav'ng sensor providing the proportional 
data. 


3.2.7 Mine Model 


For the purpose ol Simula ;ion, it was necessary to develop 
a mine model in which to evaluale the V^S control schemes. To achieve 
this, several 1 pass mine surveys were digitized and analyzed to provide 
an information base with which to generate subsequent simulated mine 
data. The statistical data obtained from these mines, including seam 
slopes, C$I variances, and seam heights, are presented in Table 3-2. 

An auto-correlation analysis of the top and bottom of each mine was per- 
formed, along with the corresponding power spectral density. The auto- 
correlation and power spectral densities of the York Canyon and Robinson 
Run mines are illustrated in Figures 3-30 through 3-33. A cross-correla- 
tion analysis was also performed between the top and bottom of the mines. 
The maximum cross-correlations (f the six mines are also shown in Table 
3-2. Figures 3-34 and 3-35 show the cross-correlations of the York Canyon 
and Robinson Run mines, respectively. It can be seen from these data 
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Figure 3-30. Autocorrelation of York Canyon Mine 
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Figure 3-31. Autocorrelation of Robinson Run Mine 
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Figure 3-33. Power Spectral Density of Robinson Run Mine 
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Table 3-2. Mine Statistics 


MINE 

SLOPE 
(IN/5 FT) 

BIAS 

(IN) 

STANDARD 

DEVIATION 

(IN) 

CROSS 

CORRELATION 
BETWEEN TOP 
AND BOTTOM 

SHOEMAKER BOTTOM 

-.0204 

50.01 

4.63 

.91 

SHOEMAKER TOP 

-.0214 

117.7 

4.88 


FEDERAL #1 BOTTOM 

-.0314 

17.24 

2.69 

.75 

FEDERAL #1 TOP 

-.0349 

116.0 

3.81 


ROBINSON RUN BOTTOM 

0.243 

-1.5 

5.18 

.85 

ROBINSON RUN TOP 

0.236 

96.06 

4.57 


MARTINKA BOTTOM 

0.0527 

30.8 

10.84 

.98 

MARTIN KA TOP 

0.057 

86.01 

10.49 


FEDERAL #2 BOTTOM 

-.0397 

51.74 

4.86 

.27 

FEDERAL #2 TOP 

-.0218 

136.41 

2.27 


YORK CANYON BOTTOM 

0.513 

38.6 

3.07 

.24 

YORK CANYON TOP 

.486 

113.7 

2.283 
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Figure 3-34. Cross Correlation Between Top and Bottom of 
York Canyon Mine 





Figure 3-35. Cross Correlation Between Top and Bottom of Robinson Run Mine 







that, statistically, there is a great deal of similarity between the 
mines considered. The CSi variances and correlation data are quite 
similar. In addition, except for the York Canyon and Federal #2 mines 
there was a significant degree of correlation between the top and bottom 
which will be of some importance when drum slaving is considered. 

For the simulation, it was necessary to have more than 
the single CSI profile for a given mine. To do this, a procedure was 
developed to generate random data with the same statistical and spatial 
correlations as the original surveyed data. The profile generation 
scheme is based on the following theorem of conditional distributions 
with multiple correlations: 

^ T 

Given: a. X = (Xp Xg, ..., Xj^) a vector formed from realiza- 
tions of a set of N jointly normally distributed random 
variables (Xp Xg, ...» Xj^) which are all N(0,1), 

b. A = (p^j) an NxN correlation matrix of the N sequence 
elements, 

*r 

c. B - (P|^+]^‘l» ^N+1 2’ ***’ bfi+1 ^ vector (Nxl) of 
desired correlations of the (N+l)‘th point and the pre- 
vious N. 

The theorem states that the density of the (N+l)'th point 
given "X, is N(Xj^^lj ¥^a”^X, 1-B^A“%) or normal with mean ¥^A"^X and 

variance 1 - ¥^A'%. 

The following procedure was used to generate multiple profiles 
from the original data: 
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1, The original data were detrended by fitting a regression line 
and subtracting it off. This results in essentially a zero mean 
process. The standard deviation was then estimated and divided 

out. 

2. Correlation analysis showed that the mine profiles are ap- 
proximately exponentially correlated. 

d 


where d is the distance between points 1 and j and 0 is the 
correlation distance. D was found to be between 25 and 50 ft for 
the mines analyzed. 

3. To generate a new point on a mine profile, all points original 
or generated within D ft of the new point are found, and become 
the elements of "X. 

4. The matrix A and vector F are formed using the spatial relation- 
ships of the selected points and the above correlation function. 

5. A is inverted, the mean and variance found, and a new point 
generated using a random number generator. 

Figures 3-36 and 3-37 show the top and bottom profiles 
for the Robinson Run mine. Five profiles were generated along with the 
original data shown. These are not 3 dimensional plots but rather the 
6 profiles biased apart for plotting purposes. It can be seen that the 
generated profiles show similar characteristics to the original data. 
Figure 3-38 presents a 3 dimensional view of a 12.5 by 50 ft section 
of this simulated Robinson Run bottom. It can be seen that the randomly 
generated profiles form a realistic surface. 
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Figure 3-36. Generated Profiles for the Robinson Run Top 
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The mine data were used in the simulation to provide current 
CID measurements and initial conditions. For a simulation run, the first 
or the original mine profile was used for initial conditions, providing 
the last cut data; and to form the conveyor track. The last cut, taken 
as the CSI on the top minus the desired bias, and the CSI on the bottom 
plus the desired bias, therefore represented a perfect cut. The track 
on which the shearer moved was formed from this perfect last cut. On 
subsequent mining passes within the same simulation run, the last cut 
was that actually cut, and the track formed from the actual cut surface. 

3.3 SIMULATION DETAILS 

3.3.1 Initialization 


The simulation structure follows the diagram shown in Fig 
ure 3-1. The b^isic iteration cycle of the simulation was 0.05 s, and 
the actuator dynamics were integrated at a 0.01 s rate. All the simu- 
lation runs were made using the same initial conditions. The shearer 
was positioned such that drum No. 1 was cutting the top, with the cut 
height initially set to a perfect cut, that is, the interface height 
minus the bias. Each mine interface profile was modified so that the 
first 30 ft were level. This allowed easy initialization, since a was 
zero in this level section. 


Fixing the front drum with a = 0, fixes the position of 
the shearer chassis and the rear drum in the mine. The initial values 
of Xy and (Xg, Zg) are then computed as follows: 


Given initial Xp-j, Zp-| and Z|, 
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thus 


X] ® Xpi - a c/2 ~ L cos 3i 
X2 = + C 


H ' 


The rear drum is then set to a perfect cut position, given Z 

-1 Z-J + b + Zrjp 

32 = - sin-^ (- 1 — 52 .) 


D2 


The initial values of Zqi, Z-|, and Zq^ are obtained from the mine profile 
data. 

The cut surfaces were initialized as discussed in Section 
3.2.7. The last cuts were assumed to be perfect, and the track formed 
from this perfect cut. The interface data is stored at 5 ft intervals 
and therefore the initial last cuts and track replicate the interface 
and this data spacing. Present cut data are stored at 3 in. intervals 
uniformly across the face. Therefor, on repeated passes, the last cut 
data has this resolution. 


3.3.2 Cut Surface Determination 


The determination of the cut surface, the surface which 
remains after the shearing drum has passed, is critical to the simulation 
of VCS operation and performance evaluation. Since CID sensors ride 
on this surface, and measure the coal thickness from the cut surface 
to the interface, the cut surface simulation affects CIO, and, therefore, 
VCS performance. Evaluation of cutting performance, i.e., rms errors, 
coal taken and left in error, is also affected by the correctness of 
the cut surface representation. 
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The cut surface algorithm starts by storing the Z drum 
coordinates at 3 in. intervals along the face. To find the cut height 
at some location x, the drum positions on either side are checked. 

Figure 3-39 shows this procedure. The cut at point x re- 
sulting from drum position N is the point where the drum radius r inter- 
sects the line X = x. The N'th point is determined by the nearest 3 
in. interval to x, and all the points in the range N - 7 to N + 7 are 
checked for their contribution at x. The cut at x is taken to be the 
highest radius intersection at that point. In Figure 3-39, this would 
be the cut from N + 1. The cut 1 coordinate is given by 

I » Zp(N + 1) + \jr^ - (x - Xjj(N + 1)/' 

3.3.3 Performance Evaluation 


A number of quantities are computed for the performance 
evaluation of the VCS. Figure 3-49 presents an example of the perform- 
ance evaluation computer output. The basic measure used is the rms of 
the deviation of the actual cut from the desired bias thickness. The 
average of the absolute value of this error is also computed, along with 
volumetric coal errors. The excess coal left is that volume of coal 
which is thicker than the desired bias, and the excess coal taken is 
the volume of coal which is thinner than desired. The volume of rock 
taken is also computed. The percentages shown after each of these quan- 
tities represent the percent of the total volume of material taken, shown 
at the top of the printout. The percent time in rock is the percentage 
of the total cutting time that the drum was in rock. 
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3.4 BASELINE SYSTEM PERFORMANCE 


3.4.1 Philosophy and Specification of Baseline System 

The baseline system involves specifying a system concept 
which will solve the required problem and provide a baseline or standard 
against which to compare the results of parameterization studies and 
performance of other system concepts. The baseline is generally the 
simplest system which will solve the problem. To this end, a baseline 
concept was specified for the VCS. The baseline sensors include a nucle- 
onic CID for a proportional measurement and a sensitized pick for a rock 
presence sensor. The last cut follower is included to maintain the cut 
within a given distance of the last cut. For these studies, this dis- 
tance was chosen to be +2 in. Finally, to cut bottom coal, drum slaving 
is used on the bottom drum. 

3.4.2 Performance of Baseline System 

3.4.2. 1 Performance as a Function of Sensor Location 


Figure 3-7 shows the geometry of the CID loca- 
tion. The CID is located at distance Kj^ behind the shearing drum. Table 
3-3 shows the effect of varying the CID location. It can be seen that 
locating the CID closer to the drum results in better cutting perform- 
ance. As Kp is increased, the delay is larger and performance degrades. 
Figures 3-40 through 3-42 show that computer runs were for the 0.75, 

1.5, and 2.5 ft displacements. It is therefore desirable to locate the 
CIO as close as possible to the cutting drum, and with current designs, 
the minimum Kp is 2.5 ft. 
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Table 3-3. Effects of CID Location 
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Figure 3-40. VCS Performance with 0.75 Feet CIO Displacement 
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With 2.5 Feet CIO Displacement 



3. 4. 2. 2 


Performance as a Function of Sensor 
Output Interval Without Noise 


The effects of varying the sensor output 
interval, in this case, for the nucleonic CID, are shown in Table 3-4. 

A measurement is outputted only at the intervals shown in the table. 
Increasing this Interval increases the delay from measurement to 
measurement and also gives a more smoothed measurement, in that the coal 
depth is averaged for a longer distance. It can be seen from Table 3-4 
that performance degrades significantly as the output interval is 
increased, showing that additional delay and smoothing adversely affects 
performance. System cutting performance with a 1, 5, and 10 s delay is 
illustrated in Figures 3-43 through 3-45, respectively^ 

3. 4. 2. 3 Performance as a Function of Sensor 

Output Interval With Noise 

Two effects are evident when the sensor output 
interval is varied with sensor noise. As seen in Section 3. 4. 2. 2, the 
VCS system perfonnance improves v/ilh reduction in the sensor output 
interval. However, the sensor noise will increase as the output inter- 
val is shortened. This is due to the smaller number of counts obtained 
in each interval thus resulting in greater dispersion of these counts. 
Since the two effects are in opposition, a clearly defined performance 
maximum should exist as a function of output interv-3'>. 

Table 3-5 shows this effect for the nucleonic 
CID, and for a 3 and 5 in. natural radiation sensor. The 5 in. sensor 
was simulated by scaling up the 3 in. sensor cal ration curve (Fig- 
ure 3-12) by the ratio of the sensor areas. It can a seen in Table 3-5 
that the nucleonic CID exhibits the best performance at 0.25 s output 
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Table 3-4. Effects of CIO Output Interval 
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Figure 3-44. VCS Performance with 5 Sec. CID Output Interval 
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Figure 3-45. VCS Performance with 10 Sec. CID Output Interval 




able 3-5. 



































Interval. The natural radiation sensors have their best performance at 
2 and 1 s for the 3 and 5 in. sensors, respectively. The nucleonic CIO 
Is better over all, having an rms error of 1.33 in. as opposed to an rms 
error of 2.09 In. for the 5 in. natural radiation sensor. Figures 3-46 
through 3-50 show the nucleonic sensor at 0.25, 1, and 2 s output inter- 
val, and the 3 and 5 in. natural radiation sensors at their maximum per- 
formance points. 


3. 4. 2. 4 VCS Performance in Different Mines 


The statistical data shown in Table 3-2 and 
discussed In Section 3.2.7 show that, statistically, the mines surveyed 
are quite similar. Several of the previous runs were made on both the 
Robinson Run and York Canyon mines. A comparison of results between the 
two mines is shown in Table 3-6. It can be seen that the VCS performance 
in the cases shown is very similar. The differences are due, essentially, 
to the differences between the mines, as shown in Table 3-2. The con- 
clusion Is therefore that VCS performance is affected very little by the 
mine in which It operates. However, mines which are statistically 
rougher will show a* slightly worse performance than will smoother mines. 
Figures 3-51 and 3-52 show VCS performance with a 2.5 ft CIO displace- 
ment and a 1 s output interval for the York Canyon and Robinson Run mines, 
respectively. 


3. 4. 2. 5 Effects of Nonlinear CID Calibration Curve 

The calibration curves for both types of propor- 
tional CIDs are not linear in their relationship between coal depth and 
number of counts received. For both sensors, the curves flatten out at 
larger coal depths. Variation in the number of counts received in a 
given Interval at the larger depths has a greater effect and produces 
larger depth errors than does the same variation at smaller coal depths. 
This effect is shown in Table 3-7. It can be seen that at a 1 s output 
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Figure 3-48. Performance with Nucleonic Sensor at 2.0 Sec. Output Interval with Noise 
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Figure 3-49. Pe»'formance with 3 In. Natural Radiation Sensor at 
a 2 Sec. Output Interval with Noise 
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Figure 3-49. Performance with 3 In. Natural Radiation Sensor at 
a 2 Sec. Output Interval with Noise (Concluded) 
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Figure 3-51. VCS Performance in the York Canyon Mine 
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Figure 3-52. VCS Performance in the Robinson Run Mine 






Effects of Nonlinear CIO Calibration Curve 











interval, there is a 61ji increase in rms error by going from a 4 to a 
6 in. bias. The increase in error is not so evident at larger output 
intervals since the dominant effect there is the delay between outputs. 

These results show that reasonable performance 
is obtained up to a 4 to 5 in. bias. Results for larger biases will 
be worse because of the sensor nonlinearities. Figure 3-53 shows the 
VCS performance with a 6 in. bias. 

3. 4. 2. 6 Performance With Various Control Laws 


The performance of various VCS control configura- 
tions is shown in Table 3-8. These control laws were discussed in detail 
in Section 3.2.6. Table 3-8 indicates the rms errors and the percent 
time in rock. Comparisons of the data show that the CID, pick, and LCF 
offer the best performance, showing a significant improvement over the 
CID sensor only. The sensor used here was the nucleonic CID. The pene- 
trometer performance appears satisfactory, but is made to look better 
because of the limiting effects of the LCF. The CID and LCF offer simi- 
lar performance to the CID, pick, and LCF at 2 and 4 in. biases, but 
at 0 and 1 in. biases, they have a larger rms error and spend more time 
in rock without the pick sensor. Figures 3-54 through 3-58 illustrate 
the results of the five control laws for a 2 in. bias. 

Table 3-9 demonstrates the effect on performance 
of introducing errors in the pick sensor data. The CID and pick combi- 
nation is affected somewhat more than the CID, pick and LCF, with the 
latter essentially unaffected by pick errors. It can also be seen in 
Table 3-9 that there is little difference between the performance with 
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rable 3-8. VCS Control Law Comparison 
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Figure 2-54. VCS Performance with CID Control Law Leaving a 2 In. Bias 
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Figure 3-56. VCS Performance with CID S Pick Control Law Leaving a 2 In. Bias 
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Figure 3-58. VCS Perforraance with CIO, Penetrometer & LCF Control Law Leaving a 2 In. Bias 






























80 and 90X correct pick data. Figures 3-59 through 3-62 show the CID, 
pick, and the CID, pick, and LCF performance at a 2 In. bias with 90 
and SOX pick errors, respectively. 

3.4.2. 7 Performance With Discrete CID Control 

The discrete sensors, when used alone, provide 
a means of taking all the coal and cutting to the CSI. The control law 
used for the pick or penetrometer-reflectometer was discussed in 
Section 3.2.6. Results of simulation studies using picks only or a 
penetrometer-reflectometer and an LCF are given in Table 3-10. Per- 
formance with the pick control shows that excellent results are obtained 
with lOOX correct picks. Performance degrades as the percent correct- 
ness decreases, but is still acceptable for values between 80 and 90 
percent. The percent TIR appears very large for this control configura- 
tion, but the percent rock, also shown in Table 3-10, Indicates that 
only a small volume of rock was actually taken. The pick sensors must 
contact rock to produce a change in their output state, and thus be 
able to control. As a result, the pick sensors are in contact with the 
rock a large percentage of the time, but take very little rock. 

The penetrometer-reflectometer control was im- 
plemented in the same manner as the picks. While the cutting error 
does not appear excessive, the performance is unacceptable. Because 
the sensor is mounted behind the drum, and only outputs a discrete measure- 
ment, the VCS operates in a severe limit cycle. The LCF limits this 
behavior, producing what appears to be acceptable performance. However, 
when the penetrometer-reflectometer is used alone, as is also shown 
In Table 3-10, the resulting errors are very large. Figures 3-63 and 3-64 
Illustrate the pick control performance for 100 and 90% correctness while 
Figures 3-65 and 3-66 show corresponding runs with the penetrometer- 
reflectometer and LCF combination. 
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CS Performance of Rock Presence Sensors Only 
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Figure 3-63. VCS Performance with 100% Correct Pick Control 





IMTERFRC 



LONGWflLL FACE (CT1 “ ^eo-cc 

Figure 3-65. VCS Performance with 100% Correct Penetrometer-Reflectometer . LCF Control 
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Figure 3-66. VCS Performance with 90% Correct Penetrometer-Reflectometer , LCF Control 



3. 4. 2, 8 Effects of Actuator Dynamics 


Table 3-11 shows a comparison between simulation 
runs made with and without actuator dynamics, to access the effects of 
these dynamics. It can be seen that the control configurations involving 
proportional CIOs are not significantly affected by the actuator dynamics 
The rms errors increase from 7 to 16?» with proportional sensors. It 
should be remembered that the VCS actuators are operated in a position 
loop which exhibited limit cycle behavior. Therefore, it is seen that 
despite this actuator loop response, the VCS performance is still accept- 
able. 


The only control configuration which was signifi- 
cantly affected by actuator dynamics was the pick control. This is due 
to the lag Introduced by the rate limited hydraulic system. While on 
a percentage basis, the performance is significantly affected by the 
actuator dynamics, the resulting error of 2.06 in. rms is not excessive. 

It may also be possible to achieve some improvement in this error by 
modifying the pick control law gains. These gains were selected to 
yield good performance with no actuator dynamics and may therefore re- 
quire adjustment when used with the actuator. 

3. 4. 2. 9 Effect of Variations in Chassis Speeds 

All of the preceding results were obtained by 
cutting with a chassis speed of 30 ft/min down the face. Results of 
varying the cutting speed for different nucleonic CIO output intervals 
are shown in Table 3-12. It can be seen that at a 0.25 s output interval, 
cutting performance is unaffected by cutting speed. However, at a 1 
s output interval, the effect of cutting speed is more noticeable. Cutting 
speed has a more pronounced effect at the longer output interval because 
the system is more sensitive to ihe distance traveled between measure- 
ments. 


Table 3-11. Effect of Actuator Dynamics on VCS Performance 
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3.5 FILTERING AND PREDICTIVE ALGORITHMS FOR VCS 


The CID trails the drum and thus does not measure coal seam 
thickness at the point the drum is cutting. Hence, it is neces- 
sary to develop techniques to predict the depth of coal at the 
drum so that the drum can be concuanded to the required height 
so as to achieve a desired coal bias. In this section, four 
filtering and predictive algorithms for VCS control, using present 
cut CID data only and using both present and past cut CID data, 
are developed. 

3.5.1 Philosophy Behind the Developed Algorithms - The philos- 
ophy behind the formulated schemes is to use CID measurements 

or estimates of the coal-shale interface (CSI) only in the immediate 
vicinity of the prediction point which is at the cutting drum. 

This is motivated by the fact that there is no common reference 
frame with respect to which the CID measurements of the CSI could 
be referenced With respect to each other over an appreciable distance 
due to the undulations in the present cut surface (PCS) on which 
the CID is riding and due to the floor undulations. If the VCS 
is doing its job then the CID measurements tend to converge to 
a constant value equal to tlie amount of coal (i.e., bias, B") 
to be left and the cut surface profile will be the same as the 
actual CSI. Since the CID rides on this cut surface, it will tend 
to measure a constant coal thickness regardless of how the CSI 
profile varies. Using measurements or estimates in the immediate 
area of the prediction point minimizes the effect of a floating 
reference frame at the expense of not being able to make an optimum 
use of all the CID measurements in the prediction of the CSI at 
the cutting drum. 

The developed algorithms require a priori knowledge of the 
coal-shale interface statistics in terms of the correlation func- 
tion and the variance as well as a description of the sensor noise 
in terms of its mean and variance. The question of any adaptive 
scheme to estimate the required statistics of the coal-shale in- 
terface from the CID measurements is ruled out due to the fact 
that there is no common reference frame, as discussed above. 

3.5.2 Schemes Using Present Cut Data - The two schemes pre- 
sented here use only the present cut CID measurements. The first 
scheme uses the present measurement to obtain an estimate of the 
CSI at the present point. Coupled with the CSI profile character- 
istics, this estimate is used to predict the CSI at the next measure- 
ment point and at the cutting drum. The estimate of the CSI at 

the next point is then combined with the CID measurement at that 
point to obtain an optimum estimate of the CSI at that point. 

Once this is obtained the CSI at the following measurement point 
and at the cutting drum is made and the process is repeated as 


the shearer proceeds down the longwall face. The second echemo 
uses the present measurement and the immediate previous measure-* 
sMnt to obtain an optimum estimate of the CSI at the present point 
end this estimate is used to predict the CSI at the cutting drum. 

In both schemes the estimates are optimal in the sense that the 
variance of the estimation error is a minimum. 

Before outlining the two schemes, the assumptions and the 
notation involved are stated as follows. 

3, 5. 2.1 Assumptions and Notation - The output of the CID 
is defined as the average coal depth taken over the sensor output 
period. Since noise is inherent in any measuring device, the 

CID output at the measurement (z. ) is assumed to be the sum 
of the average coal depth and the CID error in the k measure- 
ment i.e., 

"k 

The measurement errors are assumed to be uncorrelated since the 

measurement error is in no way influenced by the previous 
or future measurements. Further, it is assumed that the error is a 

0 th • 

zero mean process with a variance . The k measurement point 

on the coal-shale interface (CSI) is assumed to be distributed with 

a mean ]1 (same for all k) and a variance a . It is also assumed 

^k 

that the correlation function for the CSI profiles is known. For 
a perfect cut the final cut surface should have the same profile 
as the CSI but B" (bias-the amount of coal to be left) apart from 
the latter. Thus the mean y is nothing but B and during the final 
pass, the CID measurements would be B" if the control system is 
doing a perfect job. The above assumptions are summarized for 
convenience: 

mean of = y = B 

var {xjj,} = E{ (x^^~y)^} = 


correlation coefficient of and 


E{nj^} = 0 


f “ 

i j 

E(n^n.) . j ^ 

(jvar{n^> = , 

i = j 


ECx^n.} = 0 for all i and j (3.1a) 
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It is to be noted that neither the CSI nor the CID measurement 
error is assumed to be necessarily gaussian distributed. 


The following notation is adopted: i^(k+j/k) denotes the pre- 

dieted value of the coal-shale interface at (k+j)^ point (j » 1,2,3,...), 
i.e., of using the best estimate of the CSI at the k*' measure- 
ment point , i.e., . 


x.x. 
J J 


is a measure of the correlation of x^ and Xj and is defined 


as 


"x., ■ - M)} EEPEODUCllilUTY OF 'mB 

. SiSa!,. rAns: IS P.> « 

where \l is the mean of x* and x.. It can be written in terms of 

^ J 

the correlation coefficient p.. of x. and x* as 

i-J ^ J 

a 0 0 0 

x.x. = . X. X. 

t J J-J J 


where 


a and a are the standard deviations of x. and x. respectively, 
x.x# II 

1 J J 


Thus, once the variance and the correlation function of the profiles 
are known, a is known for all x. and x>. 

X.X. I 1 

The prediction and estimation errors are defined as 


x(k+j/k) = x^^j - x(k+j/k) 


3«5.2.2 Scheme I -" Predictor-Corrector Approach ~ The esti- 
mate of Xj^, the actual coal depth at the x measurement point, 
is obtained by using the predicted value of x^ based on the pre- 
vious CID measurement and then correcting the prediction with 

the new information in the current (k*"") CID measurement, i.e., 
with the difference between the actual current observation and 
the best estimate of the current observation, - 5(k/k-l). 

Thus the algorithm is developed in terms of this predictor- 
corrector approach and the estimator is of the form 
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(3.2) 


Xj^ ■ it(k/k“l) f ■" x(k/k“D) 


where is a weighting factor to be determined such that the 

variance of the estimation error, E{(Xj^ - is minimized. 

Similarly, the coal depth at (k+j)*^^ point (j « 1,2,...) is jjre- 
dieted around |J, the mean of the coal-shale interface using Xj^, 

the current best estimate: 

x(k+j/k) - M + otJjNsL - M) (3.3) 


where is to be determined such that the variance of the pre- 

diction error, Ei.(xj^^j - x(k+j/k))^i, is minimized. The estimator 

and the predictor are developed such that they are unbiased, that 

is 


E(xj^)“ U 
E(x(k+j/k))= ]X 

Thus we have an unbiased minimum variance estimator. The estima- 
tion scheme outlined above is illustrated in Figure 3-67. 



Figure 3-67. Block Diagram of Scheme 1 




The optimal value of g^, that is, the one 
the cost function which is the error variance, 
it found using variational calculus approach. 


that minimizes 

V - 

Let 


«k ■ *S * 


(3.4) 


where g|^ is the optimal value of g^^ and then set 
•a 0 

e»0 

and make this result independent of Using (3.2) and (3.4), 


he 


dV 

ae 


V " E{(xj^ - x(k/k“l) - (gg + “ x(k/k-l)))^} 

“ 0 = E{Ag|^(zj^-x(k/k-l))(xj^-x(k/k-l)-g|^(Zj^'x(k/k-l)))} 


e-0 


which results in the optimal value of gj^ of (3.2) as 

E{(zj^ - x(k/k-l))(xj^-x(k/k-l))} 

^ E{(zj^-x(k/k-l ))^} 

Using (3.1) and (3.1a), (3.5) may be written as 


^ _ var{x(k/k-l)> 

gW *» II ■ I... — II. ■ ■ 

var{x(k/k*-l)} + 0^^ 

where 

var{x(k/k-l)} - “ x(k/k-l))^} 

is the error variance in prediction. 


(3.6) 


. , , t >iv rVY OF 
U F-vF' OF IS FOOB 
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In order to determine the error variance, var{x(k/k-l)) , it i 
required to find of 

aKk+j/k-l) « )i + - y) j * 0,1,2,... 

Again, using the variational calculus approach, is found 

such that it minimizes the prediction error variance 
Vj ■ E{(X|^^j - x(k+j/k-l))^> . 

Let 

-<j) - «0)* . eAo. 

“k-1 “ ®k-l * ^“k-1 

where o;^£^*i8 the optimal value of and then set 


!Zi 

' d e 


= 0 


e*=0 


Using (3.7) and (3.8) 

'^i ’ ^ * ^®k-i* ^ 


bv 


1 
be 


“ 0 *= 


e«0 


(j) 


which results in the optimal value of (3.7) as 




Vi 


E{(\_i~y)^'} 


*k-l 


(3.7) 


(3.8) 


(3.9) 


For j “ 0, define 


Now expressions for the numerator and denominator of (3.9) are ob- 
tained as follows: 


Oa 

Vl*k+j 




- E{(Xj^^j-M)(x(k-l/k-2)-M+gj^_^(Zj^_j-x(k-l/k-2)))} using (3.2) 


E((x^^.-y)(a^.2<V2-»^>^8k-l<Vr»^-V2^W'^»^>. ,, „ 

using (3s3) 


“k-2^^"®k-l^ 


0/s 

*k-2’‘k+j 


^ «k-l 


\-l\+j 


Similarly, 

*k-2*k+j ^“3^^ 

and, in general. 


Oa 

\-3*k+j 


^ «k-2 


\-2*k+j 


Oa 

X . x, . . 
1 k+j 


«,• i(l-g ) O'' 

1-1 X. ,x, 

1-1 k+j 


g. a 

X.X, 

1 k+j 


(3.11) 


for i = 1,2, .. . ,(k-l) 
j = 0,1,2, .. . 

Thus the numerator of (3.9) is computed recursively using (3.11). 
It may be pointed out that (3.11) need be computed only for x. and 

*k+j correlation distance and hence the number of 

Qt^'s and g^'s that are to be stored for evaluating (3.11) are less 

than or equal to 2(k-l). 


The denominator of (3.9) can be evaluated as follows: 


» E{(x. - M)^} 


E{«x(k/k-l) - y) + - $(k/k-l))) 2 } (3 2) 


E{(a, 


'k-l^Vl - ^ «k^ - «kVl<Vl - (3. 


3 ) 


* CX,^ 1 + St? (O^ — Qt 0«. 4- n^\ j. 

Vi I" \ 

'®kVl"i^_j " «k“k-l‘’S|^ * ^“k-A''5 , - ^ 


*k-l*k ^ \-l 


Using (3.10) it can be seen that the terms in the braces cancel 
out> resulting in 


.2,^2 


*k ^ *k-l ^ \ \-l\ " 


(3.12) 


The variance of the prediction error may be computed as follows: 
var{x(k/k-l)} = eUx,^ - x(k/k-l))^} 

E{((xj^ - y) - - y))^> using (3.3) 

' • ^“^-‘“vi-k * <-A-^ 


■ ®k-i*^x V using (3.10) 
*k ^ \-l*k 


(3.13) 
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(3.14) 


Thus (3.12), using (3.6) and (3.13) becomes 
8jj var ljc(k/k-l)} 

The algorithm developed above is summarized as follows: 


ALGORITHM I k > 1 

Filtering: 

Xj^ » x(k/k-l) + - x(k/k-D) 

Prediction: 

x(k+j/k-l) =» M + 


j = 0,1,2,. 


*k 


var(x(k/k-l)} 
var{x(k/k-l)} + o 


n 


var(x(k/k-l)> -a - «. 

*k ^ \-l\ 


„(j) _ °ViV3 


k -1 „2 

X, 


k-1 


j “ 0,1,2,.. 


Vi Vi 


OE: 


*^X.X, = ot. ,(l-g.)a/* + g.o . n 1 o 


k+j 


i * 1,2, . .. ,(k-l) 


X 8,, var{x(k/k-l)> 

*k '■ *k-l\ ^ 
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A few remarks are in order. 

a. If *k+j ^®*^*-**®*' than the correlation 

distance of the coal-shale interface then a„ „ 0 since 

\-l*k+j 

Pir 1 ® hence Oa 0, a,, , 0 and 51(k+j/k-l), 

K-i.K+J *k~l*k+j 

the predicted value of j » becomes which is as expected. 

2 

b. If the noise variance, 0, then gj^ 1 and hence 

^k ** *k should be. 

2 

c. If the noise variance, -► «, then gj^ ■+> 1 and hence 

» x(k/k-i) which is true, .s-ince a highly noise corrupted 
observation contains practically no new information. 

d. If the variance of the coal-shale interface is large 

2 

i.e., if a •> “> then gj^ •+• i and hence = Zy^ 

^k 

2 

e. If the variance of the CSI is very small, i.e,, cr -► 0 

then a ->• 0 and hence g,, •+ 0 and a. -*■ 0 and the predicted 
*i*k 

value of is the mean of CSI profile, y and ® y . 

A few sample calculations are given illustrating the developed 
algorithm. 

Sample Calculations : 


Initial measurement : Zq 

2 2 2 

Initial estimate ~ ''' ^ ■r. 

o o U Xq Xq n 


a 


(i) %^i+i /vi+i 


2 

Oa 

X 


0 


a 


Po (j+1) 


a„ o 
^0 


a 


«0 “«0 


( 0 ) 


XqXi 


; a 


( 1 ) 


XoX.2 


, etc . 
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Prediction: 


x(j+l/0) - M + 

e.g. j - 0: x(l/0) - M + «q(xq - M) 

j - 1: x(2/0) - M + a^^^CxQ - p), etc 


varixd/O)} ■ - o«(T 

xi 0 XqX^ 


g 


1 


var{x(l/0)} 

var{x(l/0)} ^ 

n 


1 = x(l/0) + 8j(z^-x(l/0)) 


al 

X, 


% X *1 varlx(l/0)> 

XgXi I 


k ■ 2 


^*i*9a* X ^ X 

1 2+j V2+j ^ l*2+j 






Oa /ak 
*l’*2+j *1 

Oa 

^1*2 

0#\ 

X, 


Heproducibiiity of the 

ORIGINAL P.\GH IS POOR 


Prediction: 

x(j+2/l) = p + a[j^(5^-p) 


e.g. 


x(2/l) » M + ^1^*1 “ 

x(3/l) •* y + - y) where a 


and 

Xj_X 


var{x(2/l)} = - a. 

*2 "^ 1^1 

g - vartx(2/l)} 

^ var{x(2/l)} + 

n 

X2 ~ x(2/l) + g2^^2 ~ 

“ “I'^x X '*' ®2 


k = 3 


X 

2’^3+j 


a. (l-g2)CTA ■*■ So X 

1 3+,i 2 3+j 


where 


*1^3+j 


Vk,^. 

0 3+J 1 3+j 


a. 


(j) 


X 

^2^3+ i 
2 
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ufeato 


.( 0 ) 


*2X3 

X. 


Prediction: 

*(3+j/2) - M + ol^^Cx, - li) 


e.g. 


x(3/2) “ M + ®2^*2 ” 

?c(A/2) * M + t» 2 *^^(x 2 - li), etc. 


where a 


(1) 


a« 


1th. 

a« 

x„ 


and so on. 

For an example of prediction at the cutting drum assume that 
the chassis velocity is 15 ft/min and that the CID output period 
is 10 seconds. Then j = 0 and the prediction at the drum is given 
by, say for k*=3, x(3/2). 

3 . 5 . 2 . 3 Scheme II - Approach Using Present and Immediate 
Previous Measurements - In scheme II the best estimate of Xj^ is 

obtained as a linear combination of the present measurement Zj^ 

and the immediate previous measurement That is, the best 

estimate is sought in the form 

\ - U + «l(*k-l 

where gj^ and g 2 are to be determined such that the error variance 
Vj " E{(x^ - S^)2) 

is minimized. Again using the variational calculus approach, 
let 


- gif + ejAgj 
g2 - 8$ ^ £2^82 

where gj and g| are the optimal value of and and then set 


9£i 


£^^-£2 -0 


&£. 


"r-^’2“0 


V 2 - E{ (X|i'’M“gf<2^..j-M)-g*(2!,^-li)-ejAg^(Zj^_^^-p)-e2Ag2(z^-y))^} 

av. 


ae. 


Ei=£2®0 


0 = EW8i(\.,-|iK*^-|J-E5^(2,,_j-p)-g|(Z|^-|l))J 




d£ 


2 ki«£2- 


.0 “ 0 « E{Ag2(Zj^-li)(Xj^-p-g*(zj^_^-p)-g*(zj^-p))} 


which result in 

gt (a^ + a^) + g*a * a 


W-1 2 Xj^ n X 


Solving for g| and g^, 


* 1 2 

St ~ TT O a 
'’1 D. n X, X, , 

K k k-1 


where = g^ + 

and gu = 0^ - fff ^ 

^ Vl \ Vi\ " \ 
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Thus 




(3.15) 


The CID coal depth measurement at (k+j)*"” point (j « 1,2,...) 

is again predicted around the mean oE the coal-shale interface 
and using Xj^, the current best estimate: 

x(k+j/k) - M + - M) j - 1,2,3,... 

where *»« determined so that the prediction error variance 

Vj ■ E{ (X|^^j“x(k4’j/k)) } is minimized. Again using the variational 
calculus approach, let 

* ^ifk 

then 
^3 ’ 


dC 


0 - E{AY^(\-y)(x^^.-M-Y^j^* (X^-y))} 


e-0 


( 1 ) 

which results in 


E{(x. -m)(x. _^,-y)} % 


hil 

E((\-y)^} 


VScii 


The numerator and denominator are evaluated as follows: 


* 0?0 




using (3.15) 
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idbikkiiiliiyiiy 






ol - E{(x^-y)^} 

k 

* -4 E{ (a^o 


t, K K— 1 


or 0$ “ »4 0^ a^ + -4 6i io^ + (o^ + a^)e } 

*k \ \*k-l ” Dj^ ^ Wl " \ " 

The variance of the prediction error can be computed to be 

var{x(k+j/k)} •= Oa 

k+j W+j 

A summary of the algorithm developed above is given below. 
ALGORITHM II 
Estimation: 

"" 1'’" * 

where 

D, « g, + a ^ 0 ^ 

k *^k n X, . n 
k~l 


and 


2 7 9 o 0 

^ \-l \ ^-I’^k ^ \ 


Prediction; 


x(k+j/k) = p + - M) j = 1,2,3. 


where 


'k 






2 

X,. 
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- JL 


"K ' i "kw" * S“’vi»,'’» * < * 


Note that: 


a. as 0, 

b. as o^ -► 00, \ M 

2 

c. as 0, ^ ^ 

Ic 

'2 

do as CT X- -► z 

^ ’'k 


which are as expected. 

Sample Calculations: 

Initial measurement: 

Initial Estimate: x^ = x. ■= z ; 


*. Prediction; 


*(j/0) « y + Yq^(xq - y) 


j = 1,2,3,. . 


where 


Yq 


(j)= Vj 


4 


k « 1 


A 1 f 

X, - M + K-{a 


Dj^^^xj^'n ^*0 - ^ 6i(*i - M)> 


”"1 “ «1 ^ '^^Xq ^ ^n 

*1 * ‘^x ‘^x ~ °l X * 

*0 ’‘l *0*1 " *1 
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b. Prediction: 


?l(j+l/l) y + ” m) j “ 1,2,3. 


Oa 

where Yfj’- 

0A 

X, 


1 r 2 1 

Ca = — {a a a + g,a } 

VjH Vi Vj+1 " ^ ^i^jn 


0A « {a^ 0^} + ^0^ 0^ + (0^ + 0^)g,} 

Di XqX^ ^ x^,x^ n n 


and so on. 


3, 5. 2. 4 Compariaon of the Two Schemes ~ The first scheme is 
intuitively appealing in the sense that it is similar to Kalman 
filter of modern control theory. Any error due to PCS undulations 
is minimized by this scheme and it can be shown for some special 
cases that it results in less prediction error variance compared 
to scheme II and it is fair to assume that this is true in general 
cases. However, since it is so vitally dependent on all the corre- 
lations within the correlation distance it is to be expected that 
any errors in the correlation function of the CSI profiles result 
in an increased prediction error variance. Scheme II does not 
suffer from this drawback as much as scheme I since it is dependent 
on only the correlations of points which are close together. 
Further, it is computationally simpler and so is easier to imple- 
ment. However, it is more sensitive to PCS undulations as the 
algorithm depends on a previous measurement. 


It may be pointed out that when any of the measurements is 
zero, say due to cutting the rock, then the drum should be lowered 
by B" and the control algorithms are st.srted afresh (i.e., set 
k ® 0) from the point where the measurement is found to be zero. 


It is instructive to compare the developed schemes with the 
baseline scheme known as the measure and cut scheme wherein the 
present CID measurement at Xj^ is taken as the depth of coal at 


the drum (which is at >:.). The prediction error variance in the 

® 2 2 

measure and cut scheme is equal to 2(l-p,,j)0 0_ • Table 

tCQ X 

k d 


n 


3-13 illustrates that the prediction error variances of schemes I 
and II are less than the measure and cut scheme. It can easily 
be shown for both schemes I and II that 

















- 


ev (predicf.ion) < ev (estimate and cut) < ev (measure and cut) 

where ev (.) denotes the error variance for (•) ai^d estimate and 
cut: refers to taking the optimum estimate of CSI at CID (i.e., 

j^) as the depth of coal at the drum. 

3.5.3 Schemes Using Present and Past Cut Data - In scheme III 
all the available (present or past cut) CID measurements that 
fall within a semicircle of specified tadius (say, two measurement 
distance) with the present CID measurement point as the center, 
are employed to obtain an optimum estimate of the CSI at the pre- 
seiit measurement point. The estimate thus computed is then used 
to predict the coal depth at the drum. Sch* me IV uses estimates 
obtained during the present and the immediate previous passes, 
zeroing in on coal (■hickness at the drum from two directions - 
one in the direction of the cut and one into the coal seam face. 

Tlie present measurement is coupled with the predicted value at 
the present point to get an optimum estimate of the CSI at that 
point. This estimate, which is along the present cut, together 
with the immediate previous estimate of the CSI along the seam 
face (and which is along tht; immediate previous cut) is used to 
predict the OSl at the next CID measurement point as well as at 
the drum. This two-dimensional version of predictor-corrector 
approach is continued as the shearer proceeds dovm the longwall 
face as well as for subsequent passes. In both the schemes, the 
estimates are optimal in Che sense that the variance of the esti- 
mation/predict ion error is a minimum. 

Before outlining the two developed schemes, the assumptions 
and the notation involved are stated: 

3. 5. 3.1 Assumptions an d Notation - The. CID output at che k*"^ 

» till • 

measurement during tl^e r pass, denoted by z(k,r), is assumed 

to be the sum of the average coal depth, x(k,ir}, and the error 
in the k * measurement, n(k,r), i.e., 


z(k,r) ~ x(k,r) + n(k,r) 


(3.16) 
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The assumptions made are suinmarlr,ed as follows: 
mean of x(k,r) » E{x(k>r)} = |i 

2 2 

variance of x(k,r) « E{(x(k,r)-p) } = ®x(k r) 
mean of n(k,r) " 0 

2 

variance of n(k,r) = 

It is further assumed that the correlation function of the CSX pro- 
files x(k^,rj^) and x(k 2 »r 2 ) is known. It is also assumed that 

the measurement errors are uncorrelated and that the actual CSl 
is uncorrelated with the CID measurement error, that is, 

Elx(kj,rj^)n(k 2 ,r 2 )>=* 0 for all kj, k 2 » r^^, and r 2 

The following notation is adopted: x(k,r) and Xp(k,r) denote, 

respectively, the best esti.raate and the predicted value of the 

CSX at the k^ measurement point during the r pass of the shearer 
down the longwall face. x.(k,r) represents the predicted value 

^ th 

of x.(k,r), the coal depth at the drum from the k measurement 
^ th 

point during the r^ pass (henceforth referred to as the (k,r) 
measurement point). 

O y. \ /L. N is a measure of the correlation of the 
X V ^ 1 > r 2^ / X V iC 2 > ^ 2 

CSX x(kj,rj) and x(k 2 ,r 2 ) and is defined as 

'’x(kj.rj)x<kj,r2) * « (x(kj .r,)-M)(x<k2,r2)-M)} 

Xt can be rewritten in terms of the correlation coefficient 
oCkj.rjHkj.rj) « 

*^x(kj,rj)x(k2,r2) ~ ^(kj,rj)(k2,r2)®x(k2,rj)*^x(k2,r2) 


;i 'j { liiiJ 

IS POOR 


3-139 


The estimation and prediction errors are respectively defined 
as 


x(k,r) =» x(k,r) - x(k,r) 

■5ip(k,r) » x(k,r) - Xp(k,r) 

3 . 5 . 3 . 2 Scheme III - Two-Dimensional Version of Scheme II - 
The best estimate of x(k,r) is obtained from all the CID measure- 
ments within a radixis of m measurements («•) from the (k,r) point. 
This is illustrated in Figure 3-68 for m 


PRKSKHT CUT « 


PAST CUT • 


\ 

\ 

• \ 



/ 

^ « 
\ x(k-l,r-X) x(k,r-l) i ,r-i) 


x(k,r~2) ® 


Figure 3-68. Scheme III 

The estimate x(k,r) ia sought as a linear combination of the measure- 
ments z. and the estimator is tf the form 


x(k,r) = li + 


I. 

i=l 


or in vector form 
ft(k,r) - ]x + 


T 

^'■(kjr)(^(k,r) - ]x) 


(3.17) 
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where Che vector ^(k,r) is determined so Chat Che estimation error 

/V 2 

variance! Vj « E{ (x(k,r)-x(k,r)) }, is minimized. For the case 

m « 2 (see Figure 3-68) Che vector ^(k,r) consists of all Che measure 
ments within the dotted semicircle. The vector has y as its 
components. 

The coal depth at the drum is predicted around y, the mean 
of the CSl, using the current best estimate at (k,r) point ob- 
tained from (3.17) 

Sj(k,r) ■ y + Yj(k,r)(5Kk,r) - y) (3.18) 

where Yj(k,r) is found so that the prediction error variance, 

V2 “ E{(Xj(k,r)-5lj(k,r))^} , is minimized. 

The optimal value of £ of (3.17) is found using variational 
calculus approach. Let 

£(k,r) * £*(k,r) + eA£ 

where jg*(k,r) is the optimal value of £(k,r) and then set 



b e 


and make this result independent of A£. Thus 


5!i 

be 


= 0 * E{A£^(£(k,r)-j£)(x(k,r)-y-£*^(k,r)(z(k,r)-]i))} 


e“0 


which results in the optimal £(k,r) of (3.17) as 

T 

£(k,r) = (E{(£(k,r)-ji)(£(k,r)-£)'’'}3 E{(£-|i)(x(k,r)-y)} (3.19) 


Again, using the variational calculus approach, let 
Yj(k,r) = + eAYj 


tnrprn uClBlfXrY 


OF 

3 POOH 


Thus 


where Yj(k,r) is the value of that minimizes 


!!2 

d e 


so that 


0 “ E{&Y,(x(k,r)-li)(x,(k,t)-Vi-Y^(k,r)(x(k,r)-Vi))} 


e“0 


E((x(k,r)-Vi)(x ,(k,r)-li)} » ^x(k,r)x,(k,r) 

Y^Ck.r) « 2 (3.20) 


E{(^(k,r)-y)n 


i(k,r) 


The numerator and denominator of (3.20) can be computed as follows: 
®8(k,r)Xj(k,r) “ 8{(*(k,t)-M)(Kj(k.r)-n)} 


jg^Ckjt) El(z(k,r)-*}i) (x^(k,r)*-ii)} using (3.17) 


^l(k,r) “ E{(x(k,r)-U)^} 

= ^*^(k,r) E{ (^(k,r)~li)(x(k,r)“U)} iising (3,17) and (3.19) 
The variance of the prediction error can be found to be 

var{x^(k,r)} = E{(x^(k,r)“x^j(k,r) )^} 


1st, 


x,(k,r) 


Yd(k>’^)‘^^(lc,r)x^j(k,r) 
using (3.18) and ( 3 , 20 ) 


(3.23) 


(3.21) 


(3.2 


The developed algorithm is summarized for convenience: 


ALGORITHM III 


Bstination: 

x(k,r) - M + £'’^(k,r)( 2 (k,r)-M) 
Prediction: (at the drum from (k,r) point) 

Xj(k,r) * p + Yj(k,r)(x(k,r)-M) 

Yj(k,r) = ^x(k,r)x^(k,r) I *^x(k,r) 


‘’S(k,r)Kj(k.r) ■ «(z(k,t)^)(*j(k,r)-M)) 

’^x(k,r) ■ £^(k,r) E{ (_z(k,r)*-ji)(x(k,r)-M) } 

£<k,r) » [El(z(k,r)-ji)(z(k,r)-]i)^}] E{(z(k,r)-Ji)(x(k,r)-M)> 

Note that during the first pass, as well as at the first measurement 
of each pass, scheme II is used for estimation and prediction. 

During the subsequent passes one can resort to the cases m - 1 
or m ■ 2. 


^ For the case m = 1, with ^(k,r) = (z(k-l,r), z(k,r), z(k,r- 
1)) , the (noninverted) symmetric matrix that appears in the ex- 
pression for ^(k,r) can be explicitly written as follows: 


E{(z(k,r)-ji) (^(k,r)-jj)^} 


^x(k-l,r)^^n ^x(k-l ,r)x(k,r) *^x(k-l ,r)x(k,r-l) 

2 2 

^x(k-l ,r)x(k,r) ^x(k,r)^^n ^x(k,r)x(k,r-l) 

2 2 

°x(k,r-l)x(k-l ,r) ^x(k,r-l)x(k,r) *^x(k,r-l)^*^n 


Also 


E{(z(k,r)-ji)(x(k,r)-}j)) } 


^x(k-l,r)x(k,r) 

„2 

'^x(k,r) 

*^x(k,r-l)x(k,r) 
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The above two expressions doterraine ^(k,r). Explicit expressions 
for (3.21) and (3.22) can be obtained in & similar manner. 

3 . 5 . 3 . 3 Scheme IV - Two-mmensional Predictor-Corrector 
Approach - The optimum estimate of x(k,r) is obtained from the. 
present GID measurement z(k,r) in ton junction with the predicted 
value at (k,r) point, i.e., Xp(k,r), which is based on x(k-'l,r), 

the best estimate of the previous measurement point along the 
cut, and x(k>^~l)s the best ostimate of (k,r~l) point out the 
seam face from the (k,r) ooint. This is illustrated in Figure 
3-69. 


t'RESEMT C!JT« 


PAST CUT » 


'Ck-ljr) 2(1 
« * - - ■ ■ * • 


5(k,. 1) 





(k-M,r) Point 


Drum Point, ;fj(k,r) 




• 9 

x(kU,r-}) 


Figure 3-69. Sclieme IV 
The estimator is of the form 

ft(k,r) ® S^(k,r)+g(k,r)(z(k,rj-Xp(k,r)) (3.24) 

where the weighting factor g(k,r) is determined so that the vari- 

2 

ance of the estimation error, = E{(x(k,r)-2(k,r)) }, is mini- 
mised. The coal depth at the next CID measurement point (3.25) 
and at the drum (3.26) are predicted around y, the mean of the CSI, 
using the current be.st estimate along the cut, x(k,r), and the. 
immediate previous best estimate along the seam face, 2(k+l,r-l): 

Sp(k+l,r) = y+a(k,r)(5i(k,r)-ji)+3(k+l,r-l)(§(k+l,r~l)-)ji) (3.25) 

S!l^(k,r) “ yi-a^(k,r)(S(k,r)-y)+3^(k+l,r-l)(ii(k+l,r~l)-)i) (3.26) 
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Here a(k,r) and 3(k+l,r-l) are weighting constants to be deter- 
mined so that the prediction error variance, ® E{(x(k+l,r) 

-XpCk-^l ,r))^} , is minimized. Similarly, aj(k,r) and 3j(k+l,r-l) 

are found such that the prediction error variance at the drum, 

2 

Vj ■ E{ (x^(k,r)-Xj(k,r)) }, is minimized. 

The optimal value of g(k,r) of (3.24) is found using variational 
calculus approach. 


Let g(k,r) * g*(k,r)+eAg 

where g*(k,r) is the value of g(k,r) that minimizes V^. Thus 

av. 


ac 


e«o 


a 0 “ E{Ag(z(k,r)-Xp(k,r))(x(k,r)-Xp(k,r)-g*(k,r)(z(k,r)- 

X (k,r))> 

P 


which results in the optimal value of g(k,r) of (3.24) as 
var{x (k,r)} 


g(k,r) * 


var{Xp(k,r)} + 0^^ 


(3.27) 


where var{x (k,r)} = E{(x(k,r)-x (k,r)) } is the variance of the 
P P 

prediction error. 

Equation (3.25) can be written in a form similar to (3.17) and, 
hence, the weighting factors a(k,r) and 3(k+l,r-l) can be computed 
in’ a manner similar to ^(k,r) of (3.17) (see (3.19)) and can be found 
to be 


a(k,r) 


■ 

®x(k,r) 

'^x(k,r)x(k+l ,r-l) 

-1 

*^x(k,r)x(k+l ,r) 

^(k+l,r-l) 


‘^X(k,r)$(k+l,r-D 

„2 

^x(k+l ,r-l) 


*^x(k+l ,r-l)x(k+l ,r) 


which results in 


1 2 

^°'x(k+l,r-l)°^x(k,r)x(k+l ,r)~^x(k,r)x(k+l,r-l)^x(k+l ,r-l)x(k+l,r) 


o(k,r) 


* 4 


and j j 

e(k+l,r-l) « Dr“^(k,r)^Si(k+l,r~l)x(k+l,r)"°'S(k,r)5i(k+l,r-lA(k,r3x(k+l,r)^ 

(3.29) 


where 


^kr ” ‘^St(k,r)^5i(k+l,r~l ) °^$(k,r;5c(k+l,r-l) 


(3.30) 


Expressions for all the correlations that appear in the right hand 
side of equations (3.28) thru (3.30) can easily be obtained 
by using (3.24) and (3.25). These expressions are included in the 
summary of the developed algorithm below. 

It may be pointed out that during the first pass Scheme I 
should be used. Also,' the same scheme is used (laterally) at 
the first raeasuremant point during «ach pass. This is explicitly 
given in the algorithm below. 

ALGORITHM IV 

k > ii ^ i 2, i > 1, i > i 

Estimation at the present measurement point, (k,r)t 

ft(k,r) => Slp(k,r)+g(k,r)(z(k,r)-Xp(k,r)) 

Prediction at the next measurement point, (k+l,r); 

^ (k+l,r) ■= ji+a(k,r)(x(k,r)-|i)+p(k+l,r-l)(x(k+l,r“l)-ji) 

P 

Prediction at the drum from (k,r) point*. 

$^(k,r) » |j+a^(k,r)(5i(k,r)-n)+H^(k+l,r-l)(5i(k+l,r-l)~y) 


g(k,r) » 


var{x (k,r)} 
P 

var{2p(k,r)}+Oj^ 


a(k,r) = 


(oi 


ii(k+l ,r-l)°Si(k,r)x(k+l ,r) *^x(k,r)x(k+l ,r-l)^x(k+l ,r~l)x(k+l ,r) 


0A. 


kr 


,) 


3(k+l,r-l) = 5^(aJ(i^^r)‘^2(k-fl,r-l)x(k+l,r)"'^x(k,r)S(k+l,r-l)^x(k,r)x(k+l,r)^ 
2 2 _ 2 

\r ~ ^!S(k,r)^8(k+l ,r-i) ^ii(k,r)fi(k+l,r-l) 


oL 


S(k+l,r) ’ 


x(k+l,r-l) 


CTa 


*'“"'-r)®<'‘**-'-l>'’S(k.r)S<M.r-t)' 

*®‘k*>.--l)»S(k*l,r-l)x(kkl.r)> 

*g^(k-H,r)a^(lj^j_^^*g2(kn.t)o^ 

/^TkTtrTPTl ITY OF THE 

S(k«.k,.(j.l, - ® 

+(l-g(k+l,r))0(k+l,r-l)CT.... , ,, ,. 

x(k+l,r-l)x(j,i) 

+g(k+l,r)cf X , . , X 

x(k+l,r)x(j,i) 

®5(k+l,r)X(j.i) “ ^l-8<*^+l>r>)(l"g(j.O){a(k,r)a(j-l,i)a^^j^ r)X(j-i i) 

*g'^’i’(l'g'k*l-'))‘“(k-')''*(k,r)x(i,i)-6<'»l.-»OS(kn.r-U«(j,n> 


+g(k+l,r)g( j,i)a 


x(k+l,r)x(j,i) 


var{Xp(k+l,r)} = a^_ 


(k+l j when r = i and vice versa) 


+ 2 a(k,r) 0 (k+l,r-l) 05 ,, ,gx„ , ,, 

x(k,r)x(k+l,r-l) 


“2e(k+l,r-l)a, 


5c(k+l , r-1 )x( k+l , r ^‘^|(k+l , r-1 ) 
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The expressions for a^(k,r) and 3j^k+l,r~l) are similar to 
those of a(k,r) and 3(k+l,r“l) appropriately modified. 

During the first pass (r“l); k ^ 1 

^iCkn.l) * Stp(k+l,l)+g(k+l,l)(»!(k+l,l)-5ip(k+l,l)) 

^^p(k+l,l) “ ii+a(k,l)(^t(k,l)-y) 

a^(k,l) “ a+«^(k,l)(^c(k,l)-u) 

var{» (k+l,l)} 

g(kn.l) ^ 2 

var{XpCk+l , 1) } 
var{Sp(k-l,D) - 

o(k,U • ^SCkillsiKliil 

®^(k,l) 

‘^ftCk.Dx.CkjD 

a^(k,l) 2 

*^'3i(k,l) 

j ^*1, 2, •••» 

^ft(j+l,l)x^(k,l) “ “^j’^^^^“®^j'*‘^*’^^'^x(j,l)x^(k,l)'*'®^^‘^^’^^^x(j+l,l)x^(k,l) 

j « 1, 2, ..., k 

'^l(k+l,D “ “^*^’^^‘^x(k.l)x(k+l,l)'^®^'^‘*‘^’^^ var{Xp(k+l,l)} 

Initial conditions: 

x(l,l) “ x(l,l) “ z(l,l) 

4(1,1) “ 4(1, l)'^^ 
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At the first measurement point of each pass (k • l)j r ^ 1 


1k(l,r+l) ■ Si:p(l,r+l)+g(l,r+l)(z(l,r+l)-Sfp(l,r+l)) 

Skpd.r+l) - n+p(l,r)(1t(l,r)-y ) 

Stjd.r) •* ^i+0j(l,r)(xd,r)-u) 

var{x (l,r+l)} 

g(l,r+l) “ ^ — 2 

var{Xpd,r+l)}-KJj^ 

vardtpd.r+1)} - 0^(1 ,r+l)"P^^’’^^'^'xd,r)x(l,rn) 

P(1,0 - 

'^Sld.r) 

‘^5Kl,r)x .d,r) 

0^(1, r) » 2 

‘^SKld+Dxd.r+Z) * P^^'^^^^“8d.i+l)cJ5i(i^i)x(i,r+2)'^®^^’^'^^^‘^xd,i+l)x(l,r+2) 

i » 1 , 2 , r 

°^Hd,i+l)Xjd,r) * 

X “ Ij 2j •••j IT 

4(1, 1*1) ‘ B(l.r)os(i_^)„(l_,^l)*ga,r*l)var{iij,(l,t*l)) 


Initial, conditions same as for the first pass. 


It is seen from Sections 3. 5. 3. 2 and 3. 5. 3. 3 that the basic dif- 
ferences between the two developed algorithms is that one employs es- 
timates in the iimiediate vicinity of the present measurement to obtain 
an optimum estimate of the coal /shale interface at the present measure- 
ment point (algorithm IV), while the other uses CID measurements to 
accomplish the same goal (algorithm III). Moreover, algorithm III uses 
past cut data that is ahead of the present measurement point (but 
obtained during the previous passes) whereas the framework of scheme IV 
does not permit such usage. 

Analytical comparison of all the developed schemes is difficult 
due to the complexity of scheme IV. However, it is expected that scheme 
IV results in the least prediction error variance and is most sensitive 
to uninodeled errors. 


3.6 VCS PERFORMANCE WITH PREDICTIVE ALGORITHMS 


3.6.1 Predictive Algorith m Simulation 

Two simulations were used to study the performance of the 
predictive algorithms: the VCS simulation previously described, and a 
simplified simulation by which parameter studies could be made. 

In the simplified simulation, a random process was generated 
with known variance and correlation properties. The mine correlation 
studies showed that an exponential correlation model could be used to ap- 
proximate the correlation of the actual mines. Therefore, such a first 
order process was generated and inputted to the predictive algorithms. 
White noise of known variance was added to the correlated noise input, to 
simulate CIO measurement noise. The output of the algorithms was then 
compared with the input to determine both estimation and prediction per- 
formance. This was determined by computing the RMS errors. The RMS of 
the input was computed and used as a measure of unaided system performance. 
The prediction refers to the performance obtained by projecting the esti- 
mated CID measurement ahead to the cutting point. Estimation refers to 
the results obtained by processing the CID measurement to remove measu*"o- 
ment noise, and using this estimate as the raw CID data would be used. 
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3.6.2 Predictive Algorithm Performance 
With Simpi 1/fed Simul atioh 


3 . 6 . 2 . 1 Effect of Errors in Correlation Distance 


Table 3-14 shows the simulation results for a 
number of parameter studies intended to determine the sensitivity of the 
algorithm to unknown statistics. The subscripts ACT and f refer* respec- 
tively, to the actual value and the value used in the filtering process. 
^GSI variance of the coal-shale interface. The entries in 

the table denote the RMS error in inches. The nominal run lised the 
nominal conditions listed in the table. For this run, it was assumed 
that the filter had complete knowledge of the process and measurement 
noise statistics. The nominal run is in agreement with the analysis of 
Section 3.5 where it was determined that prediction has better performance 
than estimation, and that both are better than the measure and cut. 


The correlation distance in the filter was varied 
holding the other parameters fixed. The actual process correlation 
distance was a constant 25 ft, while the correlation parameter in the 
filter was set to 13.2, 25, and 50 ft. The results shown in Table 3-14 
indicate that the algorithms are relatively insensitive to uncertainty 
in knowledge of the correlation distances. 

3. 6.2.2 Effects of Errors in CSI Variance 

Filter knowledge of the .31 variance was also 
studied. The filter CSI variance was set to half and twice the nominal 


Table 3-14. Predictive Algorithms - Effects of Parameter 
Variations Sin^lified Simulation 




value. Again, there is little effect, except that slightly worse per- 
formance was obtained where the filter variance was higher than the 
nominal value. In all cases, the prediction and estimation results were 
better than measure and cut results. 


3. 6. 2. 3 Comparison of Schemes 




Parameter variation results for schemes 1 and 2 
are shown in Table 3-14. It can be seen that scheme 2 is slightly worse 
than scheme 1 in both estimation and prediction. Scheme 2, however, is 
still better than the measure and cut results. Parameter variations yield 
the same conclusion as with scheme 1 namely, that the algorithm is essen- 
tially insensitive to parameter variations, and performs as predicted 
analytically. 


3.6.3 Predictive Algorithm Performance With VCS Simulation 

3.6.3. 1 Effects of Measurement Type 

Two possible procedures can be used to derive the 
measurement to be inputted to the predictive algorithms. One procedure 

is to simply input the CID data directly. However, such a procedure intro- 

* 

duces a modeling error because the filters assume that the only randomness 
in the measurements is the CSI and measurement noise. The CID measurement 
also has randomness due to the roughness of the cut surface. A quantity 
called the total measurement is computed as the sum of the CID measurement 
plus the CID relative position sensor. This has the effect of referencing 
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the measurement to the shearer chassis and eliminating the cut roughness 
error. The total measurement and CID measurement results* shown in rows 
2 and 3 of Table 3-15, are essentially the same as the nominal run with 
the simplified simulation. These results show that the total measurement 
yields better performance than the CID measurements. It should be noted 
that the mine bottom was artificially flattened to remove the effects of 
chassis undulation. The relationship between schemes 1 and 2 is the 
same as previously, in that scheme 1 performs better than scheme 2. 

Figures 3-70 and 3-71 show the performance of scheme 1 when total measure- 
ments are employed for both estimation and prediction. Figures 3-72 and 
3-73 indicate VCS performance with scheme 2 using total measurements for 
both estimation and prediction. 

3. 6. 3. 2 Effects of Bottom Undulation 

Another source of modeling errors is the undu- 
lation of the mine bottom which introduces an additional source of 
randomness which affects algorithm performance. Runs with the actual 
bottom, i.e., with bottom undulations, are shown in rows 4 and 5 of 
Table 3-15. It can be seen that under these conditions, the CIO 
measurements yield better performance than the total measurements be- 
cause the undulation causes larger errors in the total computation from 
the skid plane than it does in the CID measurement. 


The effect of undulation has changed the per- 
formance of the algorithms relative to each other and relative to the 
measure and cut. It can be seen that estimation now performs better 




ib 
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Table 3-15. Prediction Algorithms - Effects of Modeling 
Errors VCS Simulation 
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ALL ENTRIES IN THIS TABLE ARE RI^S CUT ERROR IN INCHES. 







Figure 3-71. VCS Performance with Scheme 1 Prediction Using Total Measurements and Flat Bottom 
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Figure 3-73. VCS Performance with Scheme 2 Prediction Using Total Measurements and Flat Bottom 



than prediction for both schemes, and scheme 2 now performs somewhat 
better than scheme 1. A larger error is made attf.'inpting prediction 
than in not attempting prediction because of the untnodeled error due to 
bottom undulation. Scheme 2 performs better than scheme 1 bacause the 
latter is more sensitive to unniodeled errors. Finally, it should be 
noted that both schemes perform worse than measure and cut in the presence 
of bottom undulations. Figures 3-74 through 3-81 present times histories 
of the prediction algorithms runs shown in rows 4 and 5 of Table 3-15. 

3. 6. 3. 3 Effects of Relative Slope Between Top and Bottom 

The prediction algorithms assume that they are 
operating on the sum of zero mean, stationary random processes. If the 
mine top is not parallel to the bottom, that is, if it has a trend or 
slope relative to the bottom, the measurements appear nonstationary, and 
this affects performance. Rows 6 and 7 of Table 3-15 show this effect 
with total and CIO measurements. The CIO measurements are better than 
the total measurements and, again, scheme 2 performs better than scheme 
1. Both schemes remain worse than measure and cut, and larger errors 
are made in prediction than in estimation. 

3. 6. 3. 4 Effects of Measurement Noise Model 

The simulation results presented in rows 1 
through 7 of Table 3-15 were obtained using normally distributed mea- 
surement noise. In fact, the measurement is not normal, but involves 
the nonlinear CIO calibration curve discussed in Section 3. 2. 4.1. Rows 
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8 and 9 of Table 3-15 show the impact on performance of the actual CID 
measurement noise. The relationships between the schemes remains the 
same, but the errors are approximately 15% larger. Row 9 therefore 
represents the performance that might be expected in an actual mine. 
Figures 3-82 through 3-85 show estimation and prediction results for 
schemes 1 and 2, respectively. 

3. 6. 3. 5 Conclusions and Recommendations 

Under ideal conditions, where the measurement 
processes are exactly as modeled, the prediction algorithms can achieve 
a 21% improvement over the use of unfiltered CID data. The simplified 
simulation also showed that precise knowledge of mine statistics was not 
essential to performance. However, when the algorithms are used in a 
realistic environment, with realistic measurement noise, their perfor- 
mance is worse than that obtained with unfiltered CID data. This is due 
to the presence of unmodeled noise processes in the measurements to be 
filtered. The schemes using present and past cut data (algorithms III & 
IV) are not simulated because of the degradation in performance by using 
schemes I and II with realistic mine conditions. 

Also under ideal conditions, scheme 1 performs 
better than scheme 2, but this is reversed in the presence of modeling 
errors. Prediction was also found to be detrimental to performance under 
these conditions because of lack of knowledge of the measurements 
structure. It must therefore be concluded that in an actual mine, pre- 
diction algorithms do not appear to be a viable means of improving system 
performance and hence are not recommended. 
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Figure 3-82. VCS Performance with Scheme 1 Estimation Under Actual Conditions 
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Figure 3-84. VCS Performance with Scheme 2 Estimation Under Actual Conditions 
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3.7 CONTROL OF BOTTOM DRUM FOR BASELINE SYSTEM 


3.7.1 Drum Control Pr( i :edures 

When It Is desired to take all the coal on the bottom, the 
discrete CIDs provide a means by which this can be accomplished. Per- 
formance with these sensors is discussed in Section 3.4.2. 7, where it 
was shown that the sen$iti 2 ed pick can successfully control cutting at 
the CSI. When it is desired to leave bottom coal, an alternate control 
procedure must be used. Since it is Impractical to mount a proportional 
CIU on the bottom drum, locking and drum slaving procedures were 
developed and investigated. 


The drum locking consists of simply setting the drum at a 
desired location and applying no further control. The drum slaving 
procedures actively control the rear drum to cut a constant seam height 
or a constant distance from the top last cut. The sensor configurations 
used in this study are shown in Figure 3-86. Drum slaving was done by 
computing the distance h indicated, where 


h = + Ji sin 3^ + r 


This represents the interface height above the front ranging arm pivot. 
If the desired seam height is H, the ranging arm angle of the trailing 
drum is given by 


0 


2 



h + r + Cg - H 
_ 


‘'I’ 


^ I il iiUu 

IS POOR 
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aving Sensor 



where Cjj Is the desired coal bias on the bottom. A delayed slaving 
procedure was also used» where a measurement h was taken and saved until 
the trailing drum reached the measurement point. This delayed data was 
then used to control the drum. 


To control the rear drum to the top last cut, the height 
h|_ (Figure 3-86) was computed as 

h|_ = + Jt sin 3-| + r 


and 



\ ^ ^ ^ 
I 


These measurements were used in the delayed mode as described above. 

This same LCF procedure was used without the delaying, using an LCF 
mounted in the middle of the shearer. Finally, a rear drum top follower 
was used, which provided a direct measurement of the bottom drum 
location relative to the top present cut, 


3.7.2 Performance with Bottom Drum Locked 

The results of drum slaving and locking studies are 
presented in Table 3-16. Rows 1 and 4 show the locking results with a 
flat mine bottom and a regular or undulating bottom. These studies were 
performed on the York Canyon mine where the top and bottom are uncor- 
related, and the Robinson Run mine where the top and bottom are highly 
correlated. The correlation of top and bottom has a significant effect 
on the performance of slaving. 


11EPKODIICIBH.ITY OF THE 

oricuinAl page is poor 



























When the bottom is flat, the RMS error with the drum 
locked Is the RMS of the CSI. The desired seam cut height is 5.27 and 
7.41 ft in the York Canyon and Robinson Run mines respectively. The 
results show that with the drum locked, the mean seam heights are approx- 
imately 1 in. in error, with a standard deviation of approximately 2 
in. The cutting performance is similar for the two mines. When tne mine 
bottom is undulating, locking performs worse in the uncorrelated mine 
than in the correlated mine. Performance in the Robinson Run mine is 
essentially the same with undulating or flat floor except that with 
undulation, the height standard deviation is larger. Performance in the 
York Canyon mine is worse with undulation than with the flat track. 

Mean seam height and standard deviation are also worse. Figures 3-87 
through 3-90 are simulation runs with both mines, with the flat and 
undulating bottom. 


3*7.3 Performance with bottom Drum Slaved 


The system results with drum slaving are shown in rows 2 
and 5 of Table 3-16. In the York Canyon mine, slaving produces worse 
results than with locking. The differences between locking and slaving 
are much more apparent with the regular mine bottom. In the Robinson 
Run or correlated mine, slaving yields similar results to locking with 
a flat bottom, but worse results with regular bottom. Slaving performs 
better when the top and bottom are correlated, but in any case, all 
results are worse than when the drum is locked. 
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Figure 3-87. VCS Performance in Robinson Run Mine with Flat Bottom and Rear Drum Locked 
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F1gu.-e 3-89. VCS Perfomance In Robinson Run Mine with Undulating Bottom and Drum Locked 
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Figure 3-90. VCS Performance in York Canyon Mine with Undulating Bottom and Drum Locked 



Slaving yields worse performance than locking because it 
tends to accentuate floor undulations. If the shearer chassis pitches 
up for instance, the drum cuts into the top. The VCS corrects for this 
by lowering the drum. The slaving will then lower the rear drum to main 
tain a constant height, but the pitching motion has already caused the 
rear drum to cut into the bottom, and this is now exagerated by the 
slaving. As a result, this type of slaving yields poor results. 

Figures 3-91 and 3-92 show slaving with an undulating bottom for the 
two mines. 


3.7.4 Performance with Bottom Drum Delayed Slaving 

The delayed slaving procedure uses height measurements 
which were made earlier by the front drum to control the rear drum. 

This has the effect of reducing the problems caused by floor undula- 
tions. With the flat bottom, the delayed slaving yields very good 
results with both mines, performing better in the Robinson Run mine 
because of the correlation of top and bottom. With the regular bottom, 
delayed slaving is better than simple slaving, but still does not 
perform as well as locking. The delayed slaving still accentuates the 
effects of bottom undulations. Time histories of delayed slaving with 
bottom undulations are shown in Figures 3-93 and 3-94. 

3.7.5 Performance With Rear Drum Top Follower 

Of the systems shown in Table 3-16, the top follower pro- 
vides the best slaving performance. The average seam heights cut in the 
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Figure 3-91. VCS Performance in Robinson Run Mine with Drum Slavi 
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Figure 3-94. VCS Performence in York Canyon Mine with Delayed Slaving and Floor Undulati 



two mines were the desired heights with standard deviations of 0.25 in. 
Performance is better in the Robinson Run mine because of the correla- 
tion between the top and bottom. The reason for the good performance 
with this sytem is that the top follower provides a measure of the 
effects of chassis undulation. By directly sensing the drum position 
relative to the mine top, and the errors due to undulation are largely 
eliminated. Figures 3-95 and 3-96 show performance with the top 
follower in the Robinson Run an<l York Canyon Mines, respectively. 

3 . ■/ . 6 Performance with .slaving and Sensitized Pick 

Table S-l? shows the results of various other slaving 
techniques discussed in Section 3.7.1. In all of these schemes, the 
rear drum was instrumented with a sensitized pick. The table shows 
that slaving to the top LCF does not yield good performance. This 
procedure was intended to provide a means of limiting the change in 
height from cut to cut on the bottom, but it suffers from the same 
problem as the other slaving procedures. 

Performance with the middle LCF appears somewhat better. 
The use of a chassis mounted sensor was studied as a possible means of 
eliminating the LCF on each drum. While slaving performance was ade- 
quate, performance of the top drum cut was degraded from 1.12 in. RMS to 
3.26 in. RMS. This degradation is due to the fact that the LCF sensor 
has the highest control authority, and therefore errors introduced by 
displacing the sensor to the center of the chassis have a pronounced 
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Figure 3-96. VCS Performance with Top Follower In York Canyon Mine 
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effect. Finally, the rear drum top follower control was used with a 
pick on the bottom drum. It can be seen that use of the pick produces 
only a slight change in performance. Figures 3-97 and 3-98 show the 
results with the front LCF and middle LCF slaving. 

3.8 CONCLUSIONS AND RECOMMEND/\TIONS 

3.8.1 Recommended Sensor Configurations 

Based on the previous performance studies, a recommended 
set of sensors can be configured. Table 3-18 lists the sensors and 
their characteristics that, from a performance standpoint, should be 
used in the VCS. For a given mining application, it may not be 
necessary to use all of these sensors. Preferred combinations of these 
sensors will be given in Section 3.8.2. 

Recommended sensors for the top drum are the nucleonic 
CID, LCF, and sensitized pick. The nucleonic CID output interval was 
chosen to be 0.25 s because this represents the minimum error point and 
also minimized the effects of cutting speed variations. The CID should 
be located as close as possible to the shearing drum, indicated here as 
2.5 ft. The LCF is assumed to be a continuous sensor with a range 
greater than the last cut distance to be maintained. Finally, two picks 
yield good performance, with at least 80 to 90% accuracy. The bottom 
drum sensors include the sensitized pick, and the rear drum top follower 
for slaving. This sensor should also have a continuous output. 
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Figure 3-98. VCS Drum Slaving Performance with Middle LCF 


Table 3-18. Reconinended VCS Sensor Configuration 



RANGE SUFFICIENT TO FOLLOW TOP 


Table 3-19 shows results with a sensor configuration 
selected from Table 3-18. The sensor configuration was a CID, LCF, and 
pick on the top and bottom slaving with the rear drum top follower. 

These represent results from the second pass of a two-pass simulation run. 
These runs are representative of typical system performance in an actual 
mine. Figures 3-99 and 3-100 show performance with the CID sensors on 
top and bottom, cutting a 4 in. bias. 

3.8.2 Recommended VCS Systems For Various Mine Applications 

The configuration of sensors recommended for various 
combinations of top and bottom mining requirements is shown in Table 
3-20. These configurations allow the system to be tailored to the 
individual mine in which it will operate. 
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Table 3-19. Results with Reconmended VCS Sensor Configuration 























INTERFACE 



."NGMflLL FACE 1-' 

Figure 3-99. VCS Performance In Robinson Run Mine with Recoraraended CID, Pick and LCF 
Configuration Leaving a 4 In. Bias (Continued) 
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Figure 3-99. VCS Performance In Robinson Run Mine with Recommended CID, Pick and LCF 
Configuration Leaving a 4 In. Bias (Continued) 
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Figure 3-99, «CS Perfor^ece In Robinson Run Mine with Reco-ended CIO. 

Pick and LCF Configuration Leaving 4 In. Bias (Concluded) 
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Figure 3-100. VCS Performance In York Canyon Mine with Reconriended CID, Pick and LCF 
Configuration Leaving a 4 In. Bias (Continued) 
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Table 3-20. Recommended Sensor Configuration for Various Applications 
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4. YAW ALIGNMENT SYSTEM 


4.1 FUNCTIONAL DESCRIPTION 


In lonqwalT mining after coal is cut along the face, the conveyer 
Is pushed forward in preparation for the next cut by the roof support 
system. A typical sequence of operations is shown in Figures 4-1 
through 4-4. In Figure 4-1 the conveyer is straight, the roof support 
shields have been pulled up and tho shearer is at the headgate. The 
tailgate end of the conveyer is then advanced in preparation for the 
next cut as shown in Figure 4-2. Figure 4-3 shows the shearer at the 
tailgate end of the conveyer after it has cut the coal face. Finally, 
in Figure 4-4 the headgate end of the conveyer is pushed forward and the 
roof supports have been advanced along the tailgate end of the conveyer. 
The shearer can now travel from the tailgate to the headgate making a 
cleanup cut at the tailgate end of the conveyer and completing the cut 
at the headgate end. After the cut, the remaining roof supports are 
advanced. In this position the conveyer, shearer and roof supports are 
back in their original positions ready for another cycle of operation. 


In the manual mode of operation, the miners push the conveyer a 
full stroke and then after cutting, advance the roof supports. Care 
must be taken to straighten the conveyer after the roof supports 
advance. This is accomplished by the miners sighting down the conveyer 
and alternately adjusting the conveyer and roof supports so they are 
straight. If care is not taken, the conveyer will rapidly get out of 
alignment resulting in the roof supports being unable to advance pro- 
perly. 
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In an automatic system, measurements are made to determine the con- 
veyer yaw alignment. Any misalignment is corrected by pushing the 
conveyer different amounts along the face. 

The measurements can be taken when the system is in the sequence of 
operations as showrt in Figure 4-4 The conveyer is relatively straight 
at this time and measurements can be made as the shearer travels from 
the tailgate to the headgate. 

Two types of measurement devices are proposed. One Is an angle 
cart system that measures the relative angle between two conveyer 
sections. The other Is a directional gyro which measures the angular 
orientation of a conveyer section with respect to a directional ref- 
erence . 

Detailed descriptions of these measurement systems along with their 
performance evaluations are presented in tie following sections. 


4.2 DESCRIPTION OF THE ANGLE CART MEASURE -lENT SYSTEM 

The angle cart measurement device consists basically of two 
resolvers as shown in Figure 4-5. One resolver is lined up with one 
conveyer section while the other is lined up with an adjacent section. 
The sum of the resolver angles and 6g gives the angle between the 
conveyer sections. As the angle cart moves along two conveyer sections, 
many measurements (200) can be made and then averaged. Hence 

200 
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Figure 4-5. Angle Cart Measurement Device 
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If the standard deviation of the error in each resolver measurement is 
then the standard deviation of the angle e is 


0 


m 


n 



“ 0.1 

r r 


The bias b In the angle cart measurements can be detemn'ned and re- 
moved from the measurements by computing e when the angle cart is on a 
single conveyer section. Hence, if there are N conveyer sections, the 
bias Is computed by 


1 

N 


N 

5; 

k=1 


where is the angle cart output when It Is on the conveyer section. 
This bias is removed by 

A 

e a e - b 

A 

where e is the corrected angle between the conveyor sections. The 
standard deviation of the error ii b is given by 

4.2.1 Basic Angle Cart Mea< . urement System 

The angle cart measur's the angle e between two conveyer 
sections. Knowing all thee angles, the location of the conveyer end 
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points, and the length L of each conveyer section allows one to deter- 
mine the YAW profile of the conveyer. From the geometry shown In Figure 
4.6 the angle cart measures the angles Op through (assunlng N eon 
veyer sections). The location of Y^ and Y|^ can be measured with respect 
to surveyor stakes. The computed Y coordinates Y^ of the end points of 
each conveyer section are found as follows 


Y^ • L sin e.| + Yq 


A A 


A A 


Yg » Y^ + L sin (e^ + 0p) 


Y^ « Y^ , + L sin ) e . 
n n - 1 1 

1=1 


1=1 


A 

where are the measured angles. Since all the angles are measured 
except 0^, It must be determined 'n anothe” manner. The last equation 
above can be written 


A 



N 


r— \ 

•-E 





Note: 



Figure 4-6. Conveyor Geometry 




V 
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where y^- Is the angular orientation of the conveyer section with 
respect to the horizontal. Ass uni ng small angles 


N N 1 

Y «Y +Ly +Ly + ... +Ly = Y +Lr^y •L^(r^e)+Y 
No I 2 No i Lj S 0 

1*1 1*1 j-1 


Assunlng that 0 yields 


1=2 


The value of can then be cotnpuied by 


* ''n ■ ''n 

®I Si— 


As a result, the yaw profile Is determined by 
n 

A r— 

1=1 

4.2.2 Derivation of Weighting and Bias Estimation System ' 

If the angles ©| and can be measured, then a method of 
weighting and bias estimation can be used to determine the yaw profile. 
The equations are derived as follows: 


E 

j”i 


+ V„ n= 1.2.3 N 




4 




Using and the angle e-| plus the angles 62 through©|^, the 
Y coordinate of the point is computed by 

N/<„\ 

''n-'- E ''0 

i«i yj«i j » 

A ^ •4' 

With angle measurement errors Yj^ will not equal Y|^. Therefore a bias b 
on all the angle measurements is computed so the end points match up. 

r\ 

t+ _ 2 

^ Nffi + iy L 

The corrected angle measurements are found by 

A ^ A A ^ 

e j = “ b 

and the yaw profile is computed using those angles. 

Using Yj^| and the angle 0 |^^,i plus the other angles ©2 through 
the Y“COordinate of the 0^^ point is computed. 
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A 

Again a bias b* Is computed 


NCN -ni 


Y - Y 
0 0 


A1so» the corrected angle measurements are found by 

AAA 

e ![ » e ^ b“ 

and the yaw profile by 

N-n / N+1 

^ r— ^ / V*— > ^ 

''n-'-E E 

1=1 yj=N-i+2 

As a result, there are two Y-coordinate values for each point depending on which 
end of the conveyer was the reference. A weighted combination of the two gives 
a good estimate of the yaw profile. Hence 
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A 




^n “ ^n '^n ^n 

where a^, is the weighting function. It is derived as follows: 


The one-sided estimates can be written 

n i 


i»l j=l 
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I 4 


‘iwi ri£~"'t iiir' 


n 1 


E 


i«l j«l 


j*1 


and 




N-n N+1 


1=1 j=N-1+2 


N-n N+1 

\*<-L E 

i=l j=N-i+2 


where 


' ''n * *2 


n 1 


N-n N+1 
1=1 j=l 1=1 j=N-i+2 


is the true ordinate and 
n 


X, . L 'j 
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and 




N+1 


X 2 • L (j-n-])cj 

j«n+2 


are the errors in Yjj and Y’, respectively. 


It is assuned that the angle measurement errors € 2 * •••» 
*N* *N+1 pairwise uncorrelated random variables with zero means and 
common variance It follows that X-| and X£ are uncorrelated and 

have zero means and respective variances: 


v^ = ^ (n-j+1) 


m 


j»l 


and 


Ml 


N+1 


V rv 

/■ ^ . . 


V 2 « Y (j-n-l)^o 
j-n+2 


2 ^ 2 
m 




^ ' I'OOii 


N-n 




7' 1^2 . 2 (N.n)(M-n.1)(2N-2nt1) 

in 6 


Ml 


For any value of a^, is an unbiased estimate of Yp, i.e., 
E(V-V This is true because 
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* (> - *n) 


and 

E(f„-) • V„ * ECXj) - Y„ ♦ 0 

A 

The chosen value of a_ is that which minimizes the variance of Y_. 

n A n 

Differentiating the variance of with respect to a^, setting the 
derivative equal to zero and solving for a^^ gives 

, . ''2 _ (N-n){N-n+l)(2N-2n+l) 

®n v ^2 ’ n( n+T7(2n+l ) + 

4.2.3 Derivation of Angle Cart Measurement Algorithm Assuming that 
a Number of End Conv ey er Sections Cannot be Measured 

Since the shearer is approximately 30 feet long and it rides 
on the conveyer sections. It is possible that not all the angles between 
the conveyer sections can be measured. If this is the case, the control 
algorithm will simply assume that these angles are zero, i.e., that the 
conveyer is straight at the ends. 


4.3 DESCRIPTION OF THE MONTE CARLO TECHNIQUE TO DETERMINE THE EFFECT OF 
ANGLE CART MEASUREMENT ERRORS 

A Monte Carlo simulation was developed to determine the effect of 
angle cart measurement errors (see Figure 4.7). This was accomplished 
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by developing a conveyer yaw geometry consisting of the true conveyer 
angles e^ and the Y-coordinates Y^. The true conveyer angles were 
sampled from a normal distribution with zero mean and standard deviation 

oy, i*e«, 

* normal (0, oj^) 

. where o ^ 1 degree 

The Y-coordinates are then computed by 

1»1 0»1 

Measurement errors were then added to the conveyer angles and using the 
measurement algorithms described in Sections 4.2.1, 4.2.2 and 4.2.3 the 
computed Y-coordinates Y^ were determined. (It was assumed that Y^ and 
Yjj were measured without error.) Next, the yaw profile error 
■ Y^j - Y^ and the RMS of the yaw profile error 



were computed. The RMS computation was made 500 times— each time using 
different noise measurement errors. 

A statistical average of these 500 RMS values is needed for 
an evaluation criterion. The statistical average used is the Mode of the 
RMS and is computed by 


•St 



where ■ 500. 

4.3.1 Baste Measurement System 


\ 


The baste measurement system ts stmulated by assumtng the 
measured conveyer angles are gtven by 

A 

e = e ^ + e ^ + b 

where 

ts the true geometrte eonveyer angle 
* the random error tn eaeh measurement 
and, b » the btas error tn the measurement devtce. 


These errors are stmulated by sampling from normal dtstrtbuttons, t.e., 

♦ 

2 

= normal (0, 



The results of the Monte Carlo simulation are shown tn 
Figure 4-8. As the standard deviation of the measurement error (o^) 
Increases, the mode of the RMS (t^) Increases Itneartly. 


It will be shown later that the maximum allowable value of 
(with no other system errors) ts approximately 0.2 feet. It can be 




seen» therefore, that the maximum allowable Is approximately 0.04 
degrees. 


4.3.2 Performance of Weighting and B<as Estimation System 

The weighting and bias estimation system uses two angle 

/\ A 

measurements (e^ and6|tj^^) that are not used in the basic system. These 
measurements have an additional error— the directional reference error for 
the first and nth conveyor sections. The measured angles are simulated by 

®1 + 

A 

®N+1 * ®N+1 ®N+1 ^ ^ ^e(N+l) 

where 

* normal (0, ) 


•e(N+l) “ 

The errors in the other angles are the same as in the basic measurement 
system, I.e., 


e . = normal (0, 


b 


normal 



The results of the Monte Carlo simulation are shown in 
Figure 4.9. In this figure Mp is plotted versus Oj^ for various values 
of 00. It can be seen that using the criterion that the maximum allow- 
able M|^ Is 0.2 feet, the maximum allowable Og Is 0.3 degrees. If Oq ” 
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0, then a larger can be tolerated than In the basic measurement 
system. 


It Is concluded that the weighting and bias estimation 
algorithm is a viable system and superior to the basic system if the 
angle cart measurement Inaccuracies are greater than 0.04 degrees but 
less than approximately 0.1 degrees. In that case the basic angle cart 
system does not provide sufficient accuracy. However, the error In the 
directional reference must be less than 0.2 degrees. 


4.3.3 Performance of Basic Meas ur ement System When a Number of 

Conveyer Sections Cannot Be Measured 

The basic measurement system measures the angles 62 through 
and then computes the angle e^. If some of the angles at the con- 
veyer ends cannot be measured because the shearer Is In the way, then 
these angles will be assumed zero. In order to check the effect of this 
condition the Monte Carlo simulation was exercised with a nunber of 
conveyer angle measurements at each end of the conveyer arbitrarily set 
to zero. 

The results are shown In Figures 4-10 and 4-11 where the 
Mode of RMS Is plotted versus the number of conveyer angles at each end 
of the conveyer that cannot be measured. Figure 4-10 was obtained using 
O'c * 0.5 degrees while Figure 4-11 Is for oj » 1.0 degrees. The con- 
veyer profile was simulated by sampling the angles between the conveyer 
sections from a normal distribution with standard deviation oj. There- 
fore, foroy * 1.0 degree the conveyer is more crooked (the angles 
between sections are larger) than foro-p = 0.5 degree. By assiining that 
some of the angles at the ends of the conveyer are zero, larger errors 









in the computed yaw profile are expected. Foray ■ 0.5 degrees little 
affect on is indicated for one or two angles that cannot be measured 
at each end of the conveyer. Foroj * 1.0 degrees the effect on is 
small only when one angle at each end of the conveyer is not measured. 
Since the magnitude of oj depends on operational errors such as the roof 
support pull up errors, the effect of not measuring some of the end con- 
veyer angles must be studied with. the yaw advancement simulation. The 
results of this study are described in Section 4.4.1. 


4.4 DESCRIPTION OF THE YAW ADVANCO-IENT SYSTEM SIMULATION USING THE 
ANGLE CART 

As described in Section 4-1, the yaw advancement system consists of 
pushing the conveyer forward, pulling the roof supports up, measuring 
the yaw profile and then pushing the conveyer forward again. This same 
sequence of operations is modeled in the yaw advancement simulation (see 
Figure 4-12). 


Conve.yer Placement 


Each section of the conveyer is pushed forward an amount given by 
the placement command a Y p^ plus a placement error Cp^, i.e.. 




aY 


pn *■ pn 


where 6p^ is sampled from a norma! 
the conveyer is pushed forward aY^ 
ram stroke The equations are 


distribution (0, Op^)‘ amount 
is limited between zero and the full 



'^7 
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If *y„ > v^. 4Y„ . .y^ 


If 4y„ < 0. jy„ . 0 

The conveyer sections are pushed from their midpoints. Defining 
'^mn Y-coordinate of the midpoint of the n^^ conveyer section, 

then the Y-coordinate after pushing is given by 


yk a yk-l + 
mn mn n 

where the superscripts define the sequence before and after pushing. 

In the simulation conveyer section 1 is placed first. Its orienta- 
tion is sampled from a normal distribution (0, oj). Then its end points 
are located; 


''o * ''ml - r ®I 


Yl » y„, sine, 


Again, the advance is limited to Y^, i.e.. 

If y k y k-1 -f- y Y k s y *<“1 + y y k * 2Y - Y^*^ 
0 0 A’ 0 0 A 1 mi 0 

If Yi*^ > + Y^, Y^*^ = and Y^*' « - Y^*^ 


The angle ei is recotnputed 


« sin 


-I »1 - ''o 


Also 

Yi • e, 

I 

All the other sections are then located In sequence by the 
following equations: 


»n • - Vi 


If Vn" > Y^K-V Y^, 


k -1 


• Sin 


4--1 '^n ^n-1 


®n " I'n “ 

If after computing any 


or 


‘n > (®IMX * ^ 


e. < -e„_. 


'n ' ''max' 


then Is restricted to its maximum value, i.e.. 


or 


and 


®n * ®max 


®n “ "®max 


Yn * ^n-1 ^ ®n 
v„ - v„., + L Sin Y„ 

V ■ ("r, * "n-l) 


Roof Support Pull up 


After the conveyer has been pushed forward, the shearer will move 
along the conveyer making its cut. This operation, however, does not 
affect the conveyer yaw profile. The roof supports are then pulled up. 
This operation affects the conveyer profile by pulling the conveyer back 
as the roof supports are pulled forward. This operation is simulated by 
subtracting a pullup error from the Y-coordinate of the midpoint of each 
conveyer section, i.e.. 


y k « Y k“l « I - I 

mn mn ' sn' 

2 

where e is sampled from a norma distribution (0, a ). 

dll d 

A yaw profile after this pul up error is computed again in the same 
manner as for the conveyer placement. 


Measure YAW Profile 

For the basic angle cart measurement system the conveyer angles 62 
through 9|^ are measured. Therefore, for all those angles errors are 
added, i.e., 


+ b 


The yaw profile is then computed as described in Section 4.2.1. 
The equations are 


A 



- L 





For the weighting and bias estimation sj'stem, two additional angles 
are measured. They are and Errors In measuring these angles 

are the same as for the other annles except that they have an additional 
error. The additional error Is n the reference direction. It Is simu- 
lated by 


®1 ' ®1 + + b + Cq, 



ts 0 + e 

Nt) N+1 


^ ‘0(N+D 


where 

®el ®e{N+l) 

are sampled from a normal distribution (0, 



). 


The remaining equations are: 
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Conveyer Placement Commands 


The conveyer placement conmands aYp^ are computed using the con- 
veyer Y-coordinates Y^. It is desired to advance the conveyer as far 
as possible and at the same time straighten it. This is accomplished by 
determining the minimum Y-coordinate, advancing that section a distance 
Y^ and all the others an appropriate amount so that it is straight. The 
conveyer placement commands that will accomplish that task are derived 
from geometry shown in Figure 4.13 




• Yfl + L,-n - 0.5 (Y„ + Y„ 
pn A mm n n- 1 ) 


where is the minimum value of Y^. 


4.4.1 Overall System Perform ci nce Using the Angle Cart 


The overall system performance using the angle cart 
measuring device was detennined by exercising the yaw advancement 
simulation. 


The conveyer initial position was a straight line. Then the 
conveyer was pushed forward (with placement errors, Op) and the roof 
supports were pulled up (with pull up errors, Og). Measurements were 
then made (with measurement errors, o^) and the yaw profile computed. 
Next, conveyer placement commands were computed in preparation for the 
next advancement. Twenty consecutive advancements were made. 



tiing Pldcement Cofimdnds 



In order to determine tolerable measurement errors, a 
placement error Op and a pull up error were selected and then twenty 
advances were made with 0 ^^ » 0. The measurement error was then 
Increased in 0*005 degree steps each time advancing the conveyer twenty 
times. As the measurement error was increased, the conveyer became more 
and more crooked so that for some value of the conveyer cannot 
advance twenty times without commanding a negative placement of the 
conveyer. Such a connand would require the conveyer to be pulled back 
which is Impossible. As a result, the next lower value of was deemed 
the maximum tolerable measurement error. Three maximuns were deter- 
mined for each placement and pull up error and then averaged. The 
results are plotted in Figures 4-14 through 4-17. 

Figure 4-14 shows the maximum tolerable measurement error 
versus the pullup error for the basic angle cart measurement system. 

Two values of the placement error were considered, ffp = 0 and 0.05 feet. 
The region below the curve is considered stable; above the curve the 
region is considered unstable. 

Figures 4-15 through 4-17 show the results for the basic 
angle cart measurement system where it is assumed that some of the 
conveyer angles at each end of the conveyer cannot be measured (J 
defines the number at each end that cannot be measured). It can be seen 
that for J * 1, system performance is degraded only slightly. For J » 

2, however, system performance is degraded by a considerable amount. 

The reason for this large effect on performance is that the conveyer 
pullup errors cause the conveyer angles to be large. Since they are 
assumed zero in the yaw profile computations when they cannot be 
measured, the pullup errors cause a large effect. 
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0.02 FEET 






In Figures 4-18 and 4-19, the performance of the weighting 
end bias estimation system 1s shown. Figure 4-18 shows the performance 
for oTp ■ 0 and Figure 4-19 shows the performance for Op » 0.05 feet. 
These figures show that the reference error (the error In determining 

V 

the reference direction of the end conveyer sections) has a large effect 
on the maximum tolerable measurement error For example, the maximum 
allowable error In the standard deviation of the reference error Is 0.2 
degree. However, If the reference error can be determined more accu- 
rately than 0.2 degree, the weighting and bias estimation system gives a 
better performance than the basic angle cart measurement system. 


4*4.2 Recommended Measurement Algorithm for the Angle Cart 

The recomnended measurement algorithm for the angle cart Is 
the Basic Measurement System. This system uses the measured location of 
the conveyer end points and the angles between all the conveyer sections 
to compute the yaw profile--the Y-coordinates of the ends of the con- 
veyer sections. From the yaw profile, conveyer placement commands are 
generated. The advantage of this system over the weighting and bias 
estimation system Is that a reference direction is not needed. The 
basic system Is therefore simpler c(nd reference direction errors are of 
no concern. 

' The angle cart eight speed resolvers have a basic accuracy 
of 2 arc minutes or 1/30 degrees. Averaging 200 measuranents gives the 
angle cart measurements an accuracy of 0.0033 degrees. This accuracy Is 
an order of magnitude better than the requirement. 
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Using the yaw advancement simulation* ten advances of the 
basic measurement system were simulated and the yaw profile plotted 
after each advance. An example of these plots Is shown In Figure 4-20 
where Op ■ 0.05 feet, ■ 0.1 feet and * 0.03 degree. It can be 
seen that after ten advances, approximately 14 feet of coal was cut. 
Without errors 20 feet of coal would have been cut. Therefore, the 
system efficiency Is 


Figures 4-21 and 4-22 show the efficiency of the basic angl 
cart measurement system as a function of the system errors. 
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igure 4-21. Efficiency of Basic Angle Cart System 




4.5 DESCRIPTION OF THE DIRECTIONAL GYRO MEASUREMENT SYSTEM 


A directional gyro mounted on the shearer or on a conveyor cart can be 
used to measure (with respect to a fixed directional reference) the angular 
orientation of each individual conveyor section. Knowing the location of 
the ends of the conveyor along with the conveyor angles allows the com- 
putation of the conveyor yaw profile. 

Four different directional gyro measurement algorithms have been de- 
vised. The advantage of each depends partly on the relative magnitude of 
the gyro errors. 

A description of the g>fi? error model and the four measurement algo- 
rithms is presented in the fol lowing sections. 

4.5.1 Mathematical Model of Directional Gyro Error 

The directional gyro is used to measure the orientation angle 
of the cart with respect to a fixed reference direction. The gyro measure- 
ment error e(t) at time t is assuned to be additive, i.e., the gyro mea- 
sured angle at time t is equal to the true angle plus c(t). The dynamic 
mathematical model for e(t) is illustrated in the system block diagram of 
Figure 4-23. 


The input u(t) is assumed to be a Gaussian white noise sto- 
chastic process with spectral intensity K, i.e., the u(t) are independent, 
zero mean Gaussian random variables such that 


E[u(ti) u(t2)] = K 6(ti - tg) 

for any tj and t£, where « denotes the Dirac delta function. 
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It is seen from Figure 4-23 that c(t) is represented as the 
sum of four effects: 

(1) An integrated v4iite noise, i.e., "random walk" process 

(2) A stationary process obtained by passing white noise through a 
linear first order lag filter with gain A and time constant t. 

(3) A drifting bias error Cj^t where is the bias rate, assumed to 
be a Gaussian random variable with mean zero and variance 

(4) A random constant bias error e^, representing a fixed error in 

the directional reference, and also assumed Gaussian with mean 

2 

zero ano variance . 

The random variables u(t), t 0, ejj and are taken to be 
statistically independent. From Figure 4-23, e(t) can be expressed as 



For proper mathematical interpretation, the integral above should be under- 
stood to be a Wiener stochastic integral, with the term u(r) dr replaced by 
dw(r), where w(r) denotes a Brownian motion stochastic process with 
variance parameter K. fhe covariance function of e(t) is given by 


ECc(tj)e(t2)3 = 






+ tj t2 + 
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for 0 £ tj £ tg. If tj > tg, then tj and tg must be interchanged in the 
above formula. 


4.5.2 Derivation of the Directional Gyro Measurement Algorithm 
Without Averaging 

The directional gyro measurement algorithm without .averaging 
uses only a single directional gyro measurement of each conveyor section. 
Those measurements are combined with the Y coordinates of the conveyor 
end points to compute the yaw profile. 


n 



k=l 




+ Y, 


✓V 







+ 



where is a weighting function derived to minimize the error in Y^. The 
derivation of b^ follows. 

Let T denote the time for the cart to traverse a conveyor 
sepent. The gyro measurement angle V|^ = yj^ + where y|^ is the true 
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angle, and the error C|^ » e((k - 1)T) for k * 1, .... N, since the gyro 
measurement is made at the beginning of the segment to obtain the si.n»!l1est 
error variance. Thus, 


''n* ■ E f'k " 'k> * ''o “ ''n " ’‘l 

k-1 






where and are the true Y coordinates at points n and N, respectively, 
and 


Hence , 




+ X, 


»r ><2 


The value of is chosen to minimize the variance of Y^, and 
is obtained by setting the derivative of this variance with respect to p^ 
equal to zero and solving for 
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MXiX;) 

E(X/) 

KA (N + 1 + n + ^ n (n - 1)(3N - n -1) + 


o^Vn(N - 1) 


"I 

1 


fi(n - 1) *f nNap 


KA (N + 1 + ^) N + N (N 




1)(2N - 1) + (N - 1)2 + n2«j^2 


The quantities li(XjX 2 ) and E(X 2 ^) were computed by summing the appropriate 
covariances given by the error covariance function in Section 4,5.1 and 
then neglecting small terms. 


4*5.3 Derivation of the Directional Gyro Measurement Algorithm with 
Averaging 

The directional gyro measurement algorithm with averaging 
uses approximately 200 directional gyro measurements of each conveyor sec- 
tion. An average of the 200 measurements is computed and then used as a 
single measurement to compute the yaw profile. The aquations are: 


200 


^An 


- 1 XT' 
“TT L 


'mi 


i=l 


where are the individual measurements and y/^^ is the average measured 
angle of the n^^ conveyor section. Then 


H 
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where is the weighting function derived to minimize the yaw profile 
error for this condition. Its derivation follows* 

Proceeding as in Section 4.5.2, tie error* e|^ is now given by 


kT 

(k = 1, ..., N) 

-1)T 


where the discrete average has been replaced by an integral for computa- 
tional simplicity. Hence, 



*1 =T r 

^ c 


nT 


; ( t ) dt 


■rf 


NT 


:(t)dt 


= 


rNT ^nT 

j j E[e(ti) e{t2)] dtj dtg 

0 0 


n TNT TNT 


J J E[e(tj) eCtg)] dtj dt£ 

^ r\ ^ r\ 


0 ' 0 
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Again, using the error covariance function in Section 4.5.1, performing the 
integrations above and neglecting small terms, it can be shown that 




where 


= k[.T 

1 

- ^ (N + ^ + n)J+ — 


o ^ 

+ N n 


[ 


N^T 


2t 


At 


Q » K I -iy- + AN^ + A(— N - ^ (A + 2t) - y (2N + 


At 

"T 




t2 

T Oi 


Or 


4*5.4 Derivation of the Directional Gyro Algorithm with Drifting 
Bias Compensation 

The directional gyro algorithm with drifting bias compen- 
sation is derived by assuming that the main gyro error is the drifting 
bias. The directional gyro at the beginning of a sequence of measurements 
is nulled with respect to a fixed directional reference. As the gyro pro- 
ceeds at a constant velocity along the conveyor taking measurements, it 
drifts at a constant rate. Hence, the error in yj = Ay; the error in y 2 ® 
2ay; and the error in y^ * oay. As a result, Ay can be computed using the 
measured and computed end points, i.e.. 
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where 



The directional gyro measurements are then corrected by 
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Yn - hay 


and the yaw profile is computed by 



4*5»5 Derivation of the Directional Gyro Algorithm With Constant 
Bias Compensation 

The derivation of the constant bias algorithm is similar to 
the drifting bias algorithm except it is assumed that the constant bias 
(error in the directional reference ej^) is the most significant. In this 
case, each gyro mesurement is in error by Ay. Therefore, the error in y^ 
Ay and it can be computed by 
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where 


Ay 


1 



L 


''n ■ C ’'k * ''o 

k«l 

The gyro measurements are corrected by 


Y„ = Yn - ay 


and the yaw profile is computed by 


/s r— ^ 

L * ''o 

k=l 


4.6 DESCRIPTION OF THE MONTE CARLO TECHNIQUE TO DETERMINE THE MEASUREMENT 
ACCURACY OF THE DIRECTIONAL GYRO 

The Monte Carlo simulation i.sed to evaluate the directional gyro 
measurement algorithms is the same as that used for the angle cart system. 
Only the measurement systems are '.hanged. The simulation of the various 
measurement systems is described in the following sections. 
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4.6,1 Directional Gyro Without Averaging 


The directional gyro measurement algorithm without averaging 
(making only a single measurement of each conveyor section) is simulated on 
the digital computer as follows. 

The noise model defined in Section 4.5.1 is simulated on the 
digital computer assuming only a single measurement is taken on each 
conveyor segment. The equations required are: 



Xq = yo - ^0 “ ^ (initial conditions) 
2 

Gjj = normal (0, 0 |j^) 


u^ = normal (0, 1) 
Vj^ = normal (0, 1) 

a ® Cl u 
n 1 n 

b. ® Co u + Co V 
n 2 n 3 *n 


independent 


N “ *n-l * 

A/t 

J-n • *■ ^ n • 

cj^ » normal (0, oj^ ) 

'n ' *n * J'n ‘b" * * 'r 

where 

V * velocity at which the iirectional gyro traverses the 
conveyor in feet/secorul 

T = time for directional c to to traverse one conveyor 
segment. 


i 

1 

•* 


A * sampling inverval in s'conds 

2 

K = noise spectral intensity in (rad/s) /(rad/s) 
A = filter gain in second 
T ® filter time constant 



The directional gyn noise error at the n-1 conveyor 
section is added to the tri ; conveyor angle at the n**^ section to 

give the measured angle; 

I'll ' ■'n + 'n-1 

The yaw profile measurement algori'hm is 
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A KT °b ^ " ^) 2 

KA(N + 1 + •^) n + “g" n (n - 1)(3N - n - 1) + n(n - 1) + nNop 

““-272 

KA(N + 1 + ^ )N + N(N - 1)(2N - 1) + N^(N - 1)^ + 0 ^^ 


''n = ''n - ®n (V * ''n> 


The Monte Carlo simulation results for the directional gyro- 
single measurement algorithm are shown in Figures 4-24 through 4-28. In 
Figure 4-24 the mode of the RMS (Mp) is plotted versus the gyro noise 
spectral intensity K for various values of velocity V. It can be seen that 
traversing the face at 30 feet/minute (V » 0.5 ft/sec) requires a noise 
spectral intensity of 1.0 x 10“^ rad/sec to produce an Mj^ = 0.2 feet. This 
value of K corresponds to a random walk error of approximately 0.6 degrees/ 
hour. Gyros with such a random walk error are good quality and relatively expensive. 

The effect of the filter (jain and filter time constant are 
shown in Figures 4-25 and 4-26. It can i)e seen that they have little 
effect for values much larger than their nominal values. 

Figure 4-27 shows that th(‘ effect of a drifting bias has 
been removed by the optimum weighting measurement algorithm. 

The effect of the raferente error is shown in Figure 4-28. 

The results show that the amount )f reference error that can be tolerated 
depends on the gyro noise spectra intensity. The yaw advancement 
simulation was used to study the iolerable limits of these parameters as 
reported in Section 4.7.1. 


I 
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Figure 4-25. Directional Gyro - Single Measurement - Monte Carlo Results 






= as FT/SEC 
= 7.272 X 10“ 





4.6.2 Directional Gyro With Averaging 

The directional gyro measurement algorithm with averaging is 
simulated on the digital computer as follows. 

The directional gyro noise (assuming that many measurements 
of each section are taken and averaged) is simulated by the following 
equations: 


Vi = KT 
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Xq = ^0 " ^0 " ^0 ~ ^ (initidl condi Lions) 
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€|j = normal (0, ) 

2 

ep = normal (0, oj^ ) 

=» normal (0, 1) 

= normal (0, 1) n im ependent 

Wj^ = normal (0, 1) 
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The quantity is the averaged gyro noise on the n ’ conveyor section and 


is added in the simulation to the true angle to obtain the averaged 


measured angle 9^. 


A ^ . 

’n ’n n 


The yaw profile measurement algorithm is given by the following equations. 
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The Monte Carlo sinulatlon results for the directional gyro 
averaging system are shown in Fig ire 4-2'i. These results are almost 
identical to the single measureme t system. Hence taking a large number of 
measurements and averaging does n< t improve the performance over taking a 
single measurement. 

‘^•6*3 Directional Gyro willi Drifting Bias Compensation 

The directional gyre measurenent algorithm with drifting bias 
compensation is simulated using thi- same equations for the gyro error as 
for the model with no averaging. Hence, the gyro measurements are given by 


J 


I 

1 
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where is given by the equations in Section 4.6.1. The remaining 
equations are 




- OAy 



The Monte Carlo simi lation results of the directional gyro 
system with the drifting bias comj onsation are shown in Figure 4-30. This 
figure shows that for the same cotditions the mode of RMS is slightly 
larger than the directional gyro system with single measurements and 
optimum weighting. Since the drifting bias compensation algorithm is much 
simpler to instrument, this systec would be preferred if the size of the 
computer, on which the algorithm s implemented, is a factor. 
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4.6.4 Directional Gyro With Constant Bias Compensation 


The directional gyro with constant bias compensation is also 
simulated the same as the model without averaging except for the following 
equations 


.A. 





where is the gyro error which is simulated by the equations given in 
Section 4.6.1. The remaining equations are given by 


^ I ‘N 

''n “ L vk + 

k»l 


A - 1 

Ay - u - — I — 


Yn - Yn - Ay 


y;. ■^''o 

k=l 

The Monte Carlo simulation results of the constant bias 
compensation system are shown in Figures 4-31 and 4-32. Figure 4-31 
shows that for this system a directional reference error (constant bias) 
has little effect on the mode of the RMS. However, now a smaller level 
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of spectral noise Intensity can be tolerated. In addition* as shown In 
Figure 4*32* a drifting bias affects the system performance. Selection 
of the directional gyro algorithm (between the single measurement system 
and the constant bias system) must consider the relative magnitude of 
the spectral noise Intensity* the drifting bias* and the directional 
reference error. 


4.7 DESCRIPTION OF THE YAW ADVANCEMENT SYSTEM 
SIMULATION USING THE DIRECTIONAL GYRO 


This simulation is identical to that described in Section 4.4 
except that the directional gyro is used for ‘‘he angle measurements 
rather than the angle cart. The sequence of operations simulated are: 
(1) the conveyer is pushed forward, (2) the roof supports are pulled up, 
(3) the angle measurements are made, and (4) the yaw profile and con- 
veyer placement conmands are computed in preparation for pushing the 
conveyer forward again. 

Two yaw profile measurement algorithms are simulated: the single 

measurement (nonaveraging) algorithm, and the constant bias compensation 
algorithm. The Monte Carlo simul tion performance of the other two 
algorithms (averaging and driftin bias compensation) were very similar 
to the single measurement algoriti n. Therefore, they are not evaluated 
with the yaw advancement simulatiin. 

The yaw profile measurement algorithm for the single measurement 
system is simulated by 

A 

Y • ~ Y • + e * 

where is the gyro error and is a function of K, the noise spectral 
intensity, the standard deviation of the drifting bias and the 
standard deviation of the gyro reference error Oj^» 


n 



k*l 
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The yaw profile measurement algorithm for the constant bias 
compensation system is simulated by 
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^•7.1 Overall System Performance Using the Directional Gyro 


The overall system perfonnance using the directional gyro 
was determined by exercising the yaw advancement simulation. 

The conveyer initial position was a straight line. Then 
the conveyer was pushed forward (with placement errors Op) and the roof 
supports Were pulled up (with pull up errors o^). Measurements were then 
made (with noise error K) and the yaw profile computed. Next, conveyer 
placement conmands were computed in preparation for the next advancement. 
Twenty consecutive advancements were made. 

In order to determine tolerable noise errors, a placement 
error o and a pull up error were selected, and then 20 advances were 
made with K ® 0. The noise was then increased in 0.2 x 10 rad/s 
steps, each time advancing the conveyer 20 times. As the noise was 
increased, the conveyer became more crooked so that for some value of K 
the conveyer could not advance 20 times without coirmanding a negative 
placement of the conveyer. Such a command would require the conveyer to 
be pulled back which is impossible. As a result, the next lower value 
of K was deaned the maximum tolerable noise error. Three K maximums 
were determined for each placenient and pull up error , and then averaged. 
The results are plotted in Figures 4-33 and 4-34. 

Figure 4-33 shows the maximum tolerable noise error versus 
the pullup error for the single measurement (no averaging) system. It 
can be seen that gyro reference errors (a^) as large as 0.2 deg have 
little effect on the maximum tolerable spectral noise intensity. Also, 
a value of K = 1 x 10'^ rad/s for the noise can be tolerated. This 
amount of noise is equivalent to a random walk standard deviation of 0.6 
deg/h. Such a gyro Is of good quality and will be relatively expensive. 
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Figure 4-34 shows the maxinium tolerable noise error versus 
the pullup error for the constant bias compensation system. It can be 
seen that the drifting bias causes a slight decrease in the tolerable 
noise level K. A value of K = 0.5 x 10“^ rad/s can be tolerated for the 
nominal bias of = 7.27 x 10” rad/s. This value is equivalent to a 
drift rate of 1.5 deg/h and represents a good quality gyro. 


4.7.2 Recommended Measurement Algorithm for Directional Gyro 

The reconnmended measurement algorithm for the directional 
gyro is the single measurement (no averaging) system. This system uses 
the measured location of the conveyer end points and the angular 
orientation of all the conveyer sections to compute the yaw profile. 

Only a single measurement of each conveyer section is made. An optimum 
weighting technique is used to minimize the yaw profile errors. 

The single measurement system performance is comparable to 
the averaging system but does not need the large number of measurements 
for averaging. The drifting bias compensation system has a slightly 
inferior performance. The constai I bias compensation system requires 
the noise spectral intensity to bv half that of the single measurement 
system. 

Using the yaw advancement simulation, ten advances of the 
directional gyro-single measurement system were simulated. The ratio of 
the amount of coal cut compared to that which could be cut without 
errors is the efficiency of the system. Figures 4-35 and 4-36 show the 
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efficiency of the directional gyro-single measurement system for various 
system errors. 


4.8 CONCLUSIONS AND RECOMMENDATIONS 

The yaw alignment system has been developed to automatically 
advance the conveyer and the pull up of the roof supports, keeping the 
coal face relatively straight. Angular measurements of the conveyer 
sections are made so that the conveyer shape (profile) can be computed. 
Conveyer placement commands are then computed to realign the conveyer. 

Studies to evaluate this system indicate that either an angle cart 
or directional gyro can be used to measure the conveyer angles. The 
angle cart measures the relative angle between two adjacent conveyer 
sections while the directional gyro measures the angular orientation of 
a conveyer section with respect to a fixed reference. 

The studies also indicate that the accuracies required of the 
directional gyro dictate that a relatively expensive gyro is required. 

In addition to this accuracy requirement, there can be some operational 
problems. If the gyro is mounted on the shearer, and measurements are 
made as the shearer is cutting coal, any delay in traversing the face 
will result in gyro errors continuing to build up. Another measurement 
pass with the shearer must then be made. If the directional gyro is not 
mounted on the shearer to avoid delays and to speed up the measurements, 
then a separate cart is needed, increasing the system cost. 

As a result, the basic angle cart system is recommended. This 
system requires no reference angle measurements--only the angle cart 
measurements and the location of the two conveyer end points. Simula- 
tion results indicate that the measurement accuracies are an order of 
magnitude better than the require lents. 
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Roof support pull up errors of 0.1 ft can be tolerated. Estimates 
of the pull up errors indicate they are within this requirement. It is 
recommended, however, that study on an operating long wall miner be con' 
ducted to detennine more accurately the roof support pull up errors. In 
addition, it is recomniended that a dynamically accurate conveyer-roof 
support simulation be designed to study the roof support pullup errors 
since the present results are overly conservative. 
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5. 


ROLL CONTROL SYSTEM 


5.1 FUNCTIONAL DESCRIPTION 

REPRODUCIBILITY OP THE 
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The function of the roll control system is to provide an addi- 
tional degree of control for the shearer. The shearer is equipped with 
hydraulic actuators which allow it to be rotated about its longitudinal 
axis, thus tilting both of the shearing drums relative to the coal seam. 
The purpose of roll control is to provide the capability to correct for 
twists and undulations in the conveyer. These twists can occur when 
coal or other debris becomes lodged under the conveyer. 


The roll control system uses a roll sensor to measure the roll 
angle of the shearer relative to the local vertical. This signal may be 
biased to allow for operation ii coal seams which are not level. The 
hydraulic actuators are commanded by the roll error angle at sensor out- 
put, through the necessary compensation. This control system is intend- 
ed to be active, that is, to provide continuous control along the face 
during the shearing operation. The additional degree of freedom offered 
by roll control can be benefici il in maintaining the Shearer within the 
coal seam. 

5.2 ACTIVE ROLL CONTROL SYSTEM IMPLEMENTATION 


Two alternative configura'-ions were considered for the implementa- 
tion of active roll control. Tiese systems differ in the treatment of 
the roll sensor and actuator in the control loop. Figure 5-1 illus- 
trates the two loop configurations used in this study. Here, the roll 
sensor is an inclinometer. The system at the top of this figure is re- 
ferred to as the open actuator loop, and the system at the bottom is 
referred to as the closed actua or loop. This nomenclature refers to 
the treatment of the hydraulic actuator relative to the control loop. 
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Figure 5-1. Active Roll System Control Loop Iroplemen tat ions 














Since the system is intended to be active during shearing opera- 
tions, it will be subject to vibrational distrubances which will affect 
its performance. The inclinometer used as a roll sensor relies on 
gravitational acceleration to detect roll angles and therefore, any 
vibrational accelerations along the sensitive axis of this device are 
interpreted as roll signals and thus represent noise in the system. 

This problem was addressed in design and specification of the control 
loops. 


The open actuator loop system uses the incl incnnenter to sense the 
difference between the roll of the conveyor 0p ami the roll of the 
shearer s. This signal drives the actuators through the compensation 
Kj, which is discussed in Section 5.4. With this system, the inclino- 
meter is actively within the control loop and, therefore any filtering 
required to reduce noise levels will potentially affect loop stability. 

The closed actuator loop system is an alternative solution intended to 
provide a means of filtering inclinometer data without affecting the con- 
trol loop stability. A position loop was closed around the actuator as was 
done for the VCS system. The inclinometer is outside of this primary 
loop, and therefore, acts as a reference to update the actuator loop. 

These two control configurations were modeled and analyzed to determine 
their performance. 

5.3 DETAILED MATHEMATICAL MODEL FOR ROLL CONTROL SYSTEM 
5.3.1 Inclinometer Sensor Model 


The sensor modeled was a Moog model 86-121 inclinometer. 
This device uses the displacement of a sliding mass to detect inclina- 
tion relative to the local vertical. The device can be viewed as a mass 
free to slide but subject to damping and restoring forces. When the 
surface on which the mass slides is tilted with respect to horizontal » 


gravitational acceleration moves the mass from its null position. Ex- 
ternal acceleration in the direction of motion of the mass also causes 
motion. Given the following definitions: 

X * mass displacement 

M « mass 

D = damping force constant 
K » restoring force constant 
g » acceleration due to gravity 
a “ disturabnce acceleration 

the equation of the mass motion is: 

•• 

Mx “ - Dx - Kx + Mg sin '.3 + Ma 

The angle o is the inclination of the mass relative to horizontal, and 
is the quantity to be sensed. Since inclination angles are small, sin 0 
may be replaced by 9 . Figure 5-2 shows a block diagram of the solution 
for this equation. 

The Moog sensor has a natural frequency of 2.4 Hz and a 
damping ratio of 0.6. The unforced version of the above equation is: 

•• 

x + Dx + Kx=0 

with 


5-4 








2.4 Hz “ 15.08 rad/s 


(D « 


^ - 227,396 


and 


2cu)„ * S " 18.096 
n M 

The inclinometer equation can now be rewriten as: 

X = - 18.095X - 227.396X + 32. 2e + a 

As shown in Figure 5-2, the output of the sensor is the displacetnent 
divided by Mg to yield an angular equivalent. 

5.3.2 Vibration Envirorment 

The inclinometer lust be mounted on the shearer with its 
sensitive axis normal to the face, if it is to sense the desired roll 
angle. As a result of this orientation, the inclinometer will also 
sense cross-axis accelerations normal to the face resulting from 
coal shearing. The precise nature of these cross-axis disturbances 
was not known, so a bandlimited white noise process was used to simulate 

this noise. The simulated noise was a zero mean process with an ex- 
ponential correlation, i.e. white noise through a first order 
low pass linear filter. The first order filter was set to a 100 Hz 
bandwidth in these studies. 

Throughout the results which follow, the cross-axis 
disturbances will be specified as RMS g levels. The square of this, 
the signal variance, represents the total g in the 100 Hz bandwidth, or 
the area under the noise spectrum, llie spectral level of the white 


•St',? 








noise passed through the low pass filter is found by dividing the 
variance by the bandwidth. 

5.3.3 Actuator - Shearer Model 


The hydraulic acttation of the roll system is performed by 
two hydraulic cylinders operatec' in parallel. These cylinders are 22.25 
in. long when retracted, and have a ram extension of 8 in. The actua- 
tors are mounted such that a 4 m. extension gives zero roll angle rela- 
tive to the skid plane. Theref >re . within the range of the ram exten- 
sion, a roll of +5 deg is possil le. Figure 5-3 shows the shearer geo- 
metry. The distance from the siiearer pivot to actuator is 47 in. When 
the actuator is at its nominal pxtension of 4 in., ■ 29.2 deg. The 
roll of the shearer © is the ch.mge in roll relative to e^. The equa- 
tions for the shearer roll angl" in term of actuator displacement are 
given as: 


e . cos-* - (*0 * - o„ 

\ ^2 / 

This is the same relationship used in the VCS, to determine ranging arm 
angles. 


The model used f( the shearer-actuator system is the same 
as the VCS actuator system show i in Figure 3-16. For this simulation, 
however, the lengths D and W ar : replaced by lengths Oj and Dg 
respectively. The simulation w is also modified such that its range of 
roll was +5 deg. The model assumes that the pair of actuators 
can be treated as one actuator, and that the roll system has 
the same dynamic response and n >n1 inearities as the VCS system. 

Finally, the second order loop simulating arm flexibility in the VCS 
model was not used in the roll studies. 



5.4 ROLL LOOP DESIGN 


5*4.1 Open Actuator Loop Design 

The system diagram In Figure 5-1 shows a compensation K^* 
This was chosen to be a simple gain so that a roll angle sensed by the 
inclinometer would drive the actuator control valve through this gain. 
The value of Kj was selected by simulation, because of the nonlinear 
nature of the actuator system. The loop response varies as a function 
of both input amplitude and frequency, sc to measure loop performance, 
sinusoidal inputs ranging from 1.25 to 5 deg peak were used at frequen- 
cies from 0.025 to 0.2 Hz. The performance measure used was the rms of 
the loop error designated as c in Figure 5-1. Figures 5-4 through 5-6 
show the loop error versus frequency for several amplitudes. Values of 
Kj used were 275, 1375, and 2750 for the Figures 5-4 through 5-6, res- 
pectively. 


It can be seen that the low frequency error is reduced by 
Increasing Kj, but higher frequency error is essentially the same. This 
is due to the rate limited nature of the hydraulic system. At low fre- 
quencies, below 0.03 Hz, or at higher frequencies at the lower ampli- 
tudes, the actuator rates required to follow the input are less than 
maximum actuator rate, and resulting errors are determined largely by 
Kj. The larger amplitude inputs such as 5 deg, and higher frequencies above 
0.05 Hz, require actuator rates beyond its capacity and cause the actuator 
to limit at its maximum rate introducing errors. 

The maximum input frequency to be expected depends on the 
chassis velocity and amount of twist in the conveyor sections. If a 
shearer velocity of 30 ft/mi n is assumed, and letting every 5 ft con- 
veyor section be alternately twisted, the resultant roll frequency would 
be 0.05 Hz. It was, therefore, decided that the roll loop should 
perform well up to this frequency. Since it is unlikely that a 5 deg 
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error will occur regularly, an Intermediate amplitude was selected as a 
reference. Figures 5-4 through 5-6 illustrate that at a 2.5 deg peak 
Input amplitude, all values of yielded essentially a constant error up 
to 0.05Hz. A gain of 2750 was chosen because of minimum error. A roll 
error of 0.05 deg is sufficient to generate a 2.4 V signal and exceed 
the control valve input threshold with this gain. Figure 5-7 shows a 
time history of the roll angle for the system in response to a 0.05 Hz 
Input of 2.5 deg peak amplitude. 


5.4.2 Closed Actuator Loop Design 


In this system, the actuator is position commanded as in 
the VCS system. An actuator position feedback is compared with the 
command and applied to the acutator control valve through the gain Kg. 
This gain was also selected by simulation to yield a response similar 
to the open actuator loop system. The calculations necessary to de- 
termine the desired actuator displacement follow those used in the VCS. 
The present shearer roll angle 0* is computed from the current Xg as: 


e' 


cos 


-1 




2 D, 




The desired value of roll is then s' + 4> and the necessary X to achieve 

O 

this is: 

Xa - 2D^ Dg cos (.^ + e'7 -X^, 


reproducibility of TJf P 

oeiginal paoe is toot 



Actuator Loop Time History with No Noise 





Figures 5-8 through 5-10 show error versus frequency plots 
for the closed actuator loop system with equal to 85, 170, and 340 
respectively. These plots are very similar to Figures 5-4 through 5-6 
except that the low frequency error is lower and less affected by gain 
changes. Again, at 2.5 deg input amplitude, the error response is good 
to approximately 0.05 Hz. The value of K 2 chosen was 170, the same 
gain as in the VCS system. This gain yields good performance without 
excessive limit cycles. Time response of the closed actuator loop 
system is shown in Figure 5-11. It is seen that this system is slightly 
more accurate than the open loop system in a noiseless environment. 

This is also shown in the frequency response diagrams of the two systems. 

5.5. ROLL SYSTEM PERFORMANCE WITH CROSS-AXIS ACCELERATION 

5.5.1 Performance with open Actuator Loop 

System performance is shown in Table 5-1, as a function of the RMS 
cross-axis acceleration level, vhich was parameterized from 0.01 to 1.0 
g's. The data in this table shows RMS loop error for a sinusoidal input 
at a frequency of 0.05 Hz (1/20 Hz) and an amplitude of 2.5 deg peak. 
These results were also at maximum flow rate. Row 1 of Table 5-1 
indicates the response of the nominal open loop system to cross-axis 
disturbances. It can be seen that even small accelerations introduce 
errors. The input signal has ai RMS value of 1.77 deg, which would be 
the RMS error if no control was applied. Accelerations of 0.05 g produce 
errors approximately 45 percent of this value, and with noise above 0.1 
g the errors are essentially as great as with no control. 
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GAIN 340 







Figure 5-11. Closed Actuator Loop Time History with No Noise 





Row 3 of Table 5>1 shows the results of the open loop with the 
bandwidth of the Inclinometer artificially reduced by a factor of 10 to 
simulate additional filtering In the loop. Performance at lower g 
levels degrades compared with the higher bandwidth, but Improves slightly 
In the 0.1 to 0.2 g range. The principal reason for the change Is the 
basic change In system accuracy resulting from the lower bandwidth. 

Figure 5*12 Illustrates the frequency response of the reduced bandwidth 
loop. It can be seen that loop errors are higher over the entire 
frequency range. No Instabilities resulted from this additional filter- 
ing. Figures 5-13 and 5-14 show the time responses of the open loop 
system at both Inclinometer bandwidth^ with a 0.05 g noise Input level. 

System performance at a lower Input frequency of 0.025 Hz (1/40 
Hz) Is demonstrated In Table 5-2. Rows 1 and 3 of this table show the 
open loop results, which Indicate that performance Is improved over that 
obtained with the 0.05 Hz Input This is comtuensurate with the frequ- 
ency response data which also showed b()tter ptjrformance at lower frequ- 
ency. Essentially, with the slower input, the system has more capacity 
left to follow the Input and deal with the noisvs. At this Input; 
frequency. It takes a 0.5 g dis urbane i to yield errors approximately 
the same as with no control. F gures 5-15 and 5-16 show open loop 
performance at both Inclinomete' bandwidths with G'.l g noise and a 
0.025 Hz Input frequency. 
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Actuator Loop Response with 0.24 Hz Bandwidth 
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Figure 5-16. Open Actuator Loop Time History with 0.24 Hz Inclinometer Bandwidth 
0.1 g Cross-Axis Acceleration and 0.025 Hz Input 



5.5.2 Performance with Closed Actuator Loop 


Rows 2 and 4 of Tables 5-1 and 5-2 also show the results 
of the closed loop system with cross-axis disturbances. Table 5-1 shows 
the results with the 0.05 Hz input. The closed loop responses are 
essentially the same as the open loop responses at the nominal 2.4 Hz 
inclinometer bandwidth. At the 0.24 Hz inclinometer bandwidth, the 
closed loop and open loop systei.is perform the same at noise levels above 
0.2 g, but the closed loop is worse below this point. The two systems 
also behave similarly at the 0.025 Hz input frequency. While errors 
are generally smaller than with the 0.05 Hz input, the relative relation- 
ships of the loop performance is as described above. It should be 
emphasized that this loop is therefore equally as sensitive to small 
disturbances as was the open loop, Figures 5-17 through 5-20 show time 
responses of the closed actuatos- loop with various input frequencies 
and noise levels. 

5.5.3 Effects of Reduce i Flow Rate 

The hydraulic system on the shearer is shared between the 
roll and VCS control systems. Therefore, the 8 gpm nominally available 
can be shared between the 2 VCS cylinders and the roll cylinders. The 
flow rate is, therefore, variable depending on the activities of these 
systems. When flew must be shared with the VCS, there is clearly less 
flow available for roll actuation, and the maximum actuator extension 
rate is proportionately lower. This has a pronounced effect on per- 
formance. Figures 5-21 and 5-22 show that performance at higher input 










TUATOR LOOP CLOSED 



Figure 5-22. Closed Actuator Loop Performance 





frequencies has degraded with both loops when the flow rate Is cut In 
half to 4 gpm. Performance with the 2.5 deg Input degrades above 0.025 
Hz» thus the effective loop bandwidth has halved because the maximum 
flow rate halved. It is, therefore, expected that noise performance 
should degrade in the same manner as was seen by varying input frequency. 

Table 5-3 shows the noise responses of both loops and both 
inclinometer bandwidths. These results were obtained with a 0.025 Hz 
input, and they resemble the results obtained with a 0.025 Hz input at 
full flow. Results with 0.24 Hz inclinometer bandwidth is somewhat 
better than the higher flow conditions. Figures 5-23 through 5-26 show 
the time histories with reduced flow rate; Figures 5-23 and 5-24 are 
for the open loop and 5-25 and 5-26 are for the closed loop system. 

, V t'V ' 

5.6 CONCLUSIONS AND RECOMMENDATIONS ViV'-;;;; V ^ 

While the closed actuator loop is more accurate under no-noise 
conditions, both loops perform similarly with cross-axis acceleration 
disturbances. The signal input level considered in these studies was 
1.77 deg RMS. The nominal closed loop system with an input frequency of 
0.05 Hz yields an error of 43 percent of the input with a 0.05 g noise 
level. The open loop system error is 45 percent under the same conditions. 
Reducing the inclinometer bandwidth or equivalently introducing additional 
filtering has little effect on noise performance. The major conclusion 
of this study is that while the closed actuator loop is slightly better 
than the open loop, both exhibit appreciable error in the presence of 
small g noise levels. 
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Table 5-3. Perfomance with Cross-Axis Disturbances, 
0.025 Hz Input and 4 GPM Flow Rate 





TIME (SEC) 

Figure 5-23. Open Actuator Loop Time History with 2.4 Hz Inclinometer Bandwidth 
0.1 g Cross-Axis Acceleration, 4 GPM Flow Rate 
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Figure 5-25. Closed Act 
0.1 g Cros: 



Figure 5-26. Closed Actuator Loop Time History with 0.24 Hz Inclinometer Bandwidth, 
0.1 g Cross-Axis Acceleration, 4 6PM Flow Rate 



6. PRBLXMZNARY DEFINITION OF CONTROLS AND DISPLAYS 


6tl Control k Diaplay Phi losophy for Longwall - Automated guid- . 
anee and control of the longwaU system does not eliminate the need 
for controls and displays. It does, In fact, greatly Increase these 
requirements. Currently, the typical longwall C&D system consists 
of two small control panels on the shearer and a few lever controls 
located on the individual chocks. Figure 6-1 presents the control 
panel of the Joy LW300 shearer. Two of these identical panels are 
located on the shearer, one on each boom, and are used by the uhearer 
operator to control boom height, haulage speed and direction, cowl 
position, and to start and stop the motors. No displays are needed 
since the operator is located in direct proximity with the equipment 
he is controlling and can visually observe the conae^aences of his 
actions . 

With the introduction of an automated control system, it becomes 
necessary to monitor the automatically performed functions. It is 
al so necessary to provide the capability for the operator to intervene 
in the automated process and take command of the control system. 

He would then function in a semi-automatic, remote control mode to 
control the shearer and chocks from a remote location. Additionally, 
it is desirable for this remote operator to be able to interrogate 
the automatic control system with regard to sensed data and operating 
parameters in the event of malfunction or degraded performance. 

Finally, it is necessary to maintain the capability for manual control 
of the system as a back-up mode of operation and for certain activities 
more conducive to manual control, such as turn around. 

REPKOTlUCIBIl.rre OF THE 
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6.1.1 Control & Ditplav Location - The principal locations 
considered appropriate for control and display mounting in a longwall 
system are Che headgate, the tailgate, the shearer, and the chocks. 

Of these, the headgate would appear to be the most logical location 
for the main control console due to its importance in Che overall 
operation of the longwall system. The tailgate control station 

is some what redundant with the headgate panel and would be considered 
most useful if a remote control method of turn around/suinp were 
employed. A control panel on the shearer itself is mandatory. Such 
a panel is needed for manual control whenever the automatic or remote 
control systems are not usable. Also, some controls are still nec- 
essary on the chock themselves so Chat they too can be manually 
advanced whenever necessary. 

6.1.2 ConmunicaCiona & Interlocks - Since operation of the 
shearer and chock support machinery will be possible from two or 
more locations a system of built-in interlocks and means of communi- 
cation must be provided to the operators. If the shearer is being operated 
in Che manual or local mode, it must be impossible for the shearer 

to also be controlled under these cl rcumstances from the Main Control 
Console, and system interlock? must prevent this possibility. Like- 
wise, a miner working among tlie chocks must have the capability of 
Isolating one or all of Che chocks from the automatic sequencer Co 
prevent sudden and unexpected movement of a chock during a maintenance 
or inspection activity. 
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Additionally, to coordinate activities between control stations, 
conmunication will be required between all personnel involved. During 
this preliminary analysis communications were determined necessary 
between the headgate, tailgate (if implemented) and shearer. It 
is also felt that a communications package be installed on every 
5th chock. The package would consist of PA type receiver and a trans- 
mitter. This would mean that a miner working along the. row of chocks 
would never be more than 12% ft from a communication link. Additionally, 
each chock should have a control valve to isolate it, or all chocks 
from the processor to allow manual operations to be performed. 

6.2 Control and Display Definition A functional requirements 
analysis was conducted on a subsystem by subsystem basis to establish 
control and display requirements for the longwall system. The results 
of the analysis is presented in Table 6-1. By subsystem, functions 
to be commanded and functions o be monitored are noted. For each 
of these functions a preferred location is given from which to control 
or monitor the function and an indication as to whether this function 
should be controlled or monitoi ed from a dedicated or shared C&D. 

The importance of the function and the necessity of seeing the data 
in "real time" determines whether a dedicated display is warranted. 

The type of control or display which can be used to satis.fy the re- 
quirements is also given. 

It must be emphasized tha ; this functional requirements analysis 
is preliminary and it is expec :ed that as the subsystems become better 
defined the control and display requirements will also become better 
defined. It appears at this stage that the majority of those functional 
requirements listed in Table 6-1 are accurate and will remain unchanged. 
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Table 6-1. C 6 D Functional Requiremerits (Concluded) 
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However, additional requirements will most likely be generated to 
complement those already defined, 

6.2.1 Processor Accessed PiaplayB •“ Of all the display requirements, 
the processor accessed displays are probably least knovm at this 
time. The ability to interrogate the engineering proto-type system 
is very deairabl,; Once the system has been, proved out, these display 
requirementa may no longer cKisf, However , during the initial ex- 
perimental system doveiopraent n great number of system parameters 
and "nice-to-know" information can be, viannliised as being logically 
displayed, on conmand, by the processor. A very preliminary list 
of displays which may be necessary is shown in Table 6-2. This list 
ie expected to grow extensively as the processor design develops. 

6*2,2 Processor Accessed Conaiands “ As in the case of the 
processor accessed displays, the cotDrn.sndB too are ill-defined at 
tttis point in the system development. Those commands that have 
been identified include the following: the ability to adjust the 

threshold triggering level from + 2" to some other figure should 
it be demonstrated that the acceptable variation between the last 
cut follower and the CID should be mote or less stringent than the 
proposed ± 2 "; the ability to bias the roll sensor in order to 
accommodate seam conditions inclined to the local vertical; 

access to the chock command logic in order to modify the control 
nigorithm; the ability to shut down an individual chock and remove 
it from automatic control should some malfunction degrade its per- 


formance , 


Table 6-2. Preliminary List of Processor. 
Accessed Displays 


CID - DRUM POSITION 

L.K. & R.H. BOOM HYDRAULIC PRESSURE 

POWER SUPPLY OUTPUT VOLTAGE & CURRENT 

LINE VOLTAGE & CURRENT 

HYDRAULIC PRESSURE 

OUTPUT OF CONTROL LAW COMMAND 

ROLL SENSOR SIGNAL OUTPUT 

SSKITI2ED PICK SIGNAL OUTPUT 

INCLINOMETER SIGNAL OUTPUT 

LAST GUT FOLLOWER SIGNAL OUTPUT 

FACE MISALIGNMENT AT VARIOUS MEASUREMENT POINTS 

CHOCKS INDIVIDUALLY CALLED UP BY NUMBER 

- CHOCK ENABLE/DISABLE STATUS 

- RAM EXTENDED/RETRACTED 

- LOADED AGAINST ROOF /UNLOADED 
MALFUNCTION ISOLATION S TROUBLESHOOTING 
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® 3 Preliminary Control and Display Concepts 


6.3.1 Headgate Control Station Panel Layout - The preliminary 
headgate control station panel layout is shown in Figure 6-2 . The 
control panel measures 19 1/4" by 11 1/4" and provides for an ap- 
proximate 15% growth factor. It is felt, however, that should the 
growth of the C&D exceed 15% a somewhat larger panel should not 
pose any particular problem. Although actual C&D components have 
not yet been selected, representative components have been used in 
the patiel layout to insure adequate spacing both in front and behind 
the panel. Controls and displays are functionally grouped by subsystem, 
machine component, and control functions. Each of the functional 
groups will be described separately. 

Operating Mode - Two toggle swit.i.es are used. One to control 
system power on/off and one to select the operating mode. Only 
the s’eTSote or automatic mode can be selected from this panel. 

Should the "local" (or manual) mode be selected at the shearer, 
all conanand functions of the headgate control station are inter- 
locked inoperable. Which ever of the three operation modes 
are being used would be indicated by the illumination of that 
display indicator. 

b. Caution and Warning - The caution and warning enunciators 
are LED's with incorporated lens covers such as the Shelly/ 
Datatron "Brite-Eyes" . These enunciators would give positive 
indication of satisfactory operation-green, cautionary condition- 
amber, or malfunction-red. Should the methane concentration 
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Preliminary Headgate Control 








rise to 2% or should a mlfunctlon occur In an important component 
threatening personnel safety, the system would be automatically 
shut down with the enunclator indicating the reason, 
c. Shearer - The shearer controls consist of two double-acting 
awltch-lights to turn the pump motor and cutter motor on and 
off, a three position toggle to select desired course of the 
shearer, and a selector switch to select shearer speed. Although 
shearer speed Is Infinitely variable from alow to fast, a selector 
switch of approximately 8 position is used here rather than a 
potentiometer, since this would cseentially yield continuous control 
and is a simpler and more reliable Implementation particularly if 
multiplexing Js to be employed. 

<1. Vertical Control fiffitem ~ The cutter drums are raised and 
lowered by means of two center-off momentary toggles. LED digital 
readouts are used to display drum height above the skid plane 
in inches. 

Shearer cowls are chaagsd from the forward to reverse position 
by me/«no of two center-oif momentary toggles. Rear illuminated 
displays indicate the cowl position. While the cowl is in tran- 
sition between points, both lights would be off. 

Shearer location is an. LED display system to indicate where 
the shearer is with respect to the headgate and tailgate. If 
the shearer is at the headgate tf-.ilgate, that rear illuminated 
display will so indicate. As the shearer begins to move from 
the headgate the first LED will flash until the leading edge 
of the shearer arrives at the 50 foot mark. At this time the 
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first LED will burn steady and the second will flash until the 
shearer reaches the 100 foot mark. The actual distance from 
the headgate will be displayed above the row of light on a LED 
digital readout. 

The CID is deployed or stowed by use of a two position toggle. 

The desired coal thickness to be left uncut is set in with 
thumb wheels with the sensed coal thickness shown via a LED 
digital readout in inches and tf'.nth of inches. 

The Sensitized Pick displays whether It Is into coal or shale 
by the illumination of one side of the rear illuminated display. 

The LED digital readout represents a RMS level readout from 
the sensor and the thumbwheel control allows a change in out- 
put level to be comcmnded. 

Tlie cutter drum commands indicate the following. using LED digital readouts 
1) the difference between the last cut follov7er measurement 
and the CID measurement j 2) the. control law connmnd to the 
lead drum, and 3) separation between the two drum. All read- 
outs are in inches. A tv/o position toggle is used to deploy 
or stow the last cut follower. 

e. Face Alignment System - Two LED digital readouts are used. 

The first one indicates the maximum deviation of the conveyor 
in inches from a true path or measure along some pre selected 
points. The second indicates the output of the inclinometer 
in degrees. The third di.splay in this group is a rear illuminated 
indicator which will show whether the roll actuator is fully 


extended or retracted. 


Pyoceaaog - The pkoceosov i«put and call up control consists 


of a keyboard and a wmeric readout. This keyboard would require 
a code book to allow eall up of functions to be displayed and 
to input conimando. A mode nwiteh Is shown which could be used 
to simplify the coding operaj'ion. Viz, with mode "C" (for 
chocks) selected;, any of the 1/5 chocks could be called up 
simply by inpucing ita number. 

As long 4a the nuMb»'i!t of functions to be called up are 
kept to a minimum “his system VKmld provide a simple method 
of processor invx’iface. However, should the number of processor 
accessed displays and command grow significantly it might be 
desirable to utili,T.e a full computer input keyboard and an 
alpha-numeric readout. Such a concept ie shown in Figure 6- 
3. 

6.3,2 Headgate Control Station Console Configuration ~ Should 
the preliminary control panel *ihown in Figure 6-2 be large enough 
to accomodate all the necessary controls and displays, a lightweight, 
easily portable control console ahouLd be possible. Such a concept 
is shovm in figure 6-4. The console W'Ould be intrinsically safe 
and contain no power supply. Conditioned power would be supplied 
either from the shearer or from a separate power supply located near 
by at the headgate. Shovild the ccntrol panel grow significantly 
in size or should the more complex computer keyboard be required, 
a mobile chassis mount might be necessary. The concept shown in 
Figure 6-5 uses large pneumatic tires for ease of movement over 





Figure 6-3. Heaclgate Control Panel with 
Alpha-Numeric Keyboard 
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I'ough surfaces. The aonsole could be constructed as either a “sit" 
down" or "staud-up'' station. Although every effort would be made 
to keep this console iiitrinsically safe, such a mobile chassis mount 
could accommodate a permissab^e unit or could possibly house a per- 
missabla power supply package in addition to the intrinsically safe 
console. 

6.3.3 Shearer Control Pa.v^T * The preliminary concept of the 
shearer control panel is ’«hawn In Figure 6'-6. The panel consists 
of the present Joy coritiol panf>.l with the addition of the operating 
mode controls & disp'lfsys and the caution and warning enunciators. 
Further cons iderai: ion vvill be given to this control panel during 
phase IX. 

6,3.4. Tailgate Contr ol St ation ■ As a result of the phase 
I study, no firm requirement can be seen for a tailgate control 
station. However, it is not felt that the concept should be def- 
initely eliminated at this time and it is proposed that it be carried 
over into pUaB& IX as an open item. 

6.3.5 Chock Control Pane l - Each chock would maintain the 
control levels supplied by the manufacturer. In addition a switch 
or a series of switches would bo provided to disable that particular 
chock from automatie control. Possibly the chock on either side 
might also be disabled from this same position, and/or all of the 
chocka on line might be disabled as a group. 
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Figure 6-6. Shearer Coiitrol Panel - Preliminary Concept 
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6.3.6 Control and Display Component Selection - As mentioned 
earlier, components were Investigated during this phase although 
component selections have not yet been made. Some decisions were 
reached about component selection and recommendations made. These 
Include the following. 

a. Utilization of LED readouts and displays should be made 
whereever possible since they are low power self -Illuminating 
devices and meet Intrinsic safety requirements. 

b. It Is not felt necessary to hermetically seal the control 
panel, but oil tight swi tehee will be used, all connections 
potted, and sjJlcone gasket material will be used on every 
panel opening. 

c. Wlien required, rear Illuminated displays will be used In pref- 
erence to non- illuminated displays which would require external 
illumination. 

Proposed Phase I I Trade Studi es - During phase II a pre- 
liminary design of the lottgwall control and display system will be 
completed. Interfaces will be defined with the communications sub- 
system, power conditioning and distribution, and i/o assenblies . 

The functional requirements will be updated and panel layout mod- 
ified. Coiiq)onent selection will be made and various analyses con- 
ducted to Insure that system safety, reliability, maintainability, 
and human engineering requirements are met. 

Additionally a number of trade studies will be conducted during 
phase II. Tnese will include: 

a. Panel Lighting - A revlev/ of the benefits, cost, power 
requirements, and human factor aspects of various panel lighting 
schemes. Comparisons will be made of area lighting, electro- 
luminescence and exterior lighting. 
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b. Processor Keyboard - A trade study will be made of an alpha- 
numeric vs an all numeric system of processor addressing. Con- 
sideration will be given to the use of a code book and the added 
volume resulting from an alphanumeric keyboard. 

c. Power Supply Location - A trade study will be conducted 
to determine whether the C&D power supply should be located 
in the console itself, in the shearer, or as a separate unit. 

d. Tailgate Control/Monitor Station - The usefulness of a 
tailgate CSD will be reviewed and a determination of its need 
will be specified. 
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7. OVERALL LONGWALL SYSTEM REQUIREMENTS 


The following paragraphs outline the requirements of the recommended 
longwall system design by major subsystems. 

7 • 1 Vertical Control System 

Shearer Traverse Velocity - 30 ft/niin nominal 

Operational Modes 

Take all coal on top and bottom. 

Leave a prescribed amount of coal on top and take all coal 
on bottom. 

Take all coal on top and leave a prescribed amount of coal 
on bottom. 

Leave a prescribed amount of coal on top and bottom 
Manual turn around at headgate and tailgate 
Manual Operation of VCS 

Maintain present cut to within a prescribed amount of previous cut. 
Bottom and top drum interchangeable depending on cutting direction. 
Sensor complement for taking all coal on top and bottom. 

Top Drum Control 

Sensitized Pick - At least two mounted 180 deg. apart in 
cutting drum 

Last Cut Follower - Range greater than 2 in. 

EEPEODUCIBITJTY OF THE 
ORIGINAL PAGtl IS POOR 
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Bottom Drum Control 

Sensitised Pick - At least two mounted 180 deg. apart 
on cutting drum 

Last Cut Follower - Required when drums are. interchanged 
Sensor complement for leaving a prescribed amount of coal on 
top and taking all coal on bottom. 

Top Drum Control 

Active Nucleonic CID mounted 2.5 ft behind cutting drum 
Sensitised Pick At least two mounted 180 deg. apart 
on cutting drum 

Last Cut Follower - Range greater than 2 in. 

Bottom Drum Control 

Sensitized Pick - At least two mounted 180 deg. apart 
on cutting drum 

Active Nucleonic CiD - Required when drums are interchanged 
Last Cut Follower - Required when drums are interchanged. 
Sensor complement when taking all the coal on top and leaving a 
prescribed amount of coal on bottom 
Top Drum Control 

Sensitized Pick - At least two mounted 180 deg. apart on 
cutting drum 

Last Cut Follower - Range greater than 2 in. 

Present Cut Follower ~ Required when interchanging drums 


Bottom Drum Control (Drum Slaving) 

Sensitized Pick - At least two mounted 180 deg. apart 
on cutting drum 

Present Cut Follower - Range compatible with seam height 

Last Cut Follower - Required when interchanging drums 
Sensor complement when leaving a prescribed amount of coal on 
top and bottom 

Top Drum Control 

Active Nucleonic CID mounted 2.5 ft behind cutting drum 

Sensitized Pick - At least two mounted 180 deg. apart 
on cutting drum 

Last Cut Follower - Range greater than 2 in. 

Present Cut Follower - Required when cutting drums are 
interchanged 
Bottom Drum Control (Drum Slaving) 

Sensitized Pick - At least two mounted 180 deg. apart 
on cutting drum 

Present Cut Follower - Range compatible with seam height 

Last Cut Follower - Required when cutting drums are 
interchanged. 

Active Nucleonic CID - Required when cutting drums are 
interchanged. 

Active Nucleonic CID 

Mounted 2.5 ft behind cutting drum 

Twenty inch separation between source and detector. 


JOR 
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CSI crystal detector 
Source - 30 mCi CS-137 

Averaging time (Output Internal) - 0.25 sec. 

Sensitized Pick - 80 to 90 percent accurate. Discrete indica- 
tion whether coal or shale. 

Natural Radiation CID 

Used if shown to be advantageous during Phase II 

Mounted 2.5 ft behind cutting drum 

5 in. detection crystal 

Averaging time (Output Internal) - 1 sec. 

7 . 2 Roll Control Loop 

Control Loop Bandwidth - 0.05 Hz at 2.5 deg. peak sinusoidal input 
Inclinometer mounted on shearer 
Inclinometer characteristics 

Bandwidth - between .5 and 2.4 Hz 
Damping ratio - between 0.4 to 0.8 
Second Order Dynamics 

-3 

Noise within control loop bandwidth - No greater than 10 g RMS 
Roll Actuation Range »■ i 5 in. 

Roll Actuation Rate - Greater than 1 deg/sec 

7*3 Yaw Alignment System 

Automatic Face Advance 
Manual Face Advance 




7-4 


Each roof support individually coinnianded to advance conveyor. 
Angle cart system implementation with basic measurement algorithm 
(See Section 4.0). 

Angle cart integrated with shearer 

Required angle measurement accuracy - 0.03 deg one sigma 

Required resolver measurement accuracy - 2 rnTn one sigma 

Sample size of resolver measurements - 200 

Allowable ram placement error - 0.02 ft one sigma 

Allowable roof support roll up error - 0.1 ft (1.2 in) one sigma 

Controls and Di splays 

Main control and display station in headgate. 

Control and displays on shearer for manual operation 
Computer access through digital address keyboard on headgate 
control and display panel. 

No present requirement for control and displays at tailgate. 
Control and display panel inoperative when shearer is under 
manual control . 

Critical functions have dedicated display and controls in head- 
gate panel (See Section 6, Figure 6-2). 

Utilization of LED readouts whereever possible 

Use oil tight switches, all connections potted, and silicone 

gasket material used on every panel opening. 

Rear illuminated displays used as much as possible. 


